REM Programa que eval£a la integral de la funci¢n fnf  en (a,b)
REM mediante la regla del trapecio y compara
REM con el resultado exacto

CLS
PRINT "        INTEGRAL MEDIANTE REGLA DEL TRAPECIO "
PRINT
DEF fnf (x) = EXP(x)
a = 0: b = 4
exacto = EXP(b) - EXP(a)
       
        REM 10 INPUT "Entrar N (par>=2) "; n
        REM IF n < 2 THEN STOP

FOR m = 2 TO 20  '****************************************
n = 2 ^ m
h = (b - a) / n
suma = fnf(a) / 2

FOR i = 1 TO n - 1
        x = a + i * h
        suma = suma + fnf(x)
NEXT i

suma = suma + fnf(b) / 2
integral = h * suma
diferencia = exacto - integral
PRINT USING "N=#######   Error=##.######"; n; diferencia
        REM GOTO 10
NEXT m   '***************************************************