formula de simpson aplicada no matlab

disp('Os extremos inferior a e superior b da integral são') , a=1 , b=2 disp('O mpasso h entre os valores tabelados de x é') , h=0.25 x=a:h:b ; n=round((b-a)/h) disp('Escreva a expressão algébrica de y=f(x)') , y=2*x+exp(x)
%Cálculo da integral usando a fórmula de simpson é o seguinte s=y(1)+y(n+1); %primeiro e ultimo ponto da tabela j=2:2:n; %elementos de posição par na tabela s=s+4*sum(y(j)); j=3:2:n-1; %elementos de posição impar na tabela s=s+2*sum(y(j)); s=h*s/3; %formula de simpson fprintf('Valor aproximado da integral é: %8.5f',s) fprintf('\n') %pular linha

disp('O erro de truncamento é calculado por: ') syms x , df=diff(2*x+exp(x),4); t=a:0.5:b; R=eval(subs(df,t))
M=max(abs(R)) %M é o maior valor da 4º derivada k=n/2 fprintf('O numero de vezes em que simpson é aplicada é:%3.0f\n',k);
Et=k*(h^5/90)*M %Formula de erro de runcamento para simpson
Erro=vpa(Et,8)
fprintf('\n') %pular linha
%O valor exato da integral pode ser calculado com a function int: syms x,
I=int(2*x+exp(x),a,b)
Vexato=vpa(I,8)
Os extremos inferior a e superior b da integral são a =
1

b =
2
O mpasso h entre os valores tabelados de x é h =
0.2500

n =
4

1

Escreva a expressão algébrica de y=f(x) y =
4.7183

5.9903

7.4817

9.2546

11.3891

Valor aproximado da integral é: 7.67087
O erro de truncamento é calculado por:
R =
2.7183

4.4817

7.3891

M =
7.3891

k =
2
O numero de vezes em que simpson é aplicada é:
Et =
1.6035e-04

Erro =
0.00016035278

I = exp(2) - exp(1) + 3

Vexato =
7.6707743

Published with MATLAB® R2014a

2

2