Análise númérica

[pic]
Engenharia de Produção
Análise Numérica – Ana Carolina
Aluno: Filipi Antunes Barbosa

Integração Numérica - Regra de Simpson 1/3 composta

Program Simpson; var j,n:integer; t,i,a,b,y0,yn,y,jun:real; function f(x:real) :real; begin f:=((x*x*x) + (3*(x*x)) + (3*x)); end; Begin textcolor(white); writeln ('Integração Numérica - Regra de Simpson 1/3 (composta)'); writeln (' '); writeln('Entre com o valor do primeiro intervalo: '); readln(a); writeln('Entre com o valor do segundo intervalo: '); readln(t); writeln('Entre com a quantidade de subintervalos: '); readln(n); begin t:=(t-a)/n; y0:=f(a); yn:=f(t); for j:=1 to n-1 do begin a:=a+t; if (j mod 2) = 0 then begin jun:=(2*f(a)) + jun; end else begin jun:=(4*f(a)) + jun; end; end; i:=(t/3)*(y0+jun+yn); writeln('O valor de integração é:',i); end; end. Tela

Program Simpson; var j,n:integer; t,i,a,b,y0,yn,y,jun:real; function f(x:real) :real; begin f:=((x*x*x*x) + (3*(x*x*x)) + (3*x)); end; Begin textcolor(white); writeln ('Integração Numérica - Regra de Simpson 1/3 (composta)'); writeln (' '); writeln('Entre com o valor do primeiro intervalo: '); readln(a); writeln('Entre com o valor do segundo intervalo: '); readln(t); writeln('Entre com a quantidade de subintervalos: '); readln(n); begin t:=(t-a)/n; y0:=f(a); yn:=f(t); for j:=1 to n-1 do begin a:=a+t; if (j mod 2) = 0 then begin jun:=(2*f(a)) + jun; end else begin jun:=(4*f(a)) + jun; end; end; i:=(t/3)*(y0+jun+yn); writeln('O valor de integração é:',i); end;
end.