choose time of lastscattering to be 300,000 years. Now integration is stopped at t=t_last_scattering. Also added function get_angles. Prac_healpix.py is used for spherical harmominc calculations. It need healpix and healpy to be installed