Step 3: Approximate derivative(s) by finite difference approximation(s). Very common (standard!) formula for u′′: u′′(tn)≈un+1−2un+un−1Δt2
Use this discrete initial condition together with the ODE at t=0 to eliminate u−1: un+1−2un+un−1Δt2=−ω2un