**** Computation of Chyuthi Samskaram = Manda Kendra = Madhyamam – Mandochha ****
**** Approximation to Theta, the True Anomaly
y = year(ctod(mdate))-1899
t = y/100
e = .01675134-0.0000004180*y
R = 206265
mdk = (apk+180)%360
a = R(2e)
b = R(1/4)e^3
c = R1
d = R2
l = R5
m = R6
mjyavk =(a-b+c)sina +(d-f+g)sin2a + (h-i)sin3a + (j-k)sin4a + (l sin5a) +(msin6a)
suy = abs(mjyavk/60)
sdiff = suy/60
Computation of Ravi Mandakarna, heliocentric distance *
**** Approximation to r, the radius vector *
smaxis = 10000
rv1 = (1+(1/2)e^2) rv2 = e (1 – (3/8)*e^2) * cos(dtor(mdk))
rv3 = .5*e^2 * ( 1- (2/3)* e^2 )* cos(dtor(2mdk)) rv4 = (3/8) e^3* cos(dtor(3mdk)) rv5 = (1/2)e^4cos(dtor(4mdk))
rmk = smaxis* (rv1 – rv2 – rv3 – rv4 – rv5)
if apk > 180
rsp = np + sdiff
else
rsp = np -sdiff
endif
rsp = rsp % 360
ncrsp = (rsp – ayd)%360
- 5/96)e**5)
sina = sin(dtor(mdk[↩] - 5/4)e^2)
f = R((11/24)e^4)
g = R((17/192)e^5)
sin2a = sin(dtor(mdk*2[↩] - 13/12)e^3)
i = R((43/64)e^5)
sin3a = sin(dtor(mdk3[↩] - 103/98)e^4) k = R((451/480)e^6) sin4a = sin(dtor(mdk4[↩]
- 1097/960)e^5)
sin5a = sin(dtor(mdk*5[↩] - 1223/960)e^6)
sin6a = sin(dtor(mdk*6[↩]
https://ajiio.in/woslXgV