467 2

467

Rozdział 4

C(1)=0.

S(J)=YP/PP0 P(2, 1)-l.

P(2. 2) = B(1)

DO 70 1 = 3, N IM1=1— J IM2-1-2 PP=0.

PXP = 0.

YP=0.

DO 50 J = l. M PJ=0.

XT= I.

DO 30 K=2_y 1M1 XT=XT*X(J)

30 XJ(K) = XT

DO 40 K = l, FM1

40 PJ=PJ+P(IMI, K)*XJ(IM1 —K +1)

PP=PP + PJ*PJ PXP=PXP-4-X(J>PJ*PJ 50 YP=YP + Y(J)* PJ B(3 M1) = — PXP/PP C(IM1) = PP/PP0 S(rMl)=YP/PP PP0=PP P(I. 1)-L

P(I, 2)=B(1M1)+P0M1, 2)

IF (L EQ. 3) GOTO 70 DO 60 J=3, £M1

60 P(I, J) =    J) + B(IM1)*P(IMI, J — 1)-C(IMI)*P(IM2, J-2)

70 P(I, I) = B(IM1)*P(IMI, IMl)-C(rMl)*P(rM2, IM2)

PP=0.

YP=0.

DO 100 J=l,

PJ=0.

XT = i.

DO 80 K=2. N XT=XT *X(J)

80 XJ(K)=XT DO 90 K = 1, N

90 PJ=PJ-f P(K, K)*XJ(N—K-f 1)

PP=PP + PJ*PJ 100 YP = YP+Y(J)*PJ