ch14ex7.m

% Chapter 14 Exercise 7.
t0=0; tf = 10;
w0 = [1;0];
[t,w]= ode23('pend',t0, tf,w0);
subplot(2,1,1);
plot(t,w(:,1));ylabel('th');
subplot(2,1,2);
plot(t,w(:,2));ylabel('th_dot');