旋转椭圆体的反磁场系数

m = a/b;
if(m==1) {m=1.0001;}
L = (1/(m²-1))*((m/(2*√(m²-1)))*ln((m+√(m²-1))/(m-√(m²-1)))-1);
M = (m/(2*(m^2-1)))*(m - (1/(2*sqrt(m^2-1)))*ln((m+sqrt(m^2-1))/(m-sqrt(m^2-1))));
μeff = 1/(L+(1-L)/μr);
μeffmax = 1/L;
print(a);
print(b);
print(m);
print(L);
print(μr);
print(μeff);
print(μeffmax);
print(M);