11 py8.py 文件源码

def plot_delta(self):
        pl.figure(figsize = (8, 8))
        pl.semilogy(self.tprime, self.deltatheta, 'r')
#        pl.ylim(0.0001, 0.1)
#        pl.ylim(0.0001, 10)
        pl.xlim(0, 100)
        pl.ylabel('$\\Delta\\theta$ (radians)')
        pl.xlabel('time (yr)')
        pl.title('Hyperion $\\theta$ versus time')
        pl.text(4.1, 3e3, 'Ellipitical orbit')