(*4.6節 円輪の引張り Sec 4-6 , 表4.1の計算*)

(*式(4.60), (4.61),(4.62)*)

(*式(4.60) Pi*r=p/k[K(p)-F(p,Pi/4)], はじめに p を与えて 荷重項kを求める．*)

(* 荷重P はPi*r/2= から求める*)

Clear["Global`*"];
r=10;EI=20;
Print["ArcSin[p]     P         a          b         MC        MB"];
Do[
thp=ith*10;
If[ith==8, thp=85];
If[ith==9, thp=87];
If[ith==10, thp=89];
If[ith==11, thp=89.9];
If[thp==0, thp=0.01];
p=N[Sin[thp*Pi/180]];
p2=p*p;
k=2/(Pi*r)*p*(EllipticK[p2]-EllipticF[Pi/4,p2]);
LoadP=EI*k^2;
a[ith]=2/(k*p)*(Sqrt[1-p2/2]-Sqrt[1-p2]);
b[ith]=2/(k*p)*((1-p2/2)*(EllipticK[p2]-EllipticF[Pi/4,p2])-(EllipticE[p2]-EllipticE[Pi/4,p2]));
MB=-2/p*(LoadP*EI*(1-p2/2))^(1/2);
MC=-2/p*(LoadP*EI*(1-p2))^(1/2);
Print["    ",thp, " , ",PaddedForm[LoadP,5]," , ",PaddedForm[a[ith],5]," , ",PaddedForm[b[ith],5]," , ",PaddedForm[MC,5]," , ",PaddedForm[MB,5]];
gr[ith]=Graphics[Circle[{0,0},{a[ith],b[ith]}]];
,{ith,0,11,1}];Show[gr[0],gr[1],gr[2],gr[3],gr[4],gr[5],gr[6],gr[7],gr[8],gr[9],gr[10],gr[11]]