File: POWER.FT of Disk: Disks/MyPDP/m8-1-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

---