Function SingularityFreeJacobiPRecursion

src/Numerical/Functions.cpp:49–52

}

double SingularityFreeJacobiPRecursion(
    double x, double y, const std::array<double, 5>& cm, double pm1, double pm2) {
  return (cm[0] * (cm[1] * y + cm[2] * x) * pm1 - cm[3] * y * y * pm2) / cm[4];
}

double SingularityFreeJacobiP(unsigned n, unsigned a, unsigned b, double x, double y) {
  if (n == 0) {