cos[2*(sqrt((6.67e-11)*(1.3e32)/(100000^3)))(t-4e16)]