%% Script2
close;
clc;
clear;

t0=[0 10];
x0=[1 1]; %1 vonal

figure
%{
[t, y] = ode45('lvspec',t0,x0);
plot(y(:,1),y(:,2))
xlim([0 7]);
ylim([0 7]);
%hold on
%}

%{
%sok vonal
X=[1 1; 0.01 0.2; 5 5; 6 1;...
    1 6; 0.1 3; 0.5 2; 0.3 1.7; 0.1 1.9; ...
    0.03 0.5; 0.03 1;0.03 1.5; 0.03 1.8];
for i=1:size(X,1)
    x0=X(i,:);
    [t, y] = ode45('lvspec',t0,x0);
    plot(y(:,1),y(:,2))
    xlim([0 7]);
    ylim([0 7]);
    hold on    
end
%}

%% Keressük meg a két egyenest 
%
%1.: y=(3-x)/2
%2.: y=2-4x
%
xe1=0:0.01:3;
xe2=0:0.01:0.5;
ye1=(repmat(3,1,numel(xe1))-xe1)/2;
ye2=repmat(2,1,numel(xe2))-4*xe2;

hold on
%figure
plot(xe1,ye1);
hold on
plot(xe2,ye2);
xlim([0 7]);
ylim([0 7]);
%} 

% %% Viszonylag meredek szakaszra van szükség
% sok szakaszt rajzol 2 egyenesre
%
N=100;
sx0=0;
sy0=1;
sx1=0.2;
sy1=0;
sx=linspace(sx0,sx1,N);
sy=linspace(sy0,sy1,N);

%figure
for i=1:N
    x0=[sx(i) sy(i)];
    [t, y] = ode45('lvspec',t0,x0);
    plot(y(:,1),y(:,2))
    xlim([0 7]);
    ylim([0 7]);
    hold on   
end

hold on
%figure
plot(xe1,ye1);
hold on
plot(xe2,ye2);
xlim([0 7]);
ylim([0 7]);
hold on
%% Látszik, hogy nagyon laposan mennek be a trajektóriák
%}

N=500;
sx0=0;
sy0=4;
sx1=4;
sy1=0;
sx=linspace(sx0,sx1,N);
sy=linspace(sy0,sy1,N);

%figure
for i=1:N
    x0=[sx(i) sy(i)];
    [t, y] = ode45('lvspec',t0,x0);
    plot(y(:,1),y(:,2))
    xlim([0 7]);
    ylim([0 7]);
    hold on   
end

hold on
%figure
plot(xe1,ye1);
hold on
plot(xe2,ye2);
xlim([0 7]);
ylim([0 7]);
hold on

%{
N=80;
t0=[0 4];
S2=[repmat(2.8,N,1) linspace(0.007,0.005,N)'];
for i=1:N
    x0=S2(i,:);
    [t, y] = ode113('lvspec_vissza',t0,x0);
    plot(y(:,1),y(:,2),'r')
    xlim([0 7]);
    ylim([0 7]);
    hold on   
end
%}









