MATLAB Mathematics - SERC - Index of

Initial Value Problems for ODEs and DAEs
function yp = odefcn(x,y)
yp = ((1 - x)*y - y^2)/epsilon;
end
end % kneeode
The derivative function is defined within nested function odefcn. The value of
epsilon used in odefcn is obtained from the outer function:
function yp = odefcn(x,y)
yp = ((1 - x)*y - y^2)/epsilon;
end
The demo solves the problem using the ode15s function, first with the default
options, and then by imposing a nonnegativity constraint. To run the demo,
type kneeode at the MATLAB command prompt.
Here is the output plot. The plot confirms correct solution behavior after
imposing constraints.
5-41