function fcall(p)global DAEif ~p.n, return, endglobal DAEDAE.f(p.v) = p.u.*(DAE.y(p.If) - p.con(:,6))./p.con(:,2);% anti-windup limitidx = find(DAE.x(p.v) <= 0 & DAE.f(p.v) < 0);DAE.f(p.v(idx)) = 0;idx = find(DAE.x(p.v) >= p.con(:,7) & DAE.f(p.v) > 0);DAE.f(p.v(idx)) = 0;DAE.x(p.v) = min(DAE.x(p.v),p.con(:,7));DAE.x(p.v) = max(p.u.*DAE.x(p.v),0);