stateestimation-self-derivatives-measurementdistribution/LineCurrent.m

function [ output_args ] = LineCurrent( Linei,Linej,Liner,Linex,Volt,VAngle )
cmpV=Volt.*exp(1j*VAngle);
GB=1./(Liner+1j*Linex);
cmpI=( cmpV(Linei)-cmpV(Linej) ).*GB;
output_args=abs(cmpI);
end