function [ t, u ] = ode_euler ( f, a, b, ua, n )

  dt = ( b - a ) / n;
  t = linspace ( a, b, n + 1 );
  u = zeros ( n + 1, 1 );
  
  u(1) = ua;  
  for i = 1 : n
    u(i+1) = u(i) + dt * f(t(i),u(i));
  end

  return
end

