script Trapezios Repetidos

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

// Programa para calculo da integral aproximada de uma função f(x) através da regra dos trapézios repetida
clear
//iniciando cronômetro
tic();
// Definindo função f(x)
function y=f(x), y=4/(1+x^2), endfunction
//deff('[y]=f(x)','y=4/(1+x^2))')
// Definindo n (numero de sub-intervalos)
n=13;
// Definindo extremos do intervalo de integração [a,b]
a=0;
b=1; // b deve ser maior que a
// Definindo pontos na partição do intervalo de integração [a,b]
h=(b-a)/n;
x0=a;
for i=1:n
 x(i)=x0+i*h;
end
for i = 1:n 
 F(i)=f(x(i));
end
// Calculando I_TR
I_TR=h*(2*sum(F)+f(x0)-f(x(n)))/2;
disp(I_TR, 'I_TR=');
// Marcando tempo no cronômetro
toc();
t=toc();