Программа 37a. Кручение прямоуг.сечения (другой способ вычисления коэффициентов)

>

restart;

>	beta:=proc(z) local s,n:

>	s:=sum(1/(2*n-1)^5*tanh((2*n-1)*z*Pi/2),n=1..100):

>	evalf(1-192/z/Pi^5*s)/3:

>

end:

>	N:=100:#Число членов ряда

>	z0:=[1,1.5,1.75,2,2.5,3,4,6,8,10,999]:

Вычисления

>	for i to 11 do

>

z:=z0[i]:

>

k:=2*n-1:

>	bet[i]:=beta(z):

>	s2:=sum(1/k^2/cosh(k*z*Pi/2),n=1..N):

>	tmax:=1-8/Pi^2*s2:

>	alpha[i]:=bet[i]/evalf(tmax):

>	tmin:=-8/Pi^2*sum((-1)^n/k^2*tanh(k*z*Pi/2),n=1..N):

>	gamma1[i]:=evalf(tmin/tmax):

>

od:

Печать результатов

>	x0:=z0,alpha,bet,gamma1:

>	txt:="b/a","alpha","beta","gamma":

>	for j to 4 do

>	printf("%5s |",txt[j]);

>	for i to 11 do

>	printf("%05.3f ",x0[j][i]);

>

od;

>	printf("\n");

>

od;

  b/a |1.000 1.500 1.750 2.000 2.500 3.000 4.000 6.000 8.000 10.000 999.000 

alpha |0.208 0.231 0.239 0.246 0.258 0.267 0.282 0.298 0.307 0.312 0.333 

 beta |0.141 0.196 0.214 0.229 0.249 0.263 0.281 0.298 0.307 0.312 0.333 

gamma |1.000 0.859 0.821 0.795 0.766 0.753 0.745 0.743 0.742 0.742 0.742 

>	save beta, "beta.m";

>

beta(1);