Oscillateur harmonique, euler-chromer méthode

Bonjour,
je cherche a résoudre un système d'oscillation simple (masse + ressort), mais j'ai un petit doute avec mon algorithme (methode euler-chromer)

Les oscillations obéissent a la loi suivante :
mx''̈ + bx' + kx = f (t) (x'=dx/dt, x''=dx'/dt)

que l'on peut séparer en deux équations en introduisant les constantes

v' + γx' + ω0^2 x = f (t)/m et v(t) = x' (γ=b/m, ω0^2=k/m)

et avec la methode euler chromer on a :

vn+1 = vn − (γ*vn + ω0^2* xn )∆t + (1/m)*fn* ∆t

et xn+1 = xn + vn+1* ∆t

Mon problème et au niveau de la définition de f(t)=A*cos(Ωt),

fn correspondant sera simplement fn=A*cos(Ω*∆t) ?

Merci de votre aide.
̇