LATTICE="inhomogeneous simple cubic lattice";
L=21; 

MODEL="boson Hubbard";
Nmax=5; 

t=1.;
U=8.11; 
mu="4.05 - (0.2*(x-(L-1)/2.)*(x-(L-1)/2.) + 0.2*(y-(L-1)/2.)*(y-(L-1)/2.) + 0.2*(z-(L-1)/2.)*(z-(L-1)/2.))";


THERMALIZATION=2000;
SWEEPS=20000;
SKIP=10;  

MEASURE[Local Density]=1;

{ T=1.; } 
