program

from numpy import* from matplotlib.pyplot import* def f(p,e): p1=(-1923.14/(5c5))(23.14*G/p)(5/3)(m1m2/(m1+m2)(1/3))*(1+(73/24)*e2+(37/96)e**4)(1-e2)(-7/2) return p1

def g(p,e): e1=(-6083.14/(15(c5))) e2=e/p e3=(23.14G/(p))(5/3) e4=m1m2/((m1+m2)(1/3)) e5=(1+((121/304)*e2))/((1-e2)(5/2)) e7=e1e2e3e4*e5

#e1=(-608*3.14*e/(15*c**5*p))*(m1*m2/(m1+m2)**(1/3))*(1+(121/304)*e**2)*(1-e**2)**(-5/2)
return e7

plist=[] elist=[] tlist=[] m1=101.9891030 m2=101.9891030 c=310**8 G=6.674310**-11 e=0.1 p=0.1 t=0.00 h=.10 tlist.append(t) elist.append(e) plist.append(p)

#while(e<1.0): while (t<=100): e=e+hf(p,e) p=p+hf(p,e) t=t+h tlist.append(t) elist.append(e) plist.append(p)

figure(1) plot (tlist,plist)

grid("true") xlabel("time period") xlim(0,7) ylim(0.00,.10) ylabel("orbital period") #show() #figure(2) #plot(tlist,elist) xlim(0,7) ylim(0.00,.10)
show()