function Gycall(p)global DAE Settingsif ~p.n, return, endDAE.J11(p.bus,:) = 0;DAE.J12(p.bus,:) = 0;DAE.J11(:,p.bus) = 0;DAE.J12(:,p.bus) = 0;DAE.J11(p.bus,p.bus) = DAE.J11(p.bus,p.bus) + speye(p.n);if Settings.pv2pq  idx = p.bus(find(~p.dq));else  idx = p.bus;endif isempty(idx),return, endn = length(idx);DAE.J21(idx,:) = 0;DAE.J22(idx,:) = 0;DAE.J12(:,idx) = 0;DAE.J22(:,idx) = 0;DAE.J22(idx,idx) = DAE.J22(idx,idx) + speye(n);