summaryrefslogtreecommitdiff
path: root/2409/CH2/EX2.14/Ex2_14.sce
diff options
context:
space:
mode:
Diffstat (limited to '2409/CH2/EX2.14/Ex2_14.sce')
-rwxr-xr-x2409/CH2/EX2.14/Ex2_14.sce17
1 files changed, 17 insertions, 0 deletions
diff --git a/2409/CH2/EX2.14/Ex2_14.sce b/2409/CH2/EX2.14/Ex2_14.sce
new file mode 100755
index 000000000..cff340896
--- /dev/null
+++ b/2409/CH2/EX2.14/Ex2_14.sce
@@ -0,0 +1,17 @@
+
+//Variable Declaration
+pi = %pi
+n=14.2171401*2*%pi/86400 //Mean motion (rad/sec)
+M=204.9779+0.001*180*5/pi //Mean anomaly(rad)
+e=9.5981*10**-3 //Eccentricity
+a=7194.9 //Semimajor axis(km)
+
+//Calculation
+
+v=(M*pi/180)+2*e*sin(M*pi/180)+(5*e**2*sin(2*M*pi)/(4*180)) //True Anomaly (radians)
+v=v*180/%pi //True anomaly(degrees)
+r=a*(1-e**2)/(1+e*cos(v)) //Magnitude of radius vector after 5s(km)
+
+//Results
+printf("The true anomaly is: %.3f degrees",v)
+printf("\nThe magnitude of radius vector 5s after epoch is: %d km",r)