2 PROGRAM INTEGRATION(INPUT,OUTPUT);
4 VAR ANFANG,ENDE,FLAECHE: REAL;
7 (************* HIER INTEGRAND DEFINIEREN **************)
8 FUNCTION F(X: REAL): REAL;
12 (******************************************************)
16 FUNCTION INTEGRAL(A,B: REAL): REAL;
20 DELTA, MITTE, EINF, ZWEIF, VIERF: REAL;
21 SIMPSON, LETZTWERT: REAL;
24 N := 1; DELTA := B-A; MITTE := A+DELTA/2;
28 SIMPSON := (EINF + 4*VIERF) * DELTA/6;
32 N := 2*N; DELTA := DELTA/2; MITTE := A + DELTA/2;
33 ZWEIF := ZWEIF + VIERF;
35 FOR J := 0 TO N-1 DO VIERF := VIERF + F(MITTE+J*DELTA);
36 SIMPSON := (EINF + 2*ZWEIF + 4*VIERF) * DELTA/6
37 UNTIL ABS(SIMPSON-LETZTWERT)<ABS(SIMPSON)*EPSILON;
44 WRITELN("B E S T I M M T E S I N T E G R A L");
46 WRITE("INTERVALLGRENZEN"); READ(ANFANG,ENDE);
47 FLAECHE := INTEGRAL(ANFANG,ENDE);
49 WRITELN("WERT DES INTEGRALS: ", FLAECHE :18:10)