%% Script collapse
clear
clc
figure

%{
%pontokból indítom a rendszert zh-k körül adott csatolás erõvel, trajekt néz
k=10; %0.33
xv=-2:0.2:2;
[X Y] = meshgrid(xv,xv);
plot(X,Y,'rx');
hold on;

for i=xv
    for j=xv
        x0=[i j];
        [t, x] = ode45(@(t,x) Collapse(t,x,k),[0 100], x0);
        hold on
        plot(x(:,1),x(:,2));
    end
end
%}

%
%a csatolás erõsödésével változik a rsz, erõs csatolásnál csak a fõátlón maradnak zh-k
%vektormezõs ábra
step=0.1;
[x1, x2] = meshgrid(-1.2:step:1.2, -1.2:step:1.2);
for k=0:0.01:10
%    k
    x1dot=x1-x1.^3+k*(x2-x1);
    x2dot=x2-x2.^3+k*(x1-x2);

    quiver(x1,x2,x1dot, x2dot, 1.5);
    s=sprintf('k %1.4f',k);
    title(s);
    xlim([-1.2 1.2]);
    ylim([-1.2 1.2]);
   waitforbuttonpress;
    % pause(0.01);
end
%}