gycall.m

电力系统分析计算程序

function Gycall(a)

global DAE

if ~a.n, return, end

DAE.Gy = DAE.Gy - sparse(a.ws,a.ws,1,DAE.m,DAE.m);