| 查看: 551 | 回复: 3 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】大家帮我看看这个matlab程序啊？【已解决】
|
|||
|
RT~ 小弟刚学数值分析，自己编的matlab程序，复合-3点高斯-勒让德求积公式，代码见下面！为什么计算pi总是不对啊？pi=4/(1+x^2)，0~1积分。是这样的，这个程序要用很大很大的n才能算稍微准确一点，但是用其他自己编的辛普森公式之类的，只用几段就可以算出来了，应该高斯公式是精度最高的啊？哪位达人抽点时间帮我看一下代码程序？有什么问题啊？？谢谢啦！金币相送！！ function I=fguass(fh,a,b,n) %fh:输入一个函数 %a:积分下限,b:积分上限,n:积分区间区间的分段数 h=(b-a)/n; x=a:h:b; x1=0.5*h-sqrt(15)/10*h+x; x2=0.5*h+x; x3=0.5*h+sqrt(15)/10*h+x; fx1=5/9*fh(x1); fx2=8/9*fh(x2); fx3=5/9*fh(x3); fx=h/2*(fx1+fx2+fx3); I=sum(fx); [ Last edited by nono2009 on 2009-11-8 at 16:33 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有14人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有11人回复
-
论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有13人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有16人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有26人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
4楼2009-11-08 16:05:30
2楼2009-11-08 02:00:05
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
3楼2009-11-08 09:31:16
50