home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / fortran / library / ssp / polyops / pqsd.for < prev    next >
Text File  |  1985-11-29  |  1KB  |  47 lines
  1. C
  2. C     ..................................................................
  3. C
  4. C        SUBROUTINE PQSD
  5. C
  6. C        PURPOSE
  7. C           PERFORM QUADRATIC SYNTHETIC DIVISION
  8. C
  9. C        USAGE
  10. C           CALL PQSD(A,B,P,Q,X,IDIMX)
  11. C
  12. C        DESCRIPTION OF PARAMETERS
  13. C           A     - COEFFICIENT OF Z IN REMAINDER (CALCULATED)
  14. C           B     - CONSTANT TERM IN REMAINDER (CALCULATED)
  15. C           P     - COEFFICIENT OF Z IN QUADRATIC POLYNOMIAL
  16. C           Q     - CONSTANT TERM IN QUADRATIC POLYNOMIAL
  17. C           X     - COEFFICIENT VECTOR FOR GIVEN POLYNOMIAL, ORDERED
  18. C                   FROM SMALLEST TO LARGEST POWER
  19. C           IDIMX - DIMENSION OF X
  20. C
  21. C        REMARKS
  22. C           NONE
  23. C
  24. C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
  25. C           NONE
  26. C
  27. C        METHOD
  28. C           POLYNOMIAL IS DIVIDED BY THE QUADRATIC Z**2-P*Z-Q GIVING
  29. C           THE LINEAR REMAINDER A*Z+B
  30. C
  31. C     ..................................................................
  32. C
  33.       SUBROUTINE PQSD(A,B,P,Q,X,IDIMX)
  34.       DIMENSION X(1)
  35. C
  36.       A=0.
  37.       B=0.
  38.       J=IDIMX
  39.     1 IF(J)3,3,2
  40.     2 Z=P*A+B
  41.       B=Q*A+X(J)
  42.       A=Z
  43.       J=J-1
  44.       GO TO 1
  45.     3 RETURN
  46.       END
  47.