add findEf.f90 findEf_Newton.f90 ne_mu.f90. Also change two places in main.f90 wrongly define Delta as a complex variable