d=0.22
r=1700
n=0.84
c=800
U=2.32

p=((pi*(d^2)*r*c)/(86400*n))^0.5


m1=cosh(p+complex(real = 0, imaginary = 1)*p)
m2=d*(((exp(1)^p)-(exp(1)^-p))*cos(p)+((exp(1)^p)+(exp(1)^-p))*sin(p)-complex(real = 0, imaginary = 1)*((exp(1)^p)-(exp(1)^-p))*cos(p)+complex(real = 0, imaginary = 1)*((exp(1)^p)+(exp(1)^-p))*sin(p))/(4*n*p)
m3=(n*p*(((exp(1)^p)-(exp(1)^-p))*cos(p)+((exp(1)^p)+(exp(1)^-p))*sin(p)+complex(real = 0, imaginary = 1)*((exp(1)^p)-(exp(1)^-p))*cos(p)+complex(real = 0, imaginary = 1)*((exp(1)^p)+(exp(1)^-p))*sin(p)))/(2*d)+3.94