File: POWER.FT of Disk: Disks/MyPDP/m8-backup-rka1-rkb1
(Source file text)
SUBROUTINE POWER(IP,X,XK,F,P) C INTERPOLATION MITTELS POLYNOM DIMENSION XK(20) DIMENSION F(20) SUM=0. DO 10 K=1,IP PROD=1. DO 20 M=1,IP IF (M-K) 15,20,15 15 PROD=PROD*(X-XK(M))/(XK(K)-XK(M)) 20 CONTINUE SUM=SUM+F(K)*PROD 10 CONTINUE P=SUM RETURN END