Most of these exercises are of an advanced nature. However, any of the exercises from Chap...

Most of these exercises are of an advanced nature. However, any of the exercises from Chapter 8 may be done using odesolve, so look there for more elementary exercises.

The system

models a chemical reaction during which two components, x2 and x3, undergo oscillatory growth and decay.

Enter the equations in odesolve, then set initial conditions x1 (0) = 500 and x2(0) = x3(0) = x4(0) = 0.

a) Set the relative tolerance to le-4 and plot x2 and x3 versus time.

b) If you attempt to plot all four components versus time, they scale badly in the odesolve window. Try it! Export the solution data to the command window workspace. Take careful notice of the name of the structure that is exported. If you have not exported anything during this session of MATLAB it will be odedatal. We will assume that is the case, but the output variable could be odedata* where * is a small integer. Enter odedatal at the MATLAB prompt to see the form of the structure and a list of its fields. Assuming that the fields are t, x1, x2, x3, and x4, set x1 = odedatal.x1/200; x2 = odedatal.x2; x3 = odedatal.x3; x4 = odedatal.x4/200; and t = odedatal. t;. Notice that we are rescaling x1 and x4. Plot x1, x2, x3, and x4 versus time in a single figure. Label your axes and provide appropriate labels, title, and legend.