function [x, t, dt]=Euler_1ODE(t_in, t_fin, N, x0, f)

t=linspace(t_in,t_fin,N); dt=t(2)-t(1);
x(1)=x0;
for i=1:N-1
    x(i+1)=x(i)+f(x(i),t(i))*dt;
end

end