disp('Utlizza il file ode23, oppure ode23s, per il problema iniziale con funzione f data in riga 3')
close all
f=inline('-20*(x-sin(t))+cos(t)','t','x')
options=odeset('RelTol',1e-3,'AbsTol',1e-3);
[T,Y]=ode23s(f,[0,15],5,options);
%[T,Y]=ode23s(f,[0,15],5,options);
plot(T,Y,'o')
pause
disp('andamento dei passi')
d=T(2:end)-T(1:end-1);
plot(T(1:end-1),d)
pause
plot(d)