电力系统分析计算程序

				function Gycall(a)								global DAE								if ~a.n, return, end								m = 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);