%% 2. gyakorló feladat
% Elsőrendű, kétváltozós diffegyenlet megoldása
fs = filesep;
file_path = ['.', fs, 'lab10_pelda2_adat.bin'];
fid = fopen(file_path);
d = fread(fid, 6, 'double');
fclose(fid);
p = d(1);
q = d(2);
y0 = [d(3), d(4)];
tspan = d(5:6);
modell = @(t,y) [(1-y(2)/q)*y(1); -(1-y(1)/p)*y(2)];
[T, Y] = ode45(modell, tspan, y0);
plot(T, Y(:,1), 'ks-')
hold on
plot(T, Y(:,2), 'rs-')
ind1 = find(Y(:,2) < Y(:,1), 1);
ind2 = find(Y(ind1:end,1) < Y(ind1:end,2), 1);