function s=trapezio_comp(nodos,func)
// s=trapezio_composto(nodos,f)
// Calcula a integral definida de f(x)dx usando a regra composta do trapézio,
// sobre os nodos "nodos".
// A função "f" deve ser passada como uma "string" contendo a função,
// escrita de acordo com as regras do Scilab; por exemplo, para a
// função f(x)=e^x-2, escreva 'exp(x)-2'.
// (C) Rudnei Dias da Cunha 2008
deff('[v]=f(x)','v='+func);
n=max(size(nodos));
s = 0;
for i=1:n-1
s = s+(nodos(i+1)-nodos(i))*(f(nodos(i))+f(nodos(i+1)));
end
s = s/2;
endfunction
Nenhum comentário:
Postar um comentário