function budworm(parametri)

%budworm([0.8,20,0.35]) VARIA x0.

global r k  

r =parametri(1);
k=parametri(2);
x0=parametri(3);

ttot=120

[tau,sol] = ode45(@budworm1,[0 ttot],[x0]);


phi=zeros(size(tau));


phi=sol(:,1);


%
%
%

    figure(2)
    plot(tau,phi,'g','LineWidth',2)
	
    xlabel('Time','FontSize', 24)
ylabel('x(t)','FontSize', 24)
 

function d_phi= budworm1(tau,phi)


global r k 

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

%
% Sistema di equazioni differenziali.
%
d_phi(1)= r*phi*(1-phi/k)-phi^2/(1+phi^2);