def mean_lon_of_perigee(jd):
 T = ((jd - jd1950) / 36525.0)
 p = (0.012, 1.65, 6190.67, 1015489.951)
 return (np.polyval(p, T) / 3600.0)

def mean_lon_of_perigee ( jd ) : T = ( ( jd - jd1950 ) / 36525.0 ) p = ( 0.012 , 1.65 , 6190.67 , 1015489.951 ) return ( np . polyval ( p , T ) / 3600.0 )
computes the mean longitude of perigee of the earths orbit at the requested julian date .