align="justify"> gam1=(alpha ** 2) - 0.5 * (ka2 ** 2)
sig2=sqrt ((alpha ** 2) - (ka2 ** 2))
sig1=sqrt ((alpha ** 2) -eps * (ka2 ** 2))
k11=(sig2*(ka2**2)*(sig1*sig2*((exp(sig1*h1)+exp(-(sig1*h1)))/2)*((exp(sig2*h1)-exp(-(sig2*h1)))/2)-(alpha**2)*((exp(sig1*h1)-exp(-(sig1*h1)))/2)*((exp(sig2*h1)+exp(-(sig2*h1)))/2)))/(- 4)
k12=(sig1*sig2*(gam1+(alpha**2))*(((exp(sig1*h1)+exp(-(sig1*h1)))/2)*((exp(sig2*h1)+exp(-(sig2*h1)))/2)- 1)-(gam1*(alpha**2)+(sig1**2)*(sig2**2))*((exp(sig1*h1)-exp(-(sig1*h1)))/2)*((exp(sig2*h1)-exp(-(sig2*h1)))/2))/2
k22=(sig1*(ka2**2)*(sig2*sig1*((exp(sig2*h1)+exp(-(sig2*h1)))/2)*((exp(sig1*h1)-exp(-(sig1*h1)))/2)-(alpha**2)*((exp(sig2*h1)-exp(-(sig2*h1)))/2)*((exp(sig1*h1)+exp(-(sig1*h1)))/2)))/(- 4)
del21=(sig1*sig2*((1/4)*(ka2**4)-((gam1+(alpha**2))**2)+((1/4)*(ka2**4)+((gam1+(alpha**2))**2))*((exp(sig1*h1)+exp(-(sig1*h1)))/2)*((exp(sig2*h1)+exp(-(sig2*h1)))/2))-(alpha**2)*((gam1**2)+(sig1**2)*(sig2**2))*((exp(sig1*h1)-exp(-(sig1*h1)))/2)*((exp(sig2*h1)-exp(-(sig2*h1)))/2))/2
del11=((exp (sig2 * h1) + exp (- (sig2 * h1)))/2)
znami=aimag ((k22/del21) * psi + 1)
end function znami
ДОДАТОК Б
Лістинг програми на мові Fortran. Побудова переміщень
Program zer1
use dfimsl/vab/omega, x1 (8), external :: f !,U3,U4(8)::omega,x1,alpha,alpha11,errabs,errel,res,res1,errest,Int1,Int2,a(8)::Int3(4)::i,j,maxfn,maxfn1,ff
!open(1,action=laquo;writeraquo;,file=laquo;xf.txtraquo;,form=laquo;formattedraquo;,position=laquo;rewindraquo;,status=laquo;oldraquo;) ! файл для перемещений(1,action=laquo;writeraquo;,file=laquo;yf.txtraquo;,form=laquo;formattedraquo;,position=laquo;rewindraquo;,status=laquo;oldraquo;) ! файл дл звернень
=0.0d0=0.01d0
do x1=0d0,10d0,0.1
omega=9.3
! omega=2.100009
do a=- 30D0,10D0,1.0
alpha=a
alpha11=alpha + 0.1
call dqdags (f, alpha, alpha11, errabs, errel, res, errest)
Int1=Int1 + res
end do
write (1, (f10.5,2x, 15 (f8.5,1x)) ) x1, Int1
Int1=0.0
end do
close (1) program zer1
f (alpha)! функція
/vab/omega, x1 (8), PARAMETER :: I=(0.0,1.0)! уявна одиниця (8) :: f, alpha, x1
complex (8) :: eps, eps1, eps2, eps3, ka2, del11, del21, k11, k12, k22, psi, gam1, sig1, sig2, elem! вони комплексні
real (8) :: mu1=1, mu2=2.6, nu1=0.3, nu2=0.3, h=0.1, h1=1, ro2=3.8, ro1=1, omega! 1 набір
! real :: mu1=1, mu2=2.6, nu1=0.3, nu2=0.3, h=0.01, h1=1, ro2=3.8, ro1=1, omega! 2 набір
! real :: mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.1, h1=1, ro2=3.8, ro1=1, omega! 3 набір
! real :: mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.01, h1=1, ro2=3.8, ro1=1, omega! 4 набір
! real :: mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.1, h1=1, ro2=1, ro1=1, omega! 5 набір
! real :: mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.01, h1=1, ro2=1, ro1=1, omega! 6 набір
eps1=(mu2 * (h ** 3))/(6 * (1-nu2))
eps2=((h ** 3)/12) * ro2 * (omega ** 2)
eps3=h * ro2 * (omega ** 2)
ka2=cdsqrt (dcmplx ((ro1 * (omega ** 2))/mu1))
eps=(1-2 * nu1)/(2-2 * nu1)
psi=eps1 * (alpha ** 4) -eps2 * (alpha ** 2) -eps3
gam1=(alpha ** 2) - 0.5 * (ka2 ** 2)
sig2=cdsqrt (dcmplx ((alpha ** 2) - (ka2 ** 2)))
sig1=cdsqrt (dcmplx ((alpha ** 2) -eps * (ka2 ** 2)))
k11=(sig2*(ka2**2)*(sig1*sig2*((cdexp(sig1*h1)+cdexp(-(sig1*h1)))/2)*((cdexp(sig2*h1)-cdexp(-(sig2*h1)))/2)-(alpha**2)*((cdexp(sig1*h1)-cdexp(-(sig1*h1)))/2)*((cdexp(sig2*h1)+cdexp(-(sig2*h1)))/2)))/(- 4)
k12=(sig1*sig2*(gam1+(alpha**2))*(((cdexp(sig1*h1)+cdexp(-(sig1*h1)))/2)*((cdexp(sig2*h1)+cdexp(-(sig2*h1)))/2)- 1)-(gam1*(alpha**2)+(sig1**2)*(sig2**2))*((cdexp(sig1*h1)-cdexp(-(sig1*h1)))/2)*((cdexp(sig2*h1)-cdexp(-(sig2*h1)))/2))/2
k22=(sig1*(ka2**2)*(sig2*sig1*((cdexp(sig2*h1)+cdexp(-(sig2*h1)))/2)*((cdexp(sig1*h1)-cdexp(-(sig1*h1)))/2)-(alpha**2)*((cd...
