Résolution équation différentielle (scilab)

Bonjour,

Je souhaiterais trouver une solution numérique à l'équation différentielle suivante :
c1.y''(t)+c2.y'(t)+c3.y(t)=sin(c4.t+c5.y(t))
y''(0)=0
y'(0)=0
y(0)=0
avec:
-> t, la variable
-> y(t), la fonction inconnue dont je recherche une solution numérique "convergente" sur l'intervalle [t1;t1+2.pi/c4]
-> c1, c2, c3, c4 (pulsation de l'excitateur), et c5 qui sont des paramètres constants

Je suis ingénieur, et j'ai besoin de cette solution pour une application concrète. Malheureusement, mes études commencent à dater et je ne sais pas du tout comment aborder un tel problème (l'ai-je déjà su !?) :
Mathématiquement, je ne sais pas comment gérer le sin(c4.t+c5.y(t))
Au niveau solveur, je dispose de scilab (que je connais à peine) et je ne sais pas quelle fonction utiliser.

Quelqu'un pourrait-il m'éclairer ? Me proposer une méthodologie, voir me proposer un script "scilab" pour résoudre mon problème ?

En vous remerciant par avance...