def dVS(th1, th2):
enpt = entrance(Rsoi['SATURN'], st0S[3:], th1, th2)
exst = swingby(Rsoi['SATURN'], np.concatenate((enpt, st0S[3:])), GMS)
return np.linalg.norm(stfS[3:] - exst[3:])
def dVN(th1, th2):
enpt = entrance(Rsoi['NEPTUNE'], st0N[3:], th1, th2)
exst = swingby(Rsoi['NEPTUNE'], np.concatenate((enpt, st0N[3:])), GMN)
return np.linalg.norm(stfN[3:] - exst[3:])
def dVJ(th1, th2):
enpt = entrance(Rsoi['JUPITER'], st0J[3:], th1, th2)
exst = swingby(Rsoi['JUPITER'], np.concatenate((enpt, st0J[3:])), GMJ)
return np.linalg.norm(stfJ[3:] - exst[3:])
def dVU(th1, th2):
enpt = entrance(Rsoi['URANUS'], st0U[3:], th1, th2)
exst = swingby(Rsoi['URANUS'], np.concatenate((enpt, st0U[3:])), GMU)
return np.linalg.norm(stfU[3:] - exst[3:])