%
%    Boundary value problem due to Troesch
%    Plot of the solution
%
load troesch.dat
%
t=troesch(:,1);
x=troesch(:,2);
%
clf
plot(t,x,'b')
axis([-0.1 1.1 -0.1 1.1])
xlabel('BVP due to Troesch, \lambda = 5','FontSize',14)
%
hold
plot(0,0,'bo','MarkerSize',9)
plot(1,1,'bo','MarkerSize',9)
%
%print -depsc troesch
%