clear; clc; D=6; r=.015; e=8.85*10^(-12); a=complex(0.5,-0.866); ca=log(r/(2*D))-(log(D/(2*r))*a); cb=(log((D^2)/(r^2))*log(r/(2*D)))-(log((2*D)/r)*log(D/(2*r))); Cyr= (2*%pi*e*(10^9)*ca)/cb; mprintf(" capacitance in microF/Km is %.3f %.3fi\n",real(Cyr),imag(Cyr)); I=%i*2*%pi*50*Cyr*110000*10^-6; mprintf(" Current is %.3f+%.3fi A",real(I),imag(I));