C		Routine de definition des equations differentielles de
C	    l'oscillateur amorti par un frottement fluide.
C		Dominique Lefebvre avril 2006
	  SUBROUTINE Derivee(x,y,dy, NbEquation)
	  
	  IMPLICIT None
		
	  REAL alpha, beta
	  
C	  alpha = omega0/Q
C	  beta = omega0*omega0
	  COMMON /Oscillateur/alpha, beta
	
	  INTEGER nbEquation
	  REAL x, y(NbEquation), dy(NbEquation)

C		Description du systeme differentiel de l'oscillateur
C		y(1) = y
C		y(2) = dy/dt
		  dy(1) = y(2)						! dy/dt
		  dy(2) = -alpha*y(2) - beta*y(1) 	! d2y/dt2

	  RETURN
	  END