ex22ch3.m

these codes are for solving OED with matlab

function ex22ch3
clf reset
hold on
for i = 1:3
    infinity = 5*i;
    if i == 1
        solinit = bvpinit(linspace(0,infinity,5),[pi 0]);
    else
 		solinit = bvpinit(sol,[0,infinity]);
    end
	sol = bvp4c(@odes,@bcs,solinit);
	plot(sol.x,sol.y(1,:),sol.x(end),sol.y(1,end),'ro');
	axis([0 15 0 4]);
	drawnow
end
hold off
fprintf('Initial slope of 2 is approximated by %g.\n',sol.y(2,1));
%===========================================================
function dydz = odes(z,y)
dydz =  [ y(2); -sin(y(1)) ];

function res = bcs(ya,yb)
res = [ ya(1); yb(1)-pi ];