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κ).