clear, clc, close all;
y00=0;
y01=1;
[x,U]=ode45(@fEXS715,[0 1],[y00;0] );
plot(x,U(:,1))
pause
hold on
yend0=U(end,1)-1;
yend1=yend0;
while abs(yend1)>1e-5    
   [x,U]=ode45(@fEXS715,[0 1],[y01;0] );
   plot(x,U(:,1))
   pause
   yend1=U(end,1)-1;
   y02=y01-(y01-y00)/(yend1-yend0)*yend1;
   y00=y01;
   y01=y02;
   yend0=yend1;
end
y01
