| 查看: 2258 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【交流】合流超比函数的 fortran程序（自己搞定了，不过也希望你来交流）
|
|||
|
哪位虫子手里有合流超比函数（Hypergeometric 1F1）的fortran程序，能不能分享一下，我自己写了一个，可是和mathematic算的结果不一样。 1F1(a,b,c)=1+a*c/(b*1!)+a(a+1)*c^2/(b(b+1)*2!) +a(a+1)(a+2)*c^3/(b(b+1)(b+2)*3!)+...... 先在这里谢谢了。。。。。。。 PS：哈哈，刚才写了一个，感觉大体上能对上了，为了以后大家用着方便，我就贴一下我写的程序，希望你来一起交流，呵呵...... fortran77格式的： program fhgfh implicit real*8(a-z) a=0.5 garm=2.5 f=4 z=confluent(a,garm,f) write(*,*) z end c confluent hypergeometric function real*8 function confluent(a,garm,z) implicit real*8(a-z) confluent=1.0 do 11 ii=1,120 confluent=confluent+lianchen(a,ii)*z**ii/ 1 (lianchen(garm,ii)*FACTORIAL(ii)) 11 CONTINUE return END REAL*8 FUNCTION lianchen(a,N) IMPLICIT REAL*8(A-Z) lianchen=a do 22 ii=1,N-1 lianchen=lianchen*(a+ii) 22 CONTINUE return END REAL*8 FUNCTION FACTORIAL(N) IMPLICIT REAL*8(A-Z) FACTORIAL=1.0 DO 101 I=1,N FACTORIAL=FACTORIAL*I 101 CONTINUE return END [ Last edited by janecool on 2011-3-25 at 20:30 ] |
回复此楼» 猜你喜欢
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有20人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- FORTRAN新手 求助主程序循环问题
已经有10人回复
- fortran中到底该如何比较两个实数是否相等
已经有8人回复
- 怎能让fortran调用c的fopen和fileio函数
已经有3人回复
- Fortran中子程序与函数的区别
已经有7人回复
- 如何改变fortran的浮点精度
已经有28人回复
- 【求助】如何写FORTRAN程序实现求平均最近邻距离
已经有34人回复
7楼2013-12-03 00:35:48
on 2011-3-26 at 10:23 ] 2楼2011-03-26 10:13:36
3楼2011-03-26 11:01:33
4楼2011-03-26 11:52:59
