Integracao em Matlab

Prévia do material em texto

(Colar cada um dos programas abaixo no editor do matlab e salvar, respectivamente, como integra01.m e integra02.m)
Calcular, pela regra dos trapézios composta, subdividindo-se o intervalo de integração em 6 sub-intervalos o valor de:
format compact
a=3.0; %limite inferior
b=3.6; %limite superior
n=6; %número de sub-intervalos
h=(b-a)/n; %tamanho do sub-intervalo
x=a:h:b; %criando o vetor x (x varia de 'a' a 'b' com passo h
y=1./x; %f(x)=1/x
I=0; %variável cumulativa para calcular a Integral
disp([' i ' ' x ' ' y ' ' c ' ' y*c ']);
for i=0:n
 j=i+1;
 if((i>0)&(i<n)) 
 c=2; %coeficiente c = 2
 else
 c=1; %coeficiente c = 1
 end
 I=I+c*y(j); %acumulando a integral
 disp([i x(j) y(j) c y(j)*c]) %linhas da tabela
end
I=I*h/2 %valor da integral
plot(x,y)
Calcular, pela regra de Simpson composta, subdividindo-se o intervalo de integração em 10 sub-intervalos o valor de:
format compact
a=0.0; %limite inferior
b=1.0; %limite superior
n=10; %número de sub-intervalos
h=(b-a)/n; %tamanho do sub-intervalo
x=a:h:b; %criando o vetor x (x varia de 'a' a 'b' com passo h
y=exp(x); %f(x)=1/x
I=0; %variável cumulativa para calcular a Integral
disp([' i ' ' x ' ' y ' ' c ' ' y*c ']);
for i=0:n
 j=i+1;
 if((i>0)&(i<n))
 if(c==4)
 c=c-2; %coeficiente c = 2
 else
 c=4; %coeficiente c = 4
 end
 else
 c=1; %coeficiente c = 1
 end
 I=I+c*y(j);
 disp([i x(j) y(j) c y(j)*c]) %linhas da tabela
end
I=I*h/3 %valor da integral
plot(x,y);
_1230968382.unknown
_1231049644.unknown