Calcula a integral definida de f(x)dx com limites de integração "a" e "b" usando a regra de Simpson.

function s=simpson(a,b,func)

    // s=simpson(a,b,f)

    // Calcula a integral definida de f(x)dx com limites de

    // integração "a" e "b" usando a regra de Simpson.

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

    s = (b-a)/6*(f(a)+4*f((a+b)/2)+f(b));

endfunction