| 查看: 572 | 回复: 5 | ||
| 当前主题已经存档。 | ||
| 本帖产生 1 个 计算强帖 ,点击这里进行查看 | ||
[交流]
【求助】求助一个积分问题!! 已有3人参与
|
||
tianlangxingaa
铁杆木虫 (著名写手)
- 计算强帖: 4
- 应助: 52 (初中生)
- 贵宾: 0.02
- 金币: 7464.2
- 红花: 19
- 帖子: 1283
- 在线: 735.1小时
- 虫号: 438857
- 注册: 2007-09-30
- 专业: 理论和计算化学
2楼2010-04-12 22:09:57
tianlangxingaa
铁杆木虫 (著名写手)
- 计算强帖: 4
- 应助: 52 (初中生)
- 贵宾: 0.02
- 金币: 7464.2
- 红花: 19
- 帖子: 1283
- 在线: 735.1小时
- 虫号: 438857
- 注册: 2007-09-30
- 专业: 理论和计算化学
★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
dongdong3881(金币+2, 计算强帖+1):多谢回复。您会编程厉害 学习 2010-04-13 08:27
小木虫(金币+0.5):给个红包,谢谢回帖交流
dongdong3881(金币+2, 计算强帖+1):多谢回复。您会编程厉害 学习 2010-04-13 08:27
|
n = 200 // 取值区间-n/2到n/2 length = 0.1 // 分成很多小片段,length是小片段的长度 number = n/length // 分成片段个数 do i = 1, number x_lower = i*length - n/2 //每个小片段的x的下限 x_upper = x_lower + length //每个小片段的x的上限 x_middle = (x_lower + x_upper) * 0.5 // x片段的中间值 y1 = x_middle - u1 // 积分表里的一项 y2 = x_middle - u2 // 积分表里的第2项 y3 = y1**2 + alpha*y2**2 // 积分项的分子 y4 = exp(-y3/(sigma**2)) // exp的数值 sum= sum + length * y4 // 积分并且求和 end do [ Last edited by tianlangxingaa on 2010-4-12 at 23:17 ] |
3楼2010-04-12 23:16:04
tianlangxingaa
铁杆木虫 (著名写手)
- 计算强帖: 4
- 应助: 52 (初中生)
- 贵宾: 0.02
- 金币: 7464.2
- 红花: 19
- 帖子: 1283
- 在线: 735.1小时
- 虫号: 438857
- 注册: 2007-09-30
- 专业: 理论和计算化学
4楼2010-04-12 23:18:57
5楼2010-04-13 14:55:51
6楼2010-04-15 19:39:33













回复此楼
