Utilice programas MATLAB para escribir fórmulas trapezoidales complejas y fórmulas complejas simples para resolver integrales.
Función y = traint(a, b, n, f)
h = (b-a)/n;
x = linspace(a, b, n+1);
y1 = h * feval(f, x);
y 1(1)= y 1(1)/2 ;
y 1(n+1)= y 1(n+1)/2;
y = suma(y 1);
% compuesto Fórmula de Simpson compuesta Fórmula de Simpson
Función y = sraint(a, b, n, f)
h = (b-a)/n;
x = linspace (a, b, 2 * n+1);
y1 = feval(f, x
y 1(2:2:2 * n)= 4 * y); 1(2:2:2 * n);
y 1(3:2:2 * n-1)= 2 * y 1(3:2:2 * n-1); p> p>
y = h/6 * suma(y 1);
& gt& gtf = inline(' 4./(1+x . * x)');
& gt& gttraint(-1, 1, 8, f)
ans =
6.2624
& gt& gtsraint(-1, 1, 8 , f)
ans =
6.2832