ex25ch2.m

these codes are for solving OED with matlab

function ex25ch2
options = odeset('Events',@event,'RelTol',1e-5,'AbsTol',1e-10);
[t,y,te,ye,ie] = ode45(@odes,[0 20*pi],[0; 1],options);
plot(t,y(:,1),te,ye(:,1),'ro')
spacing = diff(te)
%=================================================
function dydt = odes(t,y)
dydt = [y(2); -sin(y(1))];

function [value,isterminal,direction] = event(t,y)
isterminal = 0;
direction = 0;
value = y(1);