function [ y ] = forward_euler(dt,y0,t)
%FORWARD EULER
N=length(t);
y=zeros(N,1); %Sólo tienes N elementos por eso el bucle termina en N-1
y(1)=y0;

for k=1:N-1

    f=my_function(t(k),y(k));
    y(k+1)=y(k)+dt*f;
    
end