Simpson

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

function S=simpson(a,b,n)
// n numero de intervalos
h=(b-a)/n
x=linspace(a,b,n+1)
S=0
for i=1:n
 A1=1/6; x1=x(i);
 A2=4/6; x2=(x(i)+x(i+1))/2;
 A3=1/6; x3=x(i+1);
 
 A =(A1*f(x1)+A2*f(x2)+A3*f(x3))* h
 S=S+A
end
endfunction
function y=f(x)
 y= cos(x^2 + 3)
endfunction