function Gycall(a)global DAEif ~a.n, return, endm = DAE.x(a.m);V = DAE.y(a.vbus);a1 = a.con(:,11);a2 = a.con(:,12);DAE.Gy = DAE.Gy + sparse(a.bus,a.vbus, ...          a.u.*a.con(:,9).*a1.*(V.^(a1-1))./(m.^a1),DAE.m,DAE.m);DAE.Gy = DAE.Gy + sparse(a.vbus,a.vbus, ...          a.u.*a.con(:,10).*a2.*(V.^(a2-1))./(m.^a2),DAE.m,DAE.m);