%%%%%   Este programa resolve uma integral definida pela   %%%%%%
%%%%%       Regra do Trapézio (Integração Numérica)        %%%%%%
function Trapezio
clear all
clc
format long
a = input('Entre com o limite inferior "a" da integral definida: ');
b = input('Entre com o limite superior "b" da integral definida: ');
n = input('Entre com o número "n" de subintervalos: ');
soma=0;
h=(b-a)/n;
for i=1:n-1
 x = a + i*h;
 soma = soma + f(x);
end
Integral = h*(f(a)+ 2*soma + f(b))/2;
disp ('O valor da integral pela Regra do Trapézio foi: ');
disp (Integral)
%
function f=f(t)
f = t;                     
%f = t^2;                     
%f = exp(-t^2/2);
%f = sqrt(t);
%f = sin(t)/t;
%f = sqrt(3*(sin(t))^2 + 5*(cos(t))^2);
%f = nthroot(t^2+1,3); %nthroot(a,n) calcula raiz n-ésima de a
%f = (cos(t))^2/t^2;
%f = exp(sin(t));