% Crtamo vektorsko polje za sustav diferencijalnih jednadzbi
% x'(t) = f(x, y)
% y'(t) = g(x, y)

a = 2;
D = 0.6;
b = 1;
C = 0.9;

f = @(x, y) a*x - a*x*y/D; % funkcija f(x, y) = ax - axyD
g = @(x, y) -b*y + b*y*x/C; % funkcija g(x, y) = -by + bxy/D

[x,y] = meshgrid(0:0.25:5); % intervali za x i y

dx = f(x, y); % x-koordinata vektora
dy = g(x, y); % y-koordinata vektora

dy_norm = dy./sqrt(dx.^2 + dy.^2); % normalizirana y-koordinata vektora
dx_norm = dx./sqrt(dx.^2 + dy.^2); % normalizirana x-koordinata vektora 

quiver(x, y, dx_norm, dy_norm); % crtanje vektorskog polja