%% Script collapse: a csatolási paramétert növelve 2 bifurkáción át a redszer a fõátlóra omlik
figure
% adott pontokból induló tajektóriák kirajzolása különbözõ k értékeknél
%
%
for k=0:0.05:1 %csatolási paramméter változtatása
Collapse=(@(t, x)[x(1)-x(1)^3+k*(x(2)-x(1)); x(2)-x(2)^3+k*(x(1)-x(2))]); 
 
 for i=-2.1:0.2:2.1 %pontokból induló trajektóriák kiszámolása
    for j=-2.1:0.2:2.1
        x0=[i j];
        [t, x] = ode45(Collapse, [0 10], x0);
        plot(x(:,1),x(:,2));
        hold on;
    end
 end
    s=sprintf('k %1.4f',k);
    title(s);
    
 %pause(0.1);
 waitforbuttonpress;
 hold off;
 end
%}

%a vektrormezõs ábrázolása adott pontokban
%
step=0.1;
[x1, x2] = meshgrid(-1.2:step:1.2, -1.2:step:1.2);
for k=0:0.01:1
%    k
    x1dot_n=x1-x1.^3+k*(x2-x1); %adott pontokban vektormezõ iránya
    x2dot_n=x2-x2.^3+k*(x1-x2);
    x1dot=x1dot_n./sqrt(x1dot_n.^2+x2dot_n.^2); %normálva vektormezõ
    x2dot=x2dot_n./sqrt(x1dot_n.^2+x2dot_n.^2);
    
    quiver(x1,x2,x1dot, x2dot);
    s=sprintf('k %1.4f',k);
    title(s);
    xlim([-1.2 1.2]);
    ylim([-1.2 1.2]);
    waitforbuttonpress;
%    pause(0.01);
end
%}