function dy=quadr(t,y)

global alfa
%alfa=1;     % transcritica
%alfa=2.95;  % estacionario
%alfa=3;     % hopf
%alfa=3.2;   % periodo-2
%alfa=3.5;   % periodo-4
%alfa=3.75;  % caos
%alfa=3.83;  % periodo-3 apos t=60

dy = (alfa*(1-y)-1)*y;

%J=(a-1)-2*ay --> J(0)=a-1, J(y*)=1-a
%y*=(a-1)/a
%a=1 --> bifurcacao transcritica

%y0=[0.1]; n=60; [t,y]=eulere('quadr',0,n,y0,n); plot(t,y); % caos pelo efeito da discretizaca
%y0=[0.1]; n=60; [t,y]=euleri('quadr',0,n,y0,n); plot(t,y); % nao ha´ sol. periodicas

%sensibilidade a condicao inicial (testar para varios alfa):
%y0=[0.1]; [t,y]=eulere('quadr',0,160,y0,160); plot(t,y); hold on
%y0=[0.101]; [t,y]=eulere('quadr',0,160,y0,160); plot(t,y,'g');
%y0=[0.10001]; [t,y]=eulere('quadr',0,160,y0,160); plot(t,y,'m');

%Diagrama cobweb
%m=length(y); p=m-1;
%A=[zeros(m,1) eye(p,m)'];
%xk=y';
%xk1=A*xk;
%stairs(xk(1:p),xk1(1:p))
