% Code revision 21 - Code probably correct -- initial guess for mu_min now very negative to allow for strong interactions reducing energy.
% Code probably needs more iterations between mu bisections

>> Master(8,0,-3,0.0000001,25,1,100);
Performing initial iterations
Beginning bisection of mus to find Ns
mu_up =   0.000000, mu_down =   0.000000, N_up =   0.036600, N_down =   0.036600, Delta_Int_error =   1.000000%
mu_up =  13.000000, mu_down = -13.000000, N_up =  13.000515, N_down =   0.000515, Delta_Int_error =   0.893105%
mu_up =   6.500000, mu_down = -19.500000, N_up =  10.000004, N_down =   0.000004, Delta_Int_error =   8.103932%
mu_up =   3.250000, mu_down = -22.750000, N_up =   7.000000, N_down =   0.000000, Delta_Int_error =   7.628391%
mu_up =   4.875000, mu_down = -24.375000, N_up =   7.000000, N_down =   0.000000, Delta_Int_error =   5.772693%
mu_up =   5.687500, mu_down = -25.187500, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   5.715128%
mu_up =   5.281250, mu_down = -25.593750, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   6.233878%
mu_up =   5.078125, mu_down = -25.796875, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   5.448684%
mu_up =   4.976562, mu_down = -25.898438, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   4.766694%
mu_up =   5.027344, mu_down = -25.949219, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   4.261217%
mu_up =   5.027344, mu_down = -25.974609, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   3.924710%
mu_up =   5.027344, mu_down = -25.987305, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   3.729372%
mu_up =   5.027344, mu_down = -25.993652, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   3.618478%
mu_up =   5.027344, mu_down = -25.996826, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   3.556423%
mu_up =   5.027344, mu_down = -25.998413, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   3.521504%
mu_up =   5.027344, mu_down = -25.999207, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   3.501982%
mu_up =   5.027344, mu_down = -25.999603, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   3.491764%
mu_up =   5.027344, mu_down = -25.999802, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   3.478129%
mu_up =   5.027344, mu_down = -25.999901, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   3.520766%
mu_up =   5.027344, mu_down = -25.999950, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   3.509227%
mu_up =   5.027344, mu_down = -25.999975, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   3.977577%
mu_up =   5.027344, mu_down = -25.999988, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   1.952098%
mu_up =   5.027344, mu_down = -25.999994, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   0.189328%
mu_up =   5.027344, mu_down = -25.999997, N_up =   8.000000, N_down =   0.000000, Delta_Int_error =   0.032483%
Converged!