function glambda(p,lambda,kg)global DAEif ~p.n, return, endjdx = find(p.u);idx = p.bus(jdx);if isempty(idx),return, endDAE.g(idx) = DAE.g(idx) - (lambda+kg*p.con(jdx,11)).*p.pg(jdx);