from numpy import * def period(M, R): '''Return the orbital period in seconds of a body in circular orbit of radius R (in metres) around a central body of mass M (in kilograms).''' G = 6.67e-11 return 2 * pi * (R ** 1.5) / sqrt(G * M)