%
% Particle density calculation for 87 Rb 133 Cs mixture
%

clear all;
tic;

% Experimental parameters

lambdaX = 1064e-9;              % axial laser wavelength
lambdaT = 1064e-9;              % transverse laser wavelength

IX = [5e5;8e5;1e6;2e6;3e6;6e6]; % axial laser intensity
IT = 2e7;                       % transverse laser intensity

omegaA = 0;                     % external harmonic potential
omegaB = 0;

rhoA = 0.5;                     % total particle density
rhoB = 0.5;

% Calculation parameters

N = 20;             % number of lattice sites

maxA = 4;           % particle number limit
maxB = 4;

tmax = 32;          % truncation states limit
wlimit = 1e-3;      % truncation weight limit
slimit = 1e-10;     % sweeping convergence limit

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

data_rbcs;
dens;

toc;
