如何编写一个GAUSS积分的程序

yanice

我只是想知道如何编写的程序可以对任意的带参数的函数，任何的上下限（当然是可以积的函数）都可以进行积分
我写的是
function result=jifen(a,b)
result=0;
xj=[-0.90617985 -0.53846931 0 0.53846931 0.90617985];
Aj=[0.23692689 0.47862867 0.56888889 0.47862867 0.23692689];
xxj=0.5*(b-a)*(xj+1)+a;
for j=1:5
    result=result+f(xxj(j))*Aj(j);
end
result=0.5*(b-a)*result;
我的只可以计算f(x)，并且只能是在文件夹里的f.m文件的函数才可以，如果要积分其他的函数，还要修改，现在我想把被积函数能够弄成任意的函数都可以，就像quad那样，请问如何修改？
多谢