%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Implementation of Euler integration for a single step
function [tf, ft] = EulerSingleStep( ode_func, t, y0 )
    ti = t(1);
    tf = t(2);
    h = tf - ti;
    k1 = feval( ode_func, ti , y0 );
    ft = y0 + h * k1;
end