The built-in plotting routines of MATLAB’s solvers were illustrated in Exercise 3. They wo...

The built-in plotting routines of MATLAB’s solvers were illustrated in Exercise 3. They work equally well with systems. For example, create

and save as heart.m.

a) Enter ode45(@heart, [0, 2*pi], [0;2]) and note that MATLAB dynamically plots both y1 and y2 versus t.

b) To change the output routine, enter options = odeset('OutputFcn', 'odephas2'), followed by ode45(@heart, [0, 2*pi], [0; 2], optionsxs). Note that MATLAB dynamically plots y2 versus y1 in the phase plane.