1 % Demo: Regla parabolica
2 % Integral(f,a,b)
3 clc, clear all
4 disp('Demo: Regla parabolica')
5 disp('======================')
6 % Integrando
7 f = @(x) exp(x^2)
8 % Limites de integracion
9 a = 0;
10 b = 1;
11 % main
12 n =100;
13 h = (b-a)/n;
14 sumai =0; % suma impares
15 sumap =0; % suma pares
16 %
17 for i = [Link]n-1
18 x = a+ i*h;
19 sumai = sumai + f(x);
20 endfor
21 %
22 for i = [Link]n-2
23 x = a + i*h;
24 sumap = sumap + f(x);
25 endfor
26
27 Inte = (h/3)*(f(a) + 4*sumai + 2*sumap + f(b));
28
29 fprintf(' Integral = %7.2f \n', Inte)
30
31
1 % Demo: Regla de los Trapecios
2 clc, clear all
3 disp('Demo: Regla de los trapecios')
4 disp('============================')
5 % Extremos :
6 a = 0; % extremo inferior
7 b = pi/2; % extremo superior
8 % Integrando :
9 f = @(x) sin(x^2)
10 % main
11 n = 100;
12 h = (b-a)/n;
13 sum = 0;
14 for i = 1:n-1
15 x = a + i*h;
16 sum = sum + f(x);
17 endfor
18 Inte = (h/2)*(f(a) + 2*sum + f(b));
19
20 fprintf(' Integral = %7.2f \n', Inte)
21
22
23
24
25
1 % Demo: Regla Parabolica
2 % Integral(f,a,b)
3 clc, clear all
4 disp('Demo: Regla Parabolica')
5 disp('======================')
6 % Data:
7 a = 0; % extremo inferior
8 b = 20; % extremo superior
9 h = 2; % tamaño de paso
10 x = a:h:b; % partición del intervalo [a,b]
11 n = 10; % número de partes [a,b]
12 f = [0 39 63 128 117 106 94 79 58 18 0 ];
13 % Data:
14 Data = [x' f']
15 % Main
16 sumai =0; % suma impares
17 sumap =0; % suma pares
18 %
19 for i = [Link]n-1
20 x = x(1) + i*h;
21 sumai = sumai + f(i);
22 endfor
23 %
24 for i = [Link]n-2
25 x = x(1) + i*h;
26 sumap = sumap + f(i);
27 endfor
28
29 Inte = (h/3)*(f(1) + 4*sumai + 2*sumap + f(n));
30
31 fprintf(' Integral = %7.2f \n', Inte)
32
33
