% Code revision 26 -- plot n_up and n_down; perturb Wjn every 100th iteration instead of every 20th. % Here, with a temperature of kT = 0.5, we are still 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. % Actually with this amount of polarization T_c < 0.1. >> Master(10,5,-3,.5,20,1,200); Performing initial iterations Beginning bisection of mus to find Ns 4 mu_up = 0.000000, mu_down = 0.000000, N_up = 0.848149, N_down = 0.848149, Delta_Int_error = 0.56%, N error 344.98% 5 mu_up = 0.000000, mu_down = 0.000000, N_up = 0.912775, N_down = 0.912775, Delta_Int_error = 0.56%, N error 330.30% 6 mu_up = 10.500000, mu_down = 10.500000, N_up = 10.685152, N_down = 10.685152, Delta_Int_error = 0.58%, N error 330.30% 7 mu_up = 10.500000, mu_down = 10.500000, N_up = 13.007988, N_down = 13.007988, Delta_Int_error = 0.33%, N error 138.53% 8 mu_up = 5.250000, mu_down = 5.250000, N_up = 9.291453, N_down = 9.291453, Delta_Int_error = 0.25%, N error 138.53% 9 mu_up = 5.250000, mu_down = 5.250000, N_up = 8.295191, N_down = 8.295191, Delta_Int_error = 0.38%, N error 77.12% 10 mu_up = 7.875000, mu_down = 2.625000, N_up = 9.916591, N_down = 5.968198, Delta_Int_error = 0.41%, N error 77.12% 11 mu_up = 7.875000, mu_down = 2.625000, N_up = 10.435446, N_down = 5.171485, Delta_Int_error = 0.24%, N error 66.76% 12 mu_up = 7.875000, mu_down = 3.150000, N_up = 10.631536, N_down = 5.243370, Delta_Int_error = 0.54%, N error 66.76% 13 mu_up = 7.875000, mu_down = 3.150000, N_up = 10.716385, N_down = 5.241715, Delta_Int_error = 0.47%, N error 80.25% 14 mu_up = 7.875000, mu_down = 3.780000, N_up = 10.756261, N_down = 5.685289, Delta_Int_error = 0.40%, N error 80.25% 15 mu_up = 7.875000, mu_down = 3.780000, N_up = 10.776214, N_down = 5.849331, Delta_Int_error = 0.26%, N error 98.26% 16 mu_up = 7.875000, mu_down = 4.536000, N_up = 10.786450, N_down = 6.461549, Delta_Int_error = 0.21%, N error 98.26% 17 mu_up = 7.875000, mu_down = 4.536000, N_up = 10.791674, N_down = 6.693993, Delta_Int_error = 0.11%, N error 55.85% 18 mu_up = 6.562500, mu_down = 2.268000, N_up = 9.798725, N_down = 5.138818, Delta_Int_error = 0.07%, N error 55.85% 19 mu_up = 6.562500, mu_down = 2.268000, N_up = 9.492768, N_down = 4.584177, Delta_Int_error = 0.13%, N error 56.43% 20 mu_up = 6.562500, mu_down = 2.721600, N_up = 9.374804, N_down = 4.681489, Delta_Int_error = 0.18%, N error 56.43% 21 mu_up = 6.562500, mu_down = 2.721600, N_up = 9.321143, N_down = 4.703239, Delta_Int_error = 0.13%, N error 68.01% 22 mu_up = 6.562500, mu_down = 3.265920, N_up = 9.294258, N_down = 5.081278, Delta_Int_error = 0.11%, N error 68.01% 23 mu_up = 6.562500, mu_down = 3.265920, N_up = 9.280071, N_down = 5.217515, Delta_Int_error = 0.04%, N error 83.33% 24 mu_up = 6.562500, mu_down = 3.919104, N_up = 9.272412, N_down = 5.754248, Delta_Int_error = 0.01%, N error 83.33% 25 mu_up = 6.562500, mu_down = 3.919104, N_up = 9.267937, N_down = 5.963560, Delta_Int_error = 0.06%, N error 102.29% 26 mu_up = 6.562500, mu_down = 4.702925, N_up = 9.265330, N_down = 6.629988, Delta_Int_error = 0.09%, N error 102.29% 27 mu_up = 6.562500, mu_down = 4.702925, N_up = 9.263219, N_down = 6.882074, Delta_Int_error = 0.18%, N error 58.09% 28 mu_up = 7.218750, mu_down = 2.351462, N_up = 9.758481, N_down = 5.265679, Delta_Int_error = 0.21%, N error 58.09% 29 mu_up = 7.218750, mu_down = 2.351462, N_up = 9.918000, N_down = 4.694841, Delta_Int_error = 0.14%, N error 58.73% 30 mu_up = 7.218750, mu_down = 2.821755, N_up = 9.979157, N_down = 4.800776, Delta_Int_error = 0.27%, N error 58.73% 31 mu_up = 7.218750, mu_down = 2.821755, N_up = 10.006193, N_down = 4.825431, Delta_Int_error = 0.24%, N error 70.83% 32 mu_up = 7.218750, mu_down = 3.386106, N_up = 10.019315, N_down = 5.213669, Delta_Int_error = 0.21%, N error 70.83% 33 mu_up = 7.218750, mu_down = 3.386106, N_up = 10.026017, N_down = 5.356090, Delta_Int_error = 0.14%, N error 86.82% 34 mu_up = 7.218750, mu_down = 4.063327, N_up = 10.029587, N_down = 5.918052, Delta_Int_error = 0.11%, N error 86.82% 35 mu_up = 7.218750, mu_down = 4.063327, N_up = 10.031452, N_down = 6.128882, Delta_Int_error = 0.03%, N error 106.57% 36 mu_up = 7.218750, mu_down = 4.875992, N_up = 10.032491, N_down = 6.825654, Delta_Int_error = 0.01%, N error 106.57% 37 mu_up = 7.218750, mu_down = 4.875992, N_up = 10.032865, N_down = 7.080417, Delta_Int_error = 0.07%, N error 60.56% 38 mu_up = 6.890625, mu_down = 2.437996, N_up = 9.794162, N_down = 5.397545, Delta_Int_error = 0.11%, N error 60.56% 39 mu_up = 6.890625, mu_down = 2.437996, N_up = 9.719982, N_down = 4.806610, Delta_Int_error = 0.09%, N error 61.15% 40 mu_up = 6.956250, mu_down = 2.925595, N_up = 9.739658, N_down = 4.919296, Delta_Int_error = 0.17%, N error 61.15% 41 mu_up = 6.956250, mu_down = 2.925595, N_up = 9.741721, N_down = 4.945138, Delta_Int_error = 0.14%, N error 73.78% 42 mu_up = 6.956250, mu_down = 3.510715, N_up = 9.741109, N_down = 5.350770, Delta_Int_error = 0.12%, N error 73.78% 43 mu_up = 6.956250, mu_down = 3.510715, N_up = 9.740384, N_down = 5.504731, Delta_Int_error = 0.06%, N error 90.45% 44 mu_up = 6.956250, mu_down = 4.212857, N_up = 9.739964, N_down = 6.086558, Delta_Int_error = 0.04%, N error 90.45% 45 mu_up = 6.956250, mu_down = 4.212857, N_up = 9.739669, N_down = 6.301357, Delta_Int_error = 0.03%, N error 111.01% 46 mu_up = 6.956250, mu_down = 5.055429, N_up = 9.739508, N_down = 7.026717, Delta_Int_error = 0.05%, N error 111.01% 47 mu_up = 6.956250, mu_down = 5.055429, N_up = 9.739279, N_down = 7.285843, Delta_Int_error = 0.14%, N error 63.08% 48 mu_up = 7.153125, mu_down = 2.527714, N_up = 9.884262, N_down = 5.536182, Delta_Int_error = 0.17%, N error 63.08% 49 mu_up = 7.153125, mu_down = 2.527714, N_up = 9.929074, N_down = 4.923741, Delta_Int_error = 0.09%, N error 63.71% 50 mu_up = 7.153125, mu_down = 3.033257, N_up = 9.946215, N_down = 5.041602, Delta_Int_error = 0.19%, N error 63.71% 51 mu_up = 7.153125, mu_down = 3.033257, N_up = 9.953906, N_down = 5.067038, Delta_Int_error = 0.16%, N error 76.88% 52 mu_up = 7.153125, mu_down = 3.639909, N_up = 9.957713, N_down = 5.496851, Delta_Int_error = 0.15%, N error 76.88% 53 mu_up = 7.153125, mu_down = 3.639909, N_up = 9.959671, N_down = 5.663118, Delta_Int_error = 0.08%, N error 94.23% 54 mu_up = 7.153125, mu_down = 4.367891, N_up = 9.960744, N_down = 6.261573, Delta_Int_error = 0.05%, N error 94.23% 55 mu_up = 7.153125, mu_down = 4.367891, N_up = 9.961289, N_down = 6.485409, Delta_Int_error = 0.02%, N error 53.40% 56 mu_up = 7.251562, mu_down = 2.183945, N_up = 10.032797, N_down = 5.002301, Delta_Int_error = 0.04%, N error 53.40% 57 mu_up = 7.251562, mu_down = 2.183945, N_up = 10.054457, N_down = 4.465379, Delta_Int_error = 0.18%, N error 54.14% 58 mu_up = 7.251562, mu_down = 2.620734, N_up = 10.062661, N_down = 4.552809, Delta_Int_error = 0.25%, N error 54.14% 59 mu_up = 7.251562, mu_down = 2.620734, N_up = 10.066331, N_down = 4.570062, Delta_Int_error = 0.22%, N error 65.18% 60 mu_up = 7.251562, mu_down = 3.144881, N_up = 10.068143, N_down = 4.941175, Delta_Int_error = 0.22%, N error 65.18% 61 mu_up = 7.251562, mu_down = 3.144881, N_up = 10.069077, N_down = 5.077439, Delta_Int_error = 0.16%, N error 79.85% 62 mu_up = 7.251562, mu_down = 3.773857, N_up = 10.069576, N_down = 5.585838, Delta_Int_error = 0.15%, N error 79.85% 63 mu_up = 7.251562, mu_down = 3.773857, N_up = 10.069839, N_down = 5.790125, Delta_Int_error = 0.08%, N error 98.09% 64 mu_up = 7.251562, mu_down = 4.528629, N_up = 10.069989, N_down = 6.428736, Delta_Int_error = 0.06%, N error 98.09% 65 mu_up = 7.251562, mu_down = 4.528629, N_up = 10.070066, N_down = 6.672515, Delta_Int_error = 0.01%, N error 55.64% 66 mu_up = 7.202344, mu_down = 2.264314, N_up = 10.034391, N_down = 5.124388, Delta_Int_error = 0.03%, N error 55.64% 67 mu_up = 7.202344, mu_down = 2.264314, N_up = 10.023694, N_down = 4.573642, Delta_Int_error = 0.14%, N error 56.34% 68 mu_up = 7.202344, mu_down = 2.717177, N_up = 10.019696, N_down = 4.672119, Delta_Int_error = 0.17%, N error 56.34% 69 mu_up = 7.202344, mu_down = 2.717177, N_up = 10.017925, N_down = 4.694659, Delta_Int_error = 0.16%, N error 67.90% 70 mu_up = 7.202344, mu_down = 3.260613, N_up = 10.017060, N_down = 5.073256, Delta_Int_error = 0.17%, N error 67.90% 71 mu_up = 7.202344, mu_down = 3.260613, N_up = 10.016615, N_down = 5.209786, Delta_Int_error = 0.12%, N error 83.20% 72 mu_up = 7.202344, mu_down = 3.912735, N_up = 10.016383, N_down = 5.745581, Delta_Int_error = 0.11%, N error 83.20% 73 mu_up = 7.202344, mu_down = 3.912735, N_up = 10.016260, N_down = 5.954736, Delta_Int_error = 0.05%, N error 102.14% 74 mu_up = 7.202344, mu_down = 4.695283, N_up = 10.016196, N_down = 6.619837, Delta_Int_error = 0.03%, N error 102.14% 75 mu_up = 7.202344, mu_down = 4.695283, N_up = 10.016161, N_down = 6.871389, Delta_Int_error = 0.03%, N error 58.03% 76 mu_up = 7.177734, mu_down = 2.347641, N_up = 9.998318, N_down = 5.255696, Delta_Int_error = 0.05%, N error 58.03% 77 mu_up = 7.177734, mu_down = 2.347641, N_up = 9.992951, N_down = 4.686163, Delta_Int_error = 0.11%, N error 58.63% 78 mu_up = 7.177734, mu_down = 2.817170, N_up = 9.990938, N_down = 4.792753, Delta_Int_error = 0.13%, N error 58.63% 79 mu_up = 7.177734, mu_down = 2.817170, N_up = 9.990043, N_down = 4.818106, Delta_Int_error = 0.12%, N error 70.72% 80 mu_up = 7.177734, mu_down = 3.380603, N_up = 9.989605, N_down = 5.206398, Delta_Int_error = 0.13%, N error 70.72% 81 mu_up = 7.177734, mu_down = 3.380603, N_up = 9.989380, N_down = 5.348721, Delta_Int_error = 0.09%, N error 86.67% 82 mu_up = 7.177734, mu_down = 4.056724, N_up = 9.989262, N_down = 5.909968, Delta_Int_error = 0.08%, N error 86.67% 83 mu_up = 7.177734, mu_down = 4.056724, N_up = 9.989199, N_down = 6.120860, Delta_Int_error = 0.03%, N error 106.39% 84 mu_up = 7.177734, mu_down = 4.868069, N_up = 9.989166, N_down = 6.816428, Delta_Int_error = 0.02%, N error 106.39% 85 mu_up = 7.177734, mu_down = 4.868069, N_up = 9.989148, N_down = 7.071024, Delta_Int_error = 0.04%, N error 60.45% 86 mu_up = 7.190039, mu_down = 2.434034, N_up = 9.998041, N_down = 5.390582, Delta_Int_error = 0.07%, N error 60.45% 87 mu_up = 7.190039, mu_down = 2.434034, N_up = 10.000716, N_down = 4.800567, Delta_Int_error = 0.09%, N error 61.05% 88 mu_up = 7.190039, mu_down = 2.920841, N_up = 10.001721, N_down = 4.913173, Delta_Int_error = 0.12%, N error 61.05% 89 mu_up = 7.190039, mu_down = 2.920841, N_up = 10.002167, N_down = 4.939165, Delta_Int_error = 0.11%, N error 73.65% 90 mu_up = 7.190039, mu_down = 3.505010, N_up = 10.002387, N_down = 5.344004, Delta_Int_error = 0.12%, N error 73.65% 91 mu_up = 7.190039, mu_down = 3.505010, N_up = 10.002500, N_down = 5.497457, Delta_Int_error = 0.08%, N error 90.29% 92 mu_up = 7.190039, mu_down = 4.206012, N_up = 10.002560, N_down = 6.078667, Delta_Int_error = 0.07%, N error 90.29% 93 mu_up = 7.190039, mu_down = 4.206012, N_up = 10.002591, N_down = 6.293220, Delta_Int_error = 0.01%, N error 110.81% 94 mu_up = 7.190039, mu_down = 5.047214, N_up = 10.002608, N_down = 7.017466, Delta_Int_error = 0.00%, N error 110.81% 95 mu_up = 7.190039, mu_down = 5.047214, N_up = 10.002617, N_down = 7.276205, Delta_Int_error = 0.06%, N error 62.97% 96 mu_up = 7.183887, mu_down = 2.523607, N_up = 9.998166, N_down = 5.529209, Delta_Int_error = 0.10%, N error 62.97% 97 mu_up = 7.183887, mu_down = 2.523607, N_up = 9.996829, N_down = 4.917728, Delta_Int_error = 0.08%, N error 63.60% 98 mu_up = 7.183887, mu_down = 3.028328, N_up = 9.996328, N_down = 5.035481, Delta_Int_error = 0.12%, N error 63.60% 99 mu_up = 7.183887, mu_down = 3.028328, N_up = 9.996106, N_down = 5.061052, Delta_Int_error = 0.10%, N error 76.74% 100 mu_up = 7.183887, mu_down = 3.633994, N_up = 9.995997, N_down = 5.489770, Delta_Int_error = 0.11%, N error 76.74%