exp(sig2*h1)-cdexp(-(sig2*h1)))/2)*((cdexp(sig1*h1)+cdexp(-(sig1*h1)))/2)))/(- 4)
del21=(sig1*sig2*((1/4)*(ka2**4)-((gam1+(alpha**2))**2)+((1/4)*(ka2**4)+((gam1+(alpha**2))**2))*((cdexp(sig1*h1)+cdexp(-(sig1*h1)))/2)*((cdexp(sig2*h1)+cdexp(-(sig2*h1)))/2))-(alpha**2)*((gam1**2)+(sig1**2)*(sig2**2))*((cdexp(sig1*h1)-cdexp(-(sig1*h1)))/2)*((cdexp(sig2*h1)-cdexp(-(sig2*h1)))/2))/2
del11=((cdexp (sig2 * h1) + cdexp (- (sig2 * h1)))/2)
! f=real (((k22/del21)/((k22/del21) * psi + 1)) * cdexp (-I * alpha * x1))! для переміщень
f=real (((1)/((k22/del21) * psi + 1)) * cdexp (-I * alpha * x1))! для напруг
end function f
ДОДАТОК B
Лістинг програми в середовищі Maple15 (побудова графіків нулів і полюсів).
(stats);
:=readdata (`c: MOmegaNuli.txt`, 1);:=readdata (` c: MNuli.txt`, 1);:=readdata (`c: MOmegaPolusa. txt`, 1);:=readdata (`c: MPolusa.txt`, 1); [scatterplot] (x, y); [scatterplot] (z, q);
ДОДАТОК Г
Лістинг програми в середовищі Maple15 (побудова графіків переміщень і напруг).
(stats);
[listplot] (readdata (`c: xf19n.txt`, float, 2), color=red, labels=[ x raquo ;, u ], labelfont=[TIMES , BOLDITALIC, 11]); [listplot] (readdata (`c: yf14.txt`, float, 2), color=red, labels=[ x raquo ;, u ], labelfont=[TIMES, BOLDITALIC, 11]);