| 查看: 2277 | 回复: 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 ] |
» 猜你喜欢
梦想
已经有7人回复
学科评审组评审是指会评吗?
已经有4人回复
面上合作单位盖章
已经有8人回复
面上函评意见出来了,像什么等级?
已经有15人回复
为什么资助数各大高校都创新高,自己申请怎么就这么难
已经有14人回复
小白求助 投论文要求的highlights应该如何写
已经有4人回复
麻烦专家们看看评委们的意见(F口面上)
已经有12人回复
怎么看青基中了没有啊
已经有6人回复
国社科又开始会评了,不知道这次命运如何
已经有11人回复
有多少人是今天查系统知道结果的?
已经有19人回复
» 本主题相关价值贴推荐,对您同样有帮助:
FORTRAN新手 求助主程序循环问题
已经有10人回复
fortran中到底该如何比较两个实数是否相等
已经有8人回复
怎能让fortran调用c的fopen和fileio函数
已经有3人回复
Fortran中子程序与函数的区别
已经有7人回复
如何改变fortran的浮点精度
已经有28人回复
【求助】如何写FORTRAN程序实现求平均最近邻距离
已经有34人回复
» 抢金币啦!回帖就可以得到:
【征稿中!】2026年机械、制造技术与材料国际会议 (MEMTM 2026)
+1/905
【IEEE列表会议征稿 】2026年第六届控制理论与应用国际会议(ICoCTA 2026)
+5/870
云南大学丁军桥团队2027年博士招生
+1/811
优质男征女友
+1/469
昆明理工大学智能车辆技术团队招聘博士后
+3/222
【以色列理工学院全额奖学金硕博士招生】王燕课题组 · 柔性生物界面实验室
+2/118
锡基钙钛矿|超低水氧惰性环境,解锁无铅光伏技术攻关路径
+1/93
天平合规性问题
+1/92
南方科技大学刘科院士团队诚聘矿物浮选分离方向研究助理
+1/83
大湾区大学刘天辉课题组招收2027级博士生(与中山大学物理学院联合培养)
+1/42
大湾区大学刘天辉课题组招收2027级博士生(与中山大学物理学院联合培养)
+1/42
山东大学-山东字节信息智能技术实验室2026年招聘博士后
+2/22
评职称需要,卖发明专利了:D有技术储备需求的设备厂商可联系
+1/19
浙江大学物理高等研究院袁怀洋课题组诚聘博士和博士后
+1/15
中国药科大学陈达威课题组招收2027年博士研究生 (药理学)
+1/11
南方科技大学环境科学与工程学院陈洪副教授课题组2026年科研助理招聘启事
+1/5
义翘神州推出全长p-Tau217蛋白,助力阿尔茨海默病下一代生物标志物检测开发
+1/5
项目助理招聘启事
+1/2
义翘神州TR-FRET激酶检测试剂盒助力高通量药物筛选
+1/1
义翘神州TR-FRET激酶检测试剂盒助力高通量药物筛选
+1/1
6楼2011-03-26 16:16:39
2楼2011-03-26 10:13:36
3楼2011-03-26 11:01:33
4楼2011-03-26 11:52:59










回复此楼
on 2011-3-26 at 16:18 ]
投票: