声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1005|回复: 0

[综合讨论] 如何编写一个GAUSS积分的程序

[复制链接]
发表于 2007-9-3 21:17 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我只是想知道如何编写的程序可以对任意的带参数的函数,任何的上下限(当然是可以积的函数)都可以进行积分
我写的是
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那样,请问如何修改?
多谢
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-6-7 16:06 , Processed in 0.046270 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表