* System of nonlinear equations * Steady-state of nonisothermal reactor * Different heats of reaction (DH) can be used. * Parameters SCALAR Aa / 448000 /; SCALAR R / 1.98725 /; SCALAR T0 / 298 /; SCALAR E / -15000 /; SCALAR theta / 3000 /; SCALAR rho / 1 /; SCALAR C0 / 0.003 /; SCALAR Cp / 1 /; SCALAR DH / -50000 /; SCALAR a ; SCALAR b ; SCALAR c ; a = (-rho) * Cp / (C0 * DH); b = theta * Aa; c = E / R; VARIABLES T s objval objective function variable; FREE VARIABLES objval; EQUATIONS f Objective function g1p g1m ; f .. objval =e=s; g1p .. T/T0 - s + b/T0*T*EXP(c/T) - b*(1+a*T0)/(a*T0)*EXP(c/T) =l= 1; g1m .. -T/T0 - s - b/T0*T*EXP(c/T) + b*(1+a*T0)/(a*T0)*EXP(c/T) =l= -1; * Bounds T.LO = 100; T.UP = 1000; * Starting point * T.L = 445.49; MODEL test /ALL/; SOLVE test USING NLP MINIMIZING objval;