home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / fortran / library / ssp / polyops / pint.for < prev    next >
Text File  |  1985-11-29  |  1KB  |  47 lines
  1. C
  2. C     ..................................................................
  3. C
  4. C        SUBROUTINE PINT
  5. C
  6. C        PURPOSE
  7. C           FIND INTEGRAL OF A POLYNOMIAL WITH CONSTANT OF INTEGRATION
  8. C           EQUAL TO ZERO
  9. C
  10. C        USAGE
  11. C           CALL PINT(Y,IDIMY,X,IDIMX)
  12. C
  13. C        DESCRIPTION OF PARAMETERS
  14. C           Y     - VECTOR OF COEFFICIENTS FOR INTEGRAL, ORDERED FROM
  15. C                   SMALLEST TO LARGEST POWER
  16. C           IDIMY - DIMENSION OF Y (EQUAL TO IDIMX+1)
  17. C           X     - VECTOR OF COEFFICIENTS FOR ORIGINAL POLYNOMIAL,
  18. C                   ORDERED 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           DIMENSION OF Y IS SET AT DIMENSION OF X PLUS ONE, AND THE
  29. C           CONSTANT TERM IS SET TO ZERO. INTEGRAL IS THEN CALCULATED
  30. C           BY DIVIDING COEFFICIENTS BY THEIR RESPECTIVE EXPONENTS.
  31. C
  32. C     ..................................................................
  33. C
  34.       SUBROUTINE PINT(Y,IDIMY,X,IDIMX)
  35.       DIMENSION X(1),Y(1)
  36. C
  37.       IDIMY=IDIMX+1
  38.       Y(1)=0.
  39.       IF(IDIMX)1,1,2
  40.     1 RETURN
  41.     2 EXPT=1.
  42.       DO 3 I=2,IDIMY
  43.       Y(I)=X(I-1)/EXPT
  44.     3 EXPT=EXPT+1.
  45.       GO TO 1
  46.       END
  47.