%solve a 1d dif equation
%xdot = a x
clear all

delta=0.01;
N=1001;

a=1;
x0=2;

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

for it=2:N

    t(it)=(it-1)*delta;
    %v=a*x(it-1);
     v=x(it-1)*(1-x(it-1));

    x(it)=x(it-1)+delta*v;
%figure(100)
% plot(t(it),x(it),'*b','Linewidth',2,'Markersize',20)
% xlabel('$t$','Interpreter','latex')
% ylabel('$x(t)$','Interpreter','latex')
% set(gca,'Fontsize',24)
% grid on
% drawnow
% hold on
% %pause(0.2)


end
%figure(100);hold off

figure(102)
%plot(t,x,'b',t,exp(a*t)*x0,'or','Linewidth',2)
plot(t,x,'','Linewidth',2)
xlabel('$t$','Interpreter','latex')
ylabel('$x$','Interpreter','latex')
set(gca,'Fontsize',24)
grid on
hold on
% fan=exp(a*t)*x0;
% 
% abs(x(end)-fan(end))/fan(end)


