C         Symbolic_Testing.f
C      SIN PROGRAM.
       DOUBLE PRECISION FUNCTION SINE(X,E)
C      THIS COMPUTES SIN(X) TO ACCURACY E
       DOUBLE PRECISION E, TERM, SUM
       REAL X
       TERM = X
       SUM = X
       		DO 20 I = 3, 100, 2
       		TERM = TERM * X**2/(I*(I-1))
       		SUM = SUM + ((-1)**(I/2)) * TERM
       IF (DABS(TERM) .LT. E) GOTO 30
  20   CONTINUE
  30   SINE = SUM
       RETURN
       END