clear all

global gamma F

ua=[0;1];
us=[1;0];

x0=[0;0]+1.e-03*ua;
x0=[1;1]+[0.01;0]
dt=0.001;

Nt=6*5000

tt(1)=0;
xt=x0;
x(1)=xt(1);y(1)=xt(2);

for it=1:Nt
    tt(it)=(it-1)*dt;

    k1=dt*fdot2D(xt);
    k2=dt*fdot2D(xt+k1/2);
    k3=dt*fdot2D(xt+k2/2);
    k4=dt*fdot2D(xt+k3);

    xt=xt+1/6*(k1+2*k2+2*k3+k4);
    x(it)=xt(1);y(it)=xt(2);

end
%%
figure(34)
plot(x,y,'r','linewidth',2)
xlabel('$x$','interpreter','latex')
ylabel('$y$','interpreter','latex')
set(gca,'Fontsize',24)
hold on



