| 查看: 1440 | 回复: 3 | |||
ghost118木虫 (小有名气)
|
[求助]
matlab龙贝格求积 已有1人参与
|
|
小弟按照书上的算法程序输入了matlab，被积函数是ln(1+x)/x,积分区间是0到1，下面是算法程序： a=0;b=1;eps=0.5*10^-6; M=1;h=b-a;err=1;J=0;R=zeros(4,4); R(1,1)=h*(1+f(b))/2; while(err>eps) J=J+1; h=h/2; x=a+h:2*h:b-h; R(J+1,1)=R(J,1)/2+h*sum(f(x)); for K=1:min(3,J) R(J+1,K+1)=R(J+1,K)+(R(J+1,K)-R(J,K))/(4^K-1); end if(J>3) err=abs(R(J+1,4)-R(J,4)); end end quad=R(J,4); 运行出来的结果和实际相差很大，而且用这个程序做书上的例题，可以做出来，是不是被积函数是反常积分的缘故呢，希望能有大神从中指点，不胜感激 [ 发自手机版 http://muchong.com/3g ] |
回复此楼» 猜你喜欢
今年审到国自然15份，谈谈感受
已经有18人回复
-
26/27博士推荐
已经有5人回复
-
评审有感
已经有19人回复
-
献血感触
已经有11人回复
-
上海大学实验技术岗位非升即走
已经有12人回复
-
如何实现卤原子转化
已经有7人回复
-
云南大学材料与能源学院解琳课题组钙钛矿博士招生
已经有5人回复
-
2026山东省优青
已经有4人回复
-
反应很差，大量原料没有反应
已经有3人回复
-
Sci. Bull. 悲剧经验
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 亲们，求助一下，这个非线性方程组怎么解？
已经有11人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68143.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.6小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2014-11-12 18:49:39
ghost118
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1773.7
- 散金: 148
- 红花: 1
- 帖子: 286
- 在线: 136.3小时
- 虫号: 1968273
- 注册: 2012-09-02
- 性别: GG
- 专业: 钢铁冶金
3楼2014-11-13 00:16:51
ghost118
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1773.7
- 散金: 148
- 红花: 1
- 帖子: 286
- 在线: 136.3小时
- 虫号: 1968273
- 注册: 2012-09-02
- 性别: GG
- 专业: 钢铁冶金
4楼2014-11-13 00:17:33
5