disp('metodo shooting-secanti per il problema dy=f(t,y) con condizioni di periodicità')
y1=input('primo valore di tentativo y1=');
y2=input('secondo valore di tentativo y2=');
[T,Y]=ode23('fun',[0,1],y1);
plot(T,Y)
pause
hold on
F1=Y(end)-y1;
while 1
    [T,Z]=ode23('fun',[0,1],y2);
    plot(T,Z)
    pause
    F2=Z(end)-y2;
    y=y2-F2*(y2-y1)/(F2-F1);
    y1=y2;
    y2=y;
    F1=F2
end
    
