1. Choose either a time step (∆t) or a positive integer (N) such that ∆t = T/N and tκ = k∆t, for k = 0,1,2,...,N.
2. Let U₀ = y(0) = A.
3. For k = 0,1,2,...,N-1, compute uvk+1 = uvk + f(tκ,uκ)∆t.
Each uvk is an approximation to the exact solution y(tκ).