Equation for Saturn in the Great Jupiter Perturbation. ky = year(ctod(mdate))+3102 shaa = ky-4660 shb = mod(shaa,918) shc = shb 12/918 she = shc30

sbj = sin(dtor(she))*10000

h = (sbj*2873)/10000
i = h/60