To use this m-file to generate Figure M5.7a, enter the following at the MATLAB command window prompt:
x0 = [124.6525;154.9880]; [t,x] = ode45('heater',0,5,x0); % MATLAB 4.x [t,x] = ode45('heater',[0 5],x0); % MATLAB 5.x