calculation of K with astropy
rho with n=1.5 fixed