1 PROGRAM NEWTON(INPUT,OUTPUT);
6 X,Y,C: ARRAY[0..NMAX] OF REAL;
7 A: ARRAY[-1..NMAX] OF REAL;
12 WRITELN("I N T E R P O L A T I O N S P O L Y N O M (NEWTON)");
13 WRITELN(" C0 + C1*X + C2*X^2 + .... + CN*X^N");
16 (***** EINGABE DER STUETZSTELLEN *****)
18 WRITELN("GEGEBENE PUNKTE:");
20 BEGIN READ(X[I],Y[I]);
21 WRITELN("P",I:1,X[I]:8:2," ,",Y[I]:8:2)
24 (***** BERECHNUNG DER DIVIDIERTEN DIFFERENZEN *****)
27 FOR J := N DOWNTO I DO Y[J] := (Y[J-1]-Y[J])/(X[J-I]-X[J]);
29 (***** KOMPOSITION DER KOEFFIZIENTEN *****)
34 A[I] := 1; C[I] := Y[I];
35 FOR J := I-1 DOWNTO 0 DO
37 A[J] := A[J-1] - A[J]*X[I-1];
38 C[J] := C[ J ] + A[J]*Y[ I ]
42 (***** AUSGABE DER KOEFFIZIENTEN *****)
45 WRITELN("DIE KOEFFIZIENTEN SIND:");
46 FOR I := 0 TO N DO WRITELN("C",I:1," =",C[I]:15:8)