% Краевая задача, пример 2   
 % dy1/dt= y2, dy2/dt=-y1, y1(0)=1, y2(pi)=-1
rp= @(t,y) [y(2);-y(1)];
bc= @(ya,yb) [ya(1)-1;yb(1)+1];
xinit= 0:pi/60:pi;
solinit = bvpinit(xinit,[1 -1]);
sol= bvp4c(rp,bc,solinit)
plot(sol.x,sol.y(1,:),sol.x,sol.y(2,:));

