function zadacha

% Инициализиране на матрицата А
A=[4,1; -3, 4]

% Вектор стълб със свободните членове
rhs=[5;1]

% Равновесната точка
a=A\(-rhs)

% Намиране на собствените стойности в D
D=eig(A)

x=a(1)-4:a(1)+4
y=a(2)-4:a(2)+4

hold on

axis([a(1)-7,a(1)+7,a(2)-9,a(2)+9])

% Функция за дясната страна
function z= ff(t,y)
    z=A*y+rhs
end

% фазови криви
[X,Y]=meshgrid(x,y);
tmax=50;
for i=1:length(x)
    for j=1:length(y)
        [U,Z]=ode45(@ff,[0,tmax],[X(i,j),Y(i,j)]);
        [U1,Z1]=ode45(@ff,[0,-tmax],[X(i,j),Y(i,j)]);
        plot(Z(:,1),Z(:,2),Z1(:,1),Z1(:,2),'b');
    end
end

% тангенциални вектори
% P и Q са  вектори v(f(x,y),g(x,y)
P=A(1,1)*X+A(1,2)*Y+rhs(1);
Q=A(2,1)*X+A(2,2)*Y+rhs(2);
D=sqrt(P .^2+Q.^2);
quiver(X,Y ,P ./D,Q./D,0.5,'r')

plot(a(1),a(2),'m*')

end