| 查看: 1329 | 回复: 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 ] |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 亲们，求助一下，这个非线性方程组怎么解？
已经有11人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68123.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 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