회전 타원체의 반자계 계수

m = a/b;
if(m==1) {m=1.0001;}
L = (1/(m^2-1))*((m/(2*sqrt(m^2-1)))*ln((m+sqrt(m^2-1))/(m-sqrt(m^2-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);