gcall.m

这是一个很适合研究和学习用的电力系统仿真软件

function gcall(a)

global Bus DAE

if ~a.n, return, end

V = DAE.V(a.bus);
B = bsvc(a);

DAE.gq = DAE.gq + sparse(a.bus,1,-B.*V.*V,Bus.n,1);