如何编写一个GAUSS积分的程序
我只是想知道如何编写的程序可以对任意的带参数的函数,任何的上下限(当然是可以积的函数)都可以进行积分我写的是
function result=jifen(a,b)
result=0;
xj=[-0.90617985 -0.53846931 0 0.53846931 0.90617985];
Aj=;
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那样,请问如何修改?
多谢
页:
[1]