function out=radiotst5(zsrc,p);
Co=1844;
Cm=1917.4;
Zm=2023.6;
Zs=560.23; %HiHi
J=.0177068;
lam=3.833e-12;
T1=exp(-(Zm-zsrc)./(2*Zs));
T2=exp(zsrc/(2*Zs));
A=sqrt(1+4*Zs/p(1));
T3=sinh(A*log(T2));
T4=sinh(A*(-1)*log(T1));
T5=T1.*sinh(A*log(T2));
T6=T2.*sinh(sqrt(1+4*Zs/p(1))*(-1)*log(T1));
T7=sinh(sqrt(1+4*Zs/p(1))*Zm/(2*Zs));
T8=J*p(1);
out =J*p(1)+((Cm-J*p(1))*T1.*sinh(A*log(T2))+(Co-J*p(1))*T2.*sinh(sqrt(1+4*Zs/p(1))*(-1)*log(T1)))/sinh(sqrt(1+4*Zs/p(1))*Zm/(2*Zs));