function Fxcall(a)global DAEif ~a.n, return, endIk = DAE.x(a.Ik);Vk = DAE.x(a.Vk) + (~a.u);Vs = DAE.y(a.vbus);pH2 = DAE.x(a.pH2) + (~a.u);pH2O = DAE.x(a.pH2O) + (~a.u);pO2 = DAE.x(a.pO2) + (~a.u);qH2 = DAE.x(a.qH2);m = DAE.x(a.m);Sn = a.con(:,2);Vn = a.con(:,3);Te = a.con(:,4);TH2 = a.con(:,5);KH2 = a.con(:,6);Kr = a.u.*a.con(:,7);iKr = 1./a.con(:,7);TH2O = a.con(:,8);KH2O = a.con(:,9);TO2 = a.con(:,10);KO2 = a.con(:,11);rHO = a.con(:,12);Tf = a.con(:,13);Uopt = a.con(:,14);r = a.u.*a.con(:,17);N0 = a.u.*a.con(:,18);E0 = a.con(:,19);RTon2F = a.con(:,20);xt = a.con(:,26);Km = a.con(:,27);Tm = a.con(:,28);Pref = a.con(:,21).*a.con(:,23);u = m <= a.con(:,29) & m >= a.con(:,30) & a.u;Vt = m.*Vk.*a.con(:,24);[Input,Umax,Umin] = fcinput(a);sq = sqrt(1-(xt.*Ik./Vs.*Vn./Sn./m).^2);dQdi = 0.5*Vs.*Vt./xt.*(-2*((xt.*Vn./Sn./m./Vs).^2).*Ik)./sq;dQdv = m.*a.con(:,24).*Vs./xt.*sq;dQdm = Vk.*a.con(:,24).*Vs./xt.*sq + ...       0.5*Vs.*Vt./xt.*(2*((xt.*Ik.*Vn./Sn./Vs).^2)./(m.^3))./sq;u1 = (Input > Umax | Input < Umin) & a.u;u2 = ~u1 & ~a.con(:,25) & a.u;DAE.Fx = DAE.Fx - sparse(a.Ik,a.Ik,1./Te,DAE.n,DAE.n);DAE.Fx = DAE.Fx - sparse(a.Vk,a.Vk,1e3,DAE.n,DAE.n);DAE.Fx = DAE.Fx - sparse(a.m,a.m,1./Tm,DAE.n,DAE.n);DAE.Fx = DAE.Fx - sparse(a.pH2,a.pH2,1./TH2,DAE.n,DAE.n);DAE.Fx = DAE.Fx - sparse(a.pH2O,a.pH2O,1./TH2O,DAE.n,DAE.n);DAE.Fx = DAE.Fx - sparse(a.qH2,a.qH2,1./Tf,DAE.n,DAE.n);DAE.Fx = DAE.Fx - sparse(a.pO2,a.pO2,1./TO2,DAE.n,DAE.n);DAE.Fx = DAE.Fx + sparse(a.Ik,a.qH2,0.5*u1.*a.con(:,16).*iKr,DAE.n,DAE.n);DAE.Fx = DAE.Fx - sparse(a.Ik,a.Vk,u2.*Sn.*Pref./(Vk.^2)./Vn./Te,DAE.n,DAE.n);DAE.Fx = DAE.Fx - sparse(a.Vk,a.Ik,1e3*r./Vn,DAE.n,DAE.n);DAE.Fx = DAE.Fx + sparse(a.Vk,a.pH2,1e3*N0.*RTon2F./pH2./Vn,DAE.n,DAE.n);DAE.Fx = DAE.Fx - sparse(a.Vk,a.pH2O,1e3*N0.*RTon2F./pH2O./Vn,DAE.n,DAE.n);DAE.Fx = DAE.Fx + sparse(a.Vk,a.pO2,1e3*0.5*N0.*RTon2F./pO2./Vn,DAE.n,DAE.n);DAE.Fx = DAE.Fx - sparse(a.pH2,a.Ik,2.*Kr./KH2./TH2,DAE.n,DAE.n);DAE.Fx = DAE.Fx + sparse(a.pH2,a.qH2,a.u./KH2./TH2,DAE.n,DAE.n);DAE.Fx = DAE.Fx + sparse(a.pH2O,a.Ik,2.*Kr./KH2O./TH2O,DAE.n,DAE.n);DAE.Fx = DAE.Fx - sparse(a.pO2,a.Ik,Kr./KO2./TO2,DAE.n,DAE.n);DAE.Fx = DAE.Fx + sparse(a.pO2,a.qH2,a.u./rHO./KO2./TO2,DAE.n,DAE.n);DAE.Fx = DAE.Fx + sparse(a.qH2,a.Ik,2.*Kr./Uopt./Tf,DAE.n,DAE.n);DAE.Fy = DAE.Fy - sparse(a.m,a.vbus,u.*Km./Tm,DAE.n,DAE.m);DAE.Gx = DAE.Gx - sparse(a.vbus,a.Vk,a.u.*dQdv,DAE.m,DAE.n);DAE.Gx = DAE.Gx - sparse(a.vbus,a.Ik,a.u.*dQdi,DAE.m,DAE.n);DAE.Gx = DAE.Gx - sparse(a.vbus,a.m,u.*dQdm,DAE.m,DAE.n);DAE.Gx = DAE.Gx - sparse(a.bus,a.Ik,a.u.*Vk.*Vn./Sn,DAE.m,DAE.n);DAE.Gx = DAE.Gx - sparse(a.bus,a.Vk,a.u.*Ik.*Vn./Sn.*a.con(:,24),DAE.m,DAE.n);