disp('metodo shooting -secanti- per il problema ddy=f(t,y)')
y0=input('valore iniziale y0=');
yf=input('valore finale yf=');
d1=input('prima derivata di tentativo=');
d2=input('seconda derivata di tentativo=');
[T,Y]=ode23('fun2',[0,1],[y0,d1]);
plot(T,Y(:,1))
pause
hold on
F1=Y(end,1)-yf;
veccolor=['b','r','g','m','k']; k=0;
while 1
    [T,Z]=ode23('fun2',[0,1],[y0,d2]);
    plot(T,Z(:,1),veccolor(rem(k,5)+1))
    k=k+1;
    pause
    F2=Z(end,1)-yf;
    d=d2-F2*(d2-d1)/(F2-F1);
    d1=d2;
    d2=d;
    F1=F2
end
    
