return column vector (t(i+1) & y(i+1)) instead of row vector

function [y,t] = euler_timestep(f,t0,tf,alpha,N)
h=(tf-t0)/N; 
t(1)=t0;    
y(1)=alpha; 
for i=1:N
    t(i+1)=t(i)+h; 
    y(i+1)=y(i)+h*f(t(i),y(i));  
end

Answer 1

[Walter Roberson]

t = zeros(N+1,1);
y = zeros(N+1,1);