static double saLonT(double w) { double t, v = 628.3319653318; t = (w - 1.75347 - math.pi) / v; v = ev(t); t += (w - saLon(t, 10)) / v; v = ev(t); t += (w - saLon(t, -1)) / v; return t; }