| 查看: 1388 | 回复: 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 ] |
回复此楼» 猜你喜欢
297，工科调剂?河南农业大学本科
已经有11人回复
-
22专硕求调剂
已经有11人回复
-
284求调剂
已经有20人回复
-
291求调剂
已经有7人回复
-
296求调剂
已经有5人回复
-
一志愿沪9，生物学326求调剂
已经有4人回复
-
生物学调剂
已经有8人回复
-
22408 312求调剂
已经有19人回复
-
26年电池方向博士申请
已经有4人回复
-
0854调剂
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 亲们，求助一下，这个非线性方程组怎么解？
已经有11人回复
ghost118
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1773.7
- 散金: 148
- 红花: 1
- 帖子: 286
- 在线: 136.3小时
- 虫号: 1968273
- 注册: 2012-09-02
- 性别: GG
- 专业: 钢铁冶金
4楼2014-11-13 00:17:33
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68130.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
5