% Crtamo vektorsko polje za diferencijalnu jednadzbu
% y'(x) = f(x,y)

f = @(x, y) y - x; % funkcija f(x, y) = y - x

[x, y] = meshgrid(-2:0.5:2); % intervali za x i y

dy = f(x, y); % y-koordinata vektora
dx = ones(size(dy)); % x-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