function bvpscripts
clear all
xlow=0;
xhigh=pi/2;
solinit = bvpinit(linspace(xlow,xhigh,10),[1 -1]);
sol = bvp4c(@bvp4ode,@bvp4bc,solinit);
xint = linspace(xlow,xhigh);
Sxint = deval(sol,xint);
plot(xint,Sxint(1,:))
xlabel('x')
ylabel('y')
% -----------------------------------------------
function dydx = bvp4ode(x,y,lambda)
dydx = [  y(2)
         1-y(1) ];
% -----------------------------------------------
function res = bvp4bc(ya,yb,lambda)
res = [  ya(1)-1
         yb(1) ];