Personel Information

BSc: Physics Department, Middle East Technical University (METU), Ankara, Turkey
MSc: Physics Department, METU
PhD: Qualified, Physics Department, METU


Monday, July 11, 2016

FACTORIAL FOR KEPLER’S MOTION


##Constants and initializations
ecc=0.1; ## Eccentricity
M=24.851090; ## Mean anomaly(degrees)
Mr=M*pi/180; ## Mean anomaly (rads)
ser=[]; ## Series in E
jser=[]; ## Bessel function
for m=1:100
for n=0:100
jser(n+1)=((-1)^n)*((m*ecc/2)^(2*n+m))/(factorial(n)*factorial(m+n));
jser(n+1)+=jser(n+1);
endfor
ser(m)=(1/m)*jser(n+1)*sin(m*Mr);
ser(m)+=ser(m);
endfor
E=M+2*ser(m);
save -text FACTORIALKEPLER.dat
printf('We find E to be %f. This is the answer of Q1(c)\n',E);

No comments:

Post a Comment