465 2

Rozdział 4

WRJTE(6, 901)

I2=4*(M —1)+1 DO 10 1=1, 12

XARG = -1 .DO - DFLOAT(I- 1)/RM D=P(XARG)— F(XARG)

10 WRITE(6. 902) XARG, D STOP

900    FORMAT(I5)

901    FORMAT(4X, 'X', 5X_? T(X)-F(X)')

902    FORMAT(F8.3, Dl2.3)

END

C

C OBLICZENIE WSPÓŁCZYNNIKÓW Cl, C2, ...,CM WIELOMIANU C INTERPOLACYJNEGO

SUBROUTINE INTERP(X, Y, C_; N)

DOUBLE PRECISION C(20), Cl, PI, SfGMA, X(20). XI, Y(20)

Cl-MD

C(1)-C1

C(21-(Y(2) —C1)/(X(2)—X(l)>

DO 100 1 = 3, N PI=1.E>0 SIGMA=0.D0 J2 = I-J XI=X(I)

DO 50 J = 2, J2 PI = PI*<XI-X<J-1))

50 SIGMA = SIGM A ł C(J)*P1 100 C<n=(Y(l)-Cl-SrGMA)/(P[*(Xr-X(I- I)))

RETURN

END

C

■C OBLICZENIE WARTOŚCI FUNKCJI F DLA XARG ZA POMOCĄ C WH-LOMIANU INTERPOLACYJNEGO REAL FUNCTION P*8(ARG)

DOUBLE PRECISION ARG, C(20), SIGMA, X(20), Y(20) COMMON C, X, Y. M SIGMA=C(M)

12 = M — 1 DO 100 1 = 1, 12

100 SIGM A = (ARG - X(M - i))*SIGM A + C(M -1)

P-SIGMA

RETURN

END