clear;
x=[[-5:0.5:-0.5],[0.5:0.5:5]]; % vynecháváme x=0 a y=0, protože tam složky
                               % vektoru síly nabývají nekonečných hodnot
y=x;

[X,Y]=meshgrid(x,y); % z vektorů x a y uděláme matice X a Y
Fx=-(Y./X.^2+1./Y);
Fy=(X./Y.^2+1./X);
% definice složek Fx a Fy vektorového pole - všimněte si,
% že veškeré operace s maticemi provádíme s operátory s tečkou (.* .^ ./)
% tečka značí tzv elementwise operaci, tedy početní operaci mezi
% jednotlivými prvky matic - [a,b].*[c,d]=[a*c,b*d]

quiver(x,y,Fx,Fy,2,'b');
grid;
fig=gcf(); fig.Position=[50 50 900 900];
xlabel('$x$','FontSize',16,'interpreter','latex');
ylabel('$y$','FontSize',16,'interpreter','latex');
daspect([1,1,1]); %nastavení stejného měřítka souřadnicových os
mtit=['$' ...
    '\vec{F}=' ...
    '-\frac{F_0}{R}\left(\frac{y}{x^2}+\frac{1}{y}\right)\:\vec{i} +' ...
    '\frac{F_0}{R}\left(\frac{x}{y^2}+\frac{1}{x}\right)\:\vec{j}' ...
    '$'];
h=title(mtit,'interpreter','latex'); h.FontSize=18;