%Γομβερτζ
clear all

dt=2.7;

Tfin=100500;

Nt=ceil(Tfin/dt);

x0=0.1;

xt=x0;
x(1)=xt;
t(1)=0;

for it=2:Nt
    t(it)=(it-1)*dt;
    dotx=vx(x(it-1));
    x(it)=x(it-1)+dt*dotx;
end

xan=exp(log(x0)*exp(-t));

figure(101)
plot(t,x,'b',t,x,'or','linewidth',2)
set(gca,'Fontsize',24)

figure(301)
hist(x)




%figure(10)
%plot(t,x,'r',t,xan,'b','linewidth',2)
%set(gca,'Fontsize',24)







