gycall.m

来自「电力系统分析计算程序」· M 代码 · 共 8 行

function Gycall(a)

global DAE

if ~a.n, return, end

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