% Code revision 26 -- Slightly wider plots; perturb Wjn every 20th iteration instead of every 10th. % Here, with a temperature of kT = 1, we are just above the critical temperature and Delta is small. The atoms form a solid lump in the centre of the trap with a halo of the excess polarization. >> Master(10,5,-3,1,20,1,100); Performing initial iterations Beginning bisection of mus to find Ns 4 mu_up = 0.000000, mu_down = 0.000000, N_up = 1.189461, N_down = 1.189461, Delta_Int_error = 0.32%, N error 339.44% 5 mu_up = 0.000000, mu_down = 0.000000, N_up = 1.208171, N_down = 1.208171, Delta_Int_error = 0.32%, N error 319.82% 6 mu_up = 10.500000, mu_down = 10.500000, N_up = 10.754984, N_down = 10.754984, Delta_Int_error = 0.39%, N error 319.82% 7 mu_up = 10.500000, mu_down = 10.500000, N_up = 13.046005, N_down = 13.046005, Delta_Int_error = 0.18%, N error 139.28% 8 mu_up = 5.250000, mu_down = 5.250000, N_up = 9.257148, N_down = 9.257148, Delta_Int_error = 0.06%, N error 139.28% 9 mu_up = 5.250000, mu_down = 5.250000, N_up = 8.246607, N_down = 8.246607, Delta_Int_error = 0.23%, N error 76.94% 10 mu_up = 7.875000, mu_down = 2.625000, N_up = 10.043655, N_down = 5.819399, Delta_Int_error = 0.26%, N error 76.94% 11 mu_up = 7.875000, mu_down = 2.625000, N_up = 10.549211, N_down = 5.043901, Delta_Int_error = 0.14%, N error 68.59% 12 mu_up = 7.875000, mu_down = 3.150000, N_up = 10.707017, N_down = 5.174351, Delta_Int_error = 0.40%, N error 68.59% 13 mu_up = 7.875000, mu_down = 3.150000, N_up = 10.762134, N_down = 5.210681, Delta_Int_error = 0.36%, N error 82.61% 14 mu_up = 7.875000, mu_down = 3.780000, N_up = 10.783346, N_down = 5.708681, Delta_Int_error = 0.32%, N error 82.61% 15 mu_up = 7.875000, mu_down = 3.780000, N_up = 10.792091, N_down = 5.880544, Delta_Int_error = 0.21%, N error 100.48% 16 mu_up = 7.875000, mu_down = 4.536000, N_up = 10.795858, N_down = 6.535807, Delta_Int_error = 0.18%, N error 100.48% 17 mu_up = 7.875000, mu_down = 4.536000, N_up = 10.797518, N_down = 6.758173, Delta_Int_error = 0.08%, N error 57.73% 18 mu_up = 6.562500, mu_down = 2.268000, N_up = 9.756283, N_down = 5.087199, Delta_Int_error = 0.06%, N error 57.73% 19 mu_up = 6.562500, mu_down = 2.268000, N_up = 9.451766, N_down = 4.513643, Delta_Int_error = 0.11%, N error 58.42% 20 mu_up = 6.562500, mu_down = 2.721600, N_up = 9.349100, N_down = 4.641443, Delta_Int_error = 0.16%, N error 58.42% 21 mu_up = 6.562500, mu_down = 2.721600, N_up = 9.310433, N_down = 4.679420, Delta_Int_error = 0.12%, N error 70.57% 22 mu_up = 6.562500, mu_down = 3.265920, N_up = 9.295520, N_down = 5.110860, Delta_Int_error = 0.37%, N error 70.57% 23 mu_up = 6.562500, mu_down = 3.265920, N_up = 9.288187, N_down = 5.264195, Delta_Int_error = 0.31%, N error 85.81% 24 mu_up = 6.562500, mu_down = 3.919104, N_up = 9.284927, N_down = 5.830530, Delta_Int_error = 0.11%, N error 85.81% 25 mu_up = 6.562500, mu_down = 3.919104, N_up = 9.283447, N_down = 6.029021, Delta_Int_error = 0.01%, N error 104.48% 26 mu_up = 6.562500, mu_down = 4.702925, N_up = 9.282769, N_down = 6.718554, Delta_Int_error = 0.05%, N error 104.48% 27 mu_up = 6.562500, mu_down = 4.702925, N_up = 9.282452, N_down = 6.951442, Delta_Int_error = 0.12%, N error 60.04% 28 mu_up = 7.218750, mu_down = 2.351462, N_up = 9.802256, N_down = 5.215701, Delta_Int_error = 0.14%, N error 60.04% 29 mu_up = 7.218750, mu_down = 2.351462, N_up = 9.956184, N_down = 4.625053, Delta_Int_error = 0.17%, N error 60.73% 30 mu_up = 7.218750, mu_down = 2.821755, N_up = 10.008347, N_down = 4.761454, Delta_Int_error = 0.30%, N error 60.73% 31 mu_up = 7.218750, mu_down = 2.821755, N_up = 10.028103, N_down = 4.801981, Delta_Int_error = 0.25%, N error 73.36% 32 mu_up = 7.218750, mu_down = 3.386106, N_up = 10.036207, N_down = 5.249621, Delta_Int_error = 0.23%, N error 73.36% 33 mu_up = 7.218750, mu_down = 3.386106, N_up = 10.039701, N_down = 5.408091, Delta_Int_error = 0.14%, N error 89.22% 34 mu_up = 7.218750, mu_down = 4.063327, N_up = 10.041253, N_down = 5.995585, Delta_Int_error = 0.11%, N error 89.22% 35 mu_up = 7.218750, mu_down = 4.063327, N_up = 10.041953, N_down = 6.199929, Delta_Int_error = 0.03%, N error 108.65% 36 mu_up = 7.218750, mu_down = 4.875992, N_up = 10.042272, N_down = 6.914856, Delta_Int_error = 0.01%, N error 108.65% 37 mu_up = 7.218750, mu_down = 4.875992, N_up = 10.042417, N_down = 7.154310, Delta_Int_error = 0.06%, N error 62.49% 38 mu_up = 6.890625, mu_down = 2.437996, N_up = 9.782922, N_down = 5.349562, Delta_Int_error = 0.08%, N error 62.49% 39 mu_up = 6.890625, mu_down = 2.437996, N_up = 9.706366, N_down = 4.740579, Delta_Int_error = 0.10%, N error 63.14% 40 mu_up = 6.956250, mu_down = 2.925595, N_up = 9.732145, N_down = 4.885745, Delta_Int_error = 0.18%, N error 63.14% 41 mu_up = 6.956250, mu_down = 2.925595, N_up = 9.737544, N_down = 4.928801, Delta_Int_error = 0.15%, N error 76.27% 42 mu_up = 6.956250, mu_down = 3.510715, N_up = 9.737870, N_down = 5.395430, Delta_Int_error = 0.57%, N error 76.27% 43 mu_up = 6.956250, mu_down = 3.510715, N_up = 9.738747, N_down = 5.557768, Delta_Int_error = 0.12%, N error 92.76% 44 mu_up = 6.956250, mu_down = 4.212857, N_up = 9.738908, N_down = 6.166680, Delta_Int_error = 0.05%, N error 92.76% 45 mu_up = 6.956250, mu_down = 4.212857, N_up = 9.738924, N_down = 6.376843, Delta_Int_error = 0.04%, N error 112.98% 46 mu_up = 6.956250, mu_down = 5.055429, N_up = 9.738917, N_down = 7.118036, Delta_Int_error = 0.05%, N error 112.98% 47 mu_up = 6.956250, mu_down = 5.055429, N_up = 9.738910, N_down = 7.364189, Delta_Int_error = 0.00%, N error 65.03% 48 mu_up = 7.153125, mu_down = 2.527714, N_up = 9.894753, N_down = 5.487640, Delta_Int_error = 0.01%, N error 65.03% 49 mu_up = 7.153125, mu_down = 2.527714, N_up = 9.940712, N_down = 4.859849, Delta_Int_error = 0.08%, N error 65.64% 50 mu_up = 7.153125, mu_down = 3.033257, N_up = 9.956344, N_down = 5.014254, Delta_Int_error = 0.07%, N error 65.64% 51 mu_up = 7.153125, mu_down = 3.033257, N_up = 9.962302, N_down = 5.059939, Delta_Int_error = 0.07%, N error 79.30% 52 mu_up = 7.153125, mu_down = 3.639909, N_up = 9.964760, N_down = 5.543080, Delta_Int_error = 0.08%, N error 79.30% 53 mu_up = 7.153125, mu_down = 3.639909, N_up = 9.965825, N_down = 5.711629, Delta_Int_error = 0.04%, N error 96.46% 54 mu_up = 7.153125, mu_down = 4.367891, N_up = 9.966300, N_down = 6.343469, Delta_Int_error = 0.04%, N error 96.46% 55 mu_up = 7.153125, mu_down = 4.367891, N_up = 9.966515, N_down = 6.559800, Delta_Int_error = 0.02%, N error 55.37% 56 mu_up = 7.251562, mu_down = 2.183945, N_up = 10.044588, N_down = 4.955545, Delta_Int_error = 0.03%, N error 55.37% 57 mu_up = 7.251562, mu_down = 2.183945, N_up = 10.067516, N_down = 4.400263, Delta_Int_error = 0.16%, N error 56.09% 58 mu_up = 7.251562, mu_down = 2.620734, N_up = 10.075302, N_down = 4.519892, Delta_Int_error = 0.23%, N error 56.09% 59 mu_up = 7.251562, mu_down = 2.620734, N_up = 10.078269, N_down = 4.555472, Delta_Int_error = 0.21%, N error 67.76% 60 mu_up = 7.231875, mu_down = 3.144881, N_up = 10.063898, N_down = 4.969444, Delta_Int_error = 0.20%, N error 67.76%