IndCalc = (a, b, n) => {
const f1 = x => (1 + x * (0.383901 + 0.017108 * x)) / (1 + 0.258952 * x);
const f2 = x => (x * (0.093842 + x * (0.002029 - x * 0.000801)));
const u0 = 0.4 * Math.PI * 1e-06;
const n2 = n * n;
const shape_ratio = 2 * a / b;
const shape_ratio_pwr2 = shape_ratio * shape_ratio;
if (shape_ratio <= 1) {
return u0 * n2 * Math.PI * a * a
* (f1(shape_ratio_pwr2) - (4 / (3 * Math.PI)) * shape_ratio)
/ b;
}
return u0 * n2 * a
* ((Math.log(4 * shape_ratio) - 0.5) * f1(1 / shape_ratio_pwr2) + f2(1 / shape_ratio_pwr2));
}