【答案】应助回帖

» 猜你喜欢

🌟 比利时新鲁汶大学（UCLouvain）诚邀CSC博士加入Pascal Gehring教授团队已经有0人回复
第一性原理计算方向2026级博士申请 PRB*1，四级484 已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有257人回复
求助VISSIM破解版软件已经有0人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有12人回复

» 本主题相关价值贴推荐，对您同样有帮助:

我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？已经有4人回复
如何改变fortran的浮点精度已经有28人回复
matlab求积分，出现奇点怎么办？已经有8人回复
【求助】如何写FORTRAN程序实现求平均最近邻距离已经有34人回复
【求助】fortran中怎么存储不确定长度的数据【已完结】已经有6人回复
【求助】Fortran程序，如何让exe文件可以单独运行？已经有6人回复

1楼2011-11-16 22:20:44

锐利的碎片

木虫 (正式写手)

star watcher

应助: 136 (高中生)
金币: 3637.1
散金: 252
红花: 22
帖子: 988
在线: 1224.9小时
虫号: 961933
注册: 2010-03-05
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★
cenwanglai(金币+2): ~~ 2011-11-17 13:34:07
suthing(金币+5): 2011-11-17 22:53:58

(%i1) integrate(x/sqrt(1-x), x, 0, 1);
(%o1) 4/3

赞一下(1人)

2楼2011-11-17 10:51:08

suthing

木虫 (正式写手)

应助: 2 (幼儿园)
金币: 6040.6
红花: 1
帖子: 354
在线: 69.7小时
虫号: 891234
注册: 2009-11-02
专业: 原子和分子物理

引用回帖:

2楼: Originally posted by 锐利的碎片 at 2011-11-17 10:51:08:
(%i1) integrate(x/sqrt(1-x), x, 0, 1);
(%o1) 4/3

Fortran中有这样的命令？

3楼2011-11-17 16:20:59

锐利的碎片

木虫 (正式写手)

star watcher

应助: 136 (高中生)
金币: 3637.1
散金: 252
红花: 22
帖子: 988
在线: 1224.9小时
虫号: 961933
注册: 2010-03-05
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★
cenwanglai: 金币+2, 呵呵~~ 2012-04-07 16:24:03

这个又不是没解析解，手动积一下很麻烦吗。

4楼2011-11-17 16:33:17

suthing

木虫 (正式写手)

应助: 2 (幼儿园)
金币: 6040.6
红花: 1
帖子: 354
在线: 69.7小时
虫号: 891234
注册: 2009-11-02
专业: 原子和分子物理

引用回帖:

4楼: Originally posted by 锐利的碎片 at 2011-11-17 16:33:17:
这个又不是没解析解，手动积一下很麻烦吗。

不麻烦，只是后面要用到复杂的公式。这发散问题解决不了，算出来的值不对

5楼2011-11-17 16:41:24

lishuxian

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 47.2
红花: 1
帖子: 131
在线: 35.5小时
虫号: 942968
注册: 2010-01-15
性别: GG
专业: 分析仪器与试剂

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
cenwanglai: 金币+3, 谢谢代码回帖~ 2012-04-07 16:24:32

可以使用标准Monte Carlo积分方法，代码如下，仅供参考：
program main
implicit none
integer :: i, N
real :: t, f

N = 1E5
f = 0.0

call random_seed()

do i = 1, N
call random_number(t)
f = f + t / sqrt(1-t)
enddo

f = f / real(N)

write (*, *) f

end program main

OneStepADay!

6楼2012-04-07 10:58:11