24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1331  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

suthing

木虫 (正式写手)

[求助] fortran中如何求奇点在边界的积分~~

公式所式,如何实现呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishuxian

银虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 suthing 的主题更新
信息提示
请填处理意见