function ChemicalOscillation(parametri)

%
%VARIA x0,y0

global a b

a =parametri(1);
b=parametri(2);
tmax=parametri(3);
x0=parametri(4);
y0=parametri(5);



[tau,sol] = ode45(@eqs_oscill,[0 tmax],[x0,y0]);


phi1=zeros(size(tau));
phi2=zeros(size(tau));


phi1=sol(:,1);
phi2=sol(:,2);

%
%
%

    subplot(2,1,1)
    plot(tau,phi1,'LineWidth',2)	
    xlabel('Time','FontSize', 24)
    ylabel('x','FontSize', 24)
    
    subplot(2,1,2)
    plot(tau,phi2,'r','LineWidth',2)	
    xlabel('Time','FontSize', 24)
    ylabel('y','FontSize', 24)
    pause
    figure(2)
    plot(phi1,phi2,'LineWidth',2)

function d_phi= eqs_oscill(tau,phi)


global a b

%
% Inizializzazione delle variabili.
%
d_phi = zeros(2,1);

%
% Sistema di equazioni differenziali.
%
d_phi(1)=a-phi(1)-4*phi(1)*phi(2)/(1+phi(1)^2); 
d_phi(2)= b*phi(1)*(1-phi(2)/(1+phi(1)^2));