function dx = oscill_modele(t,x,flag,A,B,C)
dx=A*x(:)+B*oscill_entree(t,x)+C;