gcall.m

用于电力系统的一个很好的分析软件

function gcall(a)

global DAE

if ~a.n, return, end

V = DAE.y(a.vbus);
B = a.u.*bsvc(a);

DAE.g = DAE.g ...
        - sparse(a.vbus,1,DAE.y(a.q),DAE.m,1) ...
        + sparse(a.q,1,B.*V.*V-DAE.y(a.q),DAE.m,1) ...
        + sparse(a.vref,1,a.con(:,8)-DAE.y(a.vref),DAE.m,1);