电力系统分析计算程序

				function gcall(a)								global DAE								if ~a.n, return, end								DAE.g = DAE.g ...				        - sparse(a.vbus,1,a.u.*DAE.x(a.ist).*DAE.y(a.vbus),DAE.m,1) ...				        + sparse(a.vref,1,a.Vref-DAE.y(a.vref),DAE.m,1);