simpson

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

function y=f(x)
 y=cos(x)+cos(3*x)
endfunction
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
 x1=x(i)
 x3=x(i+1)
 x2=x1+h/2
 x2=x1+(x3-x1)/2
 x2=(x1+x3)/2
 A1 =1/6; A2 =4/6; A3=1/6
 //A1 =-1/6; A2 =8/9; A3=0.277777777777777777777777777777777
 dS =(A1*f(x1)+A2*f(x2)+A3*f(x3))*h
 S=S+dS
end
endfunction