ex13ch2.m

these codes are for solving OED with matlab

function ex13ch2
%options = odeset('Stats','on');
options = odeset('Stats','on','RelTol',1.e-8,'AbsTol',1.e-8);
tout = linspace(2,10,9);
tspan = [1 tout];
y0 = 1;
[t1,y1] = ode45(@ex13ch2f,tspan,y0,options);
disp(t1)
disp(y1)
plot(t1,y1)
legend('Solution of dy/dt = -1/t^2 + 10(y-1/t) by ode45')
[t2,y2] = ode15s(@ex13ch2f,tspan,y0,options);
disp(t2)
disp(y2)
figure
plot(t2,y2)
legend('Solution of dy/dt = -1/t^2 + 10(y-1/t) by ode15s')
%===========================================================
function yprime = ex13ch2f(t, y, prob)
yprime = -1 / t^2 + 10 * (y - 1 / t);