Equation for Sheeghra Anomaly —

* SHEEGRA KRIYA ***

sk = ecd-rsp1

if sign(sk) = -1
sk = sk +360
endif

if sk < 180
skh = sk/2
else
skh = (360-sk)/2
endif

str13 = tan(dtor(skh))*10000

shkd = skh – int(skh)
shke = rtod(shkd)
shker = round(shke,0)

kd = kvmkz – rmk1
kt = rmk1+kvmkz
ks = (str13*kd)/kt

KS1 = KS/10000

STR15 = RTOD(ATAN(KS1))