d

function dy = lotka(t, y)

%zmienne
a = 1.3;
b = 1.5;
c = 1.2;
d = 1.6;

%set up differential equations
dy = zeros(2,1);
%dx po dt
dy(1) = a * y(1) - c * y(1) * y(2);
%dy po dt
dy(2) = b * y(2) - d * y(1) * y(2);
%end of function
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%execution

clc


%1st way of
ode = @(t,y)lotka(t,y)
[t,y] = ode45 (ode, tspan, y0)
plot (t,y)

%2nd way of
tspan = [0, 10]
y0 = [10;10]
ode45(@lotka, tspan, y0)