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

MODEL="boson Hubbard";
Nmax=5; 

t=1.;
U=60; 
mu="30.0 - (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=100000;
SWEEPS=2000000;
SKIP=1000;  

MEASURE[Local Density]=1;

{ T=1.; } 
