电力系统分析计算程序

				function gcall(a)								global DAE								if ~a.n, return, end								V = DAE.y(a.vbus);				B = a.u.*bsvc(a);								DAE.g = DAE.g ...				        - sparse(a.vbus,1,DAE.y(a.q),DAE.m,1) ...				        + sparse(a.q,1,B.*V.*V-DAE.y(a.q),DAE.m,1) ...				        + sparse(a.vref,1,a.con(:,8)-DAE.y(a.vref),DAE.m,1);