% Code revision 25 -- Much improved plotting (two graphs in same window; pause after every plot). Avoid the pause by holding down an arrow key.

% Here, with a temperature of kT = 10, we are above the critical temperature and Delta -> 0. The atoms form a solid lump in the centre of the trap.

>> Master(10,8,-3,10,20,1,100);
Performing initial iterations
Beginning bisection of mus to find Ns
mu_up =   0.000000, mu_down =   0.000000, N_up =   6.504206, N_down =   6.504206, Delta_Int_error =  3.05%, N error 162.04%
mu_up =   0.000000, mu_down =   0.000000, N_up =   6.504357, N_down =   6.504357, Delta_Int_error =  2.97%, N error 103.57%
mu_up =  10.500000, mu_down =  10.500000, N_up =  11.313020, N_down =  11.313020, Delta_Int_error =  2.77%, N error 103.57%
mu_up =  10.500000, mu_down =  10.500000, N_up =  11.883251, N_down =  11.883251, Delta_Int_error =  2.75%, N error 53.13%
mu_up =   5.250000, mu_down =   5.250000, N_up =   9.477299, N_down =   9.477299, Delta_Int_error =  2.76%, N error 53.13%
mu_up =   5.250000, mu_down =   5.250000, N_up =   9.191218, N_down =   9.191218, Delta_Int_error =  2.67%, N error 26.51%
mu_up =   7.875000, mu_down =   2.625000, N_up =  10.394765, N_down =   7.949800, Delta_Int_error =  2.62%, N error 26.51%
mu_up =   7.875000, mu_down =   2.625000, N_up =  10.539324, N_down =   7.803847, Delta_Int_error =  2.67%, N error 25.96%
mu_up =   7.875000, mu_down =   3.150000, N_up =  10.557003, N_down =   8.021367, Delta_Int_error =  1.00%, N error 25.96%
mu_up =   7.875000, mu_down =   3.150000, N_up =  10.559434, N_down =   8.047194, Delta_Int_error =  7.91%, N error 31.46%
mu_up =   7.875000, mu_down =   3.780000, N_up =  10.559756, N_down =   8.336746, Delta_Int_error =  0.40%, N error 31.46%
mu_up =   7.875000, mu_down =   3.780000, N_up =  10.559798, N_down =   8.371389, Delta_Int_error =  1.84%, N error 18.42%
mu_up =   6.562500, mu_down =   1.890000, N_up =   9.938260, N_down =   7.526268, Delta_Int_error =  2.29%, N error 18.42%
mu_up =   6.562500, mu_down =   1.890000, N_up =   9.863635, N_down =   7.428180, Delta_Int_error =  2.47%, N error 18.38%
mu_up =   6.562500, mu_down =   2.268000, N_up =   9.854325, N_down =   7.582613, Delta_Int_error =  2.54%, N error 18.38%
mu_up =   6.562500, mu_down =   2.268000, N_up =   9.853138, N_down =   7.600589, Delta_Int_error =  2.55%, N error 22.25%
mu_up =   6.562500, mu_down =   2.721600, N_up =   9.852985, N_down =   7.804625, Delta_Int_error =  2.55%, N error 22.25%
mu_up =   6.562500, mu_down =   2.721600, N_up =   9.852965, N_down =   7.828722, Delta_Int_error =  2.53%, N error 26.96%
mu_up =   6.562500, mu_down =   3.265920, N_up =   9.854086, N_down =   8.076599, Delta_Int_error =  1.00%, N error 26.96%
mu_up =   6.562500, mu_down =   3.265920, N_up =   9.853092, N_down =   8.106189, Delta_Int_error =  4.19%, N error 32.68%
mu_up =   6.562500, mu_down =   3.919104, N_up =   9.852978, N_down =   8.407505, Delta_Int_error =  3.68%, N error 32.68%
mu_up =   6.562500, mu_down =   3.919104, N_up =   9.852964, N_down =   8.443608, Delta_Int_error =  3.28%, N error 19.13%
mu_up =   7.218750, mu_down =   1.959552, N_up =  10.163568, N_down =   7.565365, Delta_Int_error =  3.02%, N error 19.13%
mu_up =   7.218750, mu_down =   1.959552, N_up =  10.200929, N_down =   7.463318, Delta_Int_error =  2.92%, N error 19.09%
mu_up =   7.218750, mu_down =   2.351462, N_up =  10.205579, N_down =   7.623723, Delta_Int_error =  2.83%, N error 19.09%
mu_up =   7.218750, mu_down =   2.351462, N_up =  10.206169, N_down =   7.642419, Delta_Int_error =  2.75%, N error 23.12%
mu_up =   7.218750, mu_down =   2.821755, N_up =  10.206244, N_down =   7.854434, Delta_Int_error =  2.71%, N error 23.12%
mu_up =   7.218750, mu_down =   2.821755, N_up =  10.206254, N_down =   7.879510, Delta_Int_error =  2.66%, N error 28.01%
mu_up =   7.218750, mu_down =   3.386106, N_up =  10.207475, N_down =   8.136819, Delta_Int_error =  1.00%, N error 28.01%
mu_up =   7.218750, mu_down =   3.386106, N_up =  10.206410, N_down =   8.167902, Delta_Int_error =  0.50%, N error 33.95%
mu_up =   7.218750, mu_down =   4.063327, N_up =  10.206275, N_down =   8.481116, Delta_Int_error =  1.76%, N error 33.95%
mu_up =   7.218750, mu_down =   4.063327, N_up =  10.206258, N_down =   8.518692, Delta_Int_error =  2.22%, N error 19.87%
mu_up =   6.890625, mu_down =   2.031664, N_up =  10.050871, N_down =   7.605998, Delta_Int_error =  2.40%, N error 19.87%
mu_up =   6.890625, mu_down =   2.031664, N_up =  10.032194, N_down =   7.499826, Delta_Int_error =  2.52%, N error 19.82%
mu_up =   6.890625, mu_down =   2.437996, N_up =  10.029866, N_down =   7.666441, Delta_Int_error =  2.56%, N error 19.82%
mu_up =   6.890625, mu_down =   2.437996, N_up =  10.029570, N_down =   7.685889, Delta_Int_error =  2.57%, N error 24.01%
mu_up =   6.890625, mu_down =   2.925595, N_up =  10.029532, N_down =   7.906203, Delta_Int_error =  2.58%, N error 24.01%
mu_up =   6.890625, mu_down =   2.925595, N_up =  10.029527, N_down =   7.932299, Delta_Int_error =  2.57%, N error 29.10%
mu_up =   6.890625, mu_down =   3.510715, N_up =  10.030247, N_down =   8.200753, Delta_Int_error =  1.00%, N error 29.10%