function dy=BVPeqn(t,y)
%Sample BVP equations:

dy=zeros(size(y));
dy(1)=y(2);
dy(2)=4*y(1);