function bvppractice
xlow=???;
xhigh=???;
solinit = bvpinit(linspace(xlow,xhigh,20),[1 0]);
sol = bvp4c(@bvp4ode,@bvp4bc,solinit);
xint = linspace(xlow,xhigh);
Sxint = deval(sol,xint);
eps=0.1;
analyt=(exp(xint/eps)-1)/(exp(1/eps)-1);
plot(xint,Sxint(1,:),xint,analyt,'r')
xlabel('x')
ylabel('y')
errorvec=abs((analyt-Sxint(1,:))./analyt)
figure, plot(xint,errorvec)
ylabel('error')
xlabel('x')
% -----------------------------------------------
function dydx = bvp4ode(x,y)
eps=0.1;
dydx = [  ???? ; ???? ];
% -----------------------------------------------
function res = bvp4bc(ya,yb)
res = [  ????  ;  ???? ];