r0=4000/(2.6*2.6*58*pi);f=1.2*10^7;

a=2.6*pi*10^(-3)*sqrt(f*5.8);
x=a-a*i;
z=x*bessela(0,x)/bessel(1,x)
z1=z*sqrt(2)
