11 py4.py 文件源码

def plot(self):
        pl.figure(figsize = (8, 8))
        pl.plot(self.t,self.omega, 'c')
        pl.ylim(-20, 60)
        pl.xlim(0, 10)
        pl.ylabel('$\\omega$ (radians/yr)')
        pl.xlabel('time (yr)')
        pl.title('Hyperion $\\omega$ versus time')
        pl.text(4.1, 55, 'Elliptical orbit')