   % pohyb ve zkrizenem elmag. poli
clear;
clc;
# vymazeme nastaveni vykreslovani
__gnuplot_raw__ "reset;";
# nastavuje opet vykresleni plnou carou
__gnuplot_set__ data style lines;
# obnoveni rozsahu os
__gnuplot_set__ xrange [:]
__gnuplot_set__ yrange [:]
% definice poèáteèních podmínek (v základních jednotkách SI)
x(1)=0;     	% pocatecni poloha x
y(1)=0;     	% pocatecni poloha y
t(1)=0;
vx=0;       	% pocatecni rychlost vx
vy=5;       	% pocatecni rychlost vy
mez=20;     	% mez integrace
h = mez/1000; 	% krok integrace
B=1;
E=.2*B;	% pomer el. a magn. pole 10,1,.2
q=1;
% -----------------------------------------------
% vlastní numerická integrace - RK metoda 2. øádu
% -----------------------------------------------
i=2;
while t<mez
k1x=q*vy*B;
k1y=q*(E-vx*B);
vpx=vx+k1x*2*h/3;
vpy=vy+k1y*2*h/3;
k2x=q*vpy*B;
k2y=q*(E-vpx*B);
x(i)=x(i-1)+vx*h+(k1x+k2x)*h^2/4;
y(i)=y(i-1)+vy*h+(k1y+k2y)*h^2/4;
vx=vx+(k1x+3*k2x)*h/4;
vy=vy+(k1y+3*k2y)*h/4;
t(i)=t(i-1)+h;
i=i+1;
endwhile  
% -----------------------------------------------
% vykreslení závislosti
clearplot;
axis('equal')
__gnuplot_set__ tics out
__gnuplot_set__ xlabel "x(m)"
__gnuplot_set__ ylabel "y(m)"
__gnuplot_set__ nokey
plot(x,y);
