vPerihelion function

double vPerihelion(
  1. double a,
  2. double e
)

Velocity at perihelion for eccentricity e and semimajor axis a.

Implementation

double vPerihelion(double a, double e) {
  return 0.0172021 * math.sqrt((1 + e) / (a * (1 - e)));
}