segundaregradeSimpson

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

function I=AreaSsimpson(a,b,x,y,n) // n= subintervalos
 h=(b-a)/n;
 I=y(1)// comeca com o primeiro elemento
 for i=2:n
 if modulo(i,3)==1 
 I=I+2*y(i);
 end
 if modulo(i,3)==2
 I=I+3*y(i);
 end
 if modulo(i,3)==0
 I=I+3*y(i);
 end
 end
 I=(3*h/8)*(I+y(n+1));
endfunction