% Code revision 28 -- change plotting order of n to make it prettier.

% Master call signature:
%    Master(N_up, N_down, U, kT, Nc, omega, MaxIterations);

>> Master(10,5,-3,.01,20,1,200);     
    
    % Plot 12.

>> Master(10,9,-3,.01,20,1,200);     

    % Plot 13. Large condensate, very hard to make Ns different by adjusting mus. Weird.

>> Master(10,9,-3,.1 ,20,1,200);     

    % Still condenses.

>> Master(10,9,-3,1,20,1,200);       

    % Plot 14. Still condenses; Ns not quite equal (but constant ratio?). 
    % Convergence slower -- almost at Tc? 

>> Master(10,9,-3,2,20,1,200);       

    % Plot 15. Delta now very small and asymmetric (i.e. T > Tc).

>> Master(10,9,-3,3,20,1,200);       

    % Delta now about 10^-6, similar shape. T > Tc.

>> Master(10,9,-5,1,20,1,200);       

    % Plot 16. Around iteration 120, mu_up ~ 9 and mu_down ~ 2, 
    % BCS state in centre and halo of excess particles.
    % After a while (iteration ~180), mu_down increases to 4 and all of a sudden n_up(x) == n_down(x) everywhere again.
    % Perturbations (every time i%100 = 0) don't seem to have much effect ( because  delta != 0  ) -- this is good

>> Master(10,9,-5,2,20,1,200);       

    % Plot 17. Compare to plot 15. Doubled interaction strength -> higher T_c. In this plot, T < Tc.

>> Master(10,9,-5,3,20,1,500);       
    
    % Still condensed (eventually). for U = -3, Tc ~ 1. For U = -5, Tc ~ 4.

>> Master(10,9,-5,4,20,1,500);
 
    % Looks as though it would condense after several thousand iterations. This must be almost on Tc.
    % Delta is quite small (5E-3) and rises very slowly (about 1% per iteration).

>> Master(10,9,-5,5,20,1,500);

    % Definitely above Tc now - Delta gradually shrinking.