Euler's Method for Initial Value Problems

function eulers_ivp_method(f::Function, T::S, t0::S, t1::S, h::S) where {S<:Real}
    collect(begin 
            T += h * f(T)
            T
        end for t in t0:h:t1)
end