The systemmodels a chemical reaction called an oregonator. Suppose that ε = 10-2 and q = 9...

The system

models a chemical reaction called an oregonator. Suppose that ε = 10-2 and q = 9 × 10-5. Put the system into normal form and write an ODE function M-file for the system that passes f as a parameter. The idea is to vary the parameter f and note its affect on the solution of the oregonator model. We will use the initial conditions x (0) = 0.2 and y (0) = 0.2 and the solution interval [0, 50].

a) Use ode45 to solve the system with f = 1/4. This should provide no difficulties.

b) Use ode45 to solve the system with f = 1. This should set off all kinds of warning messages.

c) Try to improve the accuracy with options = odeset('RelTol', 1e-6) and using options as the options parameter in ode45. You should find that this slows computation to a crawl as the system is very stiff.

d) The secret is to use ode 16s instead of ode45. Then you can note the oscillations in the reaction.