*) k2
program MyDip
fr (alpha)! функція
(4), PARAMETER :: i=(0.0,1.0)! уявна одиниця (4) :: fr, alpha
complex (4) :: eps, eps1, eps2, eps3, ka2, del11, del21, k11, k12, k22, psi, gam1, sig1, sig2! вони комплексні
! real :: mu1=1, mu2=2.6, nu1=0.16, nu2=0.25, h=0.1, h1=1, ro2=3.9, ro1=1, omega! старі параметри
! real :: 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 набір
common/global/omega
eps1=(mu2 * (h ** 3))/(6 * (1-nu2))
eps2=((h ** 3)/12) * ro2 * (omega ** 2)
eps3=h * ro2 * (omega ** 2)
ka2=sqrt ((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=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)
fr=real (k22/del21)
function fr
fi (alpha)! функція
(4), PARAMETER :: i=(0.0,1.0)! уявна одиниця (4) :: fi, alpha
complex (4) :: eps, eps1, eps2, eps3, ka2, del11, del21, k11, k12, k22, psi, gam1, sig1, sig2! вони комплексні
! real :: 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 набір
common/global/omega
eps1=(mu2 * (h ** 3))/(6 * (1-nu2))
eps2=((h ** 3)/12) * ro2 * (omega ** 2)
eps3=h * ro2 * (omega ** 2)
ka2=sqrt ((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=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...
