clear;
x=[[-5:0.5:-0.5],[0.5:0.5:5]];
y=x;
[X,Y]=meshgrid(x,y);

gda().font_size=5; gdf().figure_position(2)=0;
gcf().figure_size=[1210,1210];
Fx=-(Y./X.^2+1./Y);
Fy=(X./Y.^2+1./X);

champ(x,y,Fx',Fy'); // POZOR aby to spravne fungovalo, matice Fx a Fy MUSI byt transponovane!!!!!!!!!!!!!!
xgrid;
xlabel('x','fontsize',5); ylabel('y','fontsize',5);
title("$\ \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} $",'fontsize',6)
gce().colored='on';
gca().isoview='on';
gcf().color_map=jetcolormap(64);
colorbar(0,max((Fx.^2+Fy.^2).^(1/2)));
gce().y_label.text=["velikost síly [N]","$\ \times F_0/R $"];
gce().y_label.font_size=6;