%a=inline('[0,-1;1,0]');      %la funzione a(x)
%b=inline('zeros(2,1)');       %la funzione b(x)
%t0=0;              %punto iniziale
%tf=5;           %punto finale
%y0=[1,0]';             %valore iniziale



%oscillante
%  a=inline('1');
%  b=inline('-.5*sin(50*x)-sin(x)+25*cos(50*x)+cos(x)');
% solution='.5*sin(50*t)+sin(t)';    

%robinson
%a=inline('-20');
% b=inline('-20*sin(x)+cos(x)');


a=inline('sqrt(-1)');
b=inline('0');
t0=0;              %punto iniziale
tf=20;           %punto finale
y0=1;           %valore iniziale
h=.2;               %passo iniziale
