Changed findphi60 module to solve an integral equation. Also added a findphiUsingphi0 module that begins the integration at n=0 instead of n=60. It seems to work a little better, but gives negative results for the spectral index.