calc4.m

可以计算电力系统概率潮流

function pci=CalC4h(sgm, u1, u2, u3, u4)

% calculate Gram Charles index 
pci=zeros(1,4);

if sgm~=0
    pci(1,1)=0;
    pci(1,2)=0;
    pci(1,3)=-u3/(sgm^3);   
    pci(1,4)=u4/(sgm^4)-3;
    
end

return