double distortionFactor(double rSquared) { double rFactor = 1.0; double factor = 1.0; for (final k in coefficients) { rFactor *= rSquared; factor += k * rFactor; } return factor; }