Calcula a integral definida de f(x)dx usando a regra composta de Simpson sobre os nodos

function s=simpson_comp(nodos,func)

    // s=simpson_composto(nodos,f)

    // Calcula a integral definida de f(x)dx usando a regra composta de Simpson,

    // 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

    n=max(size(nodos));

    s = 0;

    for i=1:n-1

        s = s+simpson(nodos(i),nodos(i+1),func);

    end

endfunction