Hi,
The Matlab Command to obtain Corona Losses in a given transmission line is as follows:

clc

clear all

close all

% To calculate the corona losses in a transmission line

f= input('Enter the operating frequency of the transmission line in hertz:');

r= input('Enter the radius of conductor in cm:');

d= input('Enter the distance between the conductors in cm:');

b= input('Enter the barometric pressure/cm of Hg:');

t= input('Enter the temperature of the surrounding medium in Celsius:');

Vph= input('Enter the operating phase voltage in KV of the transmission line:');

m= input('Enter the surface factor of the conductor:');

g= 21.21; % KV/Phase

p= 3.92*b/(273+t); % Air density factor

Vc= m*g*r*p*log(d/r) % Critical disruptive voltage

Pc= 241*((f+25)/p)*(r/d)^(1/2)*(Vph-Vc)^2*10^(-5);

if Vph > Vc

    display Pc

    fprintf('\n The power loss due to corona in KW/phase/Km is %f KW/Phase/Km \n', Pc)

else

    display Pc=0

end