11 py10.py 文件源码

def plot(self):
        pl.figure(figsize = (8, 8))
        pl.plot(self.theta,self.omega, 'k.')
#        pl.ylim(-4, 4)
#        pl.xlim(0, 8)
        pl.ylabel('$\\omega$ (radians/yr)')
        pl.xlabel('$\\theta$ (radians)')
        pl.title('Hyperion $\\omega$ versus $\\theta$')
        pl.text(-0.75, 56, 'Elliptical orbit')