Version 1.5 ; modified the variable names; Using gap equation to determin the gap; Modified particle number determine skeme, change the 'Fermi energy' of Bogoliubov qusiparticles to 0