home *** CD-ROM | disk | FTP | other *** search
/ Programmer 7500 / MAX_PROGRAMMERS.iso / CLIPPER / MISC / LAGRINT.ZIP / LAGINT.FOR
Encoding:
Text File  |  1989-06-19  |  1.2 KB  |  40 lines
  1.       SUBROUTINE LAGINT (X,Y,N,XINT,YOUT)
  2. C
  3. C   FROM APPLIED NUMERICAL ANALYSIS, THIRD EDITION (GERALD / WHEATLEY)
  4. C   CHAPTER 3, PROGRAM 1
  5. C
  6. C  --------------------------------------------------------------
  7. C
  8. C   THIS SUBROUTINE PERFORMS LAGRANGIAN INTERPOLATION WITHIN A SET
  9. C   OF (X,Y) PAIRS TO GIVE THE Y VALUE CORRESPONDING TO XINT.  THE
  10. C   DEGREE OF THE INTERPOLATING POLYNOMIAL IS ONE LESS THAN THE
  11. C   NUMBER OF POINTS SUPPLIED.
  12. C
  13. C  --------------------------------------------------------------
  14. C
  15. C   PARAMETERS ARE:
  16. C   X      ARRAY OF VALUES OF THE INDEPENDENT VARIABLE
  17. C   Y      ARRAY OF FUNCTION VALUES CORRESPONDING TO X
  18. C   N      NUMBER OF POINTS
  19. C   XINT   THE X-VALUE FOR WHICH ESTIMATE OF Y IS DESIRED
  20. C   YOUT   THE Y VALUE RETURNED TO CALLER
  21. C
  22. C  --------------------------------------------------------------
  23. C
  24.       DIMENSION X(N),Y(N)
  25. C
  26. C  --------------------------------------------------------------
  27. C
  28.       YOUT=0.0
  29.       DO 20 I=1,N
  30.          TERM = Y(I)
  31.          DO 10 J= 1,N
  32.             IF (I.NE.J) THEN
  33.                TERM = TERM * (XINT - X(J))/(X(I)-X(J))
  34.             END IF
  35.    10    CONTINUE
  36.          YOUT = YOUT + TERM
  37.    20 CONTINUE
  38.       RETURN
  39.       END
  40.