| 查看: 1134 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
δ函数或者comb怎么编程
|
||
| 如何在matlab上编写二维梳状函数 |
回复此楼» 猜你喜欢
英文综述是否需要润色及查重
已经有5人回复
-
救命帖
已经有5人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有35人回复
-
最失望的一年
已经有18人回复
-
求推荐英文EI期刊
已经有5人回复
★ ★ ★ ★ ★ ★
华丽的飘过: 金币+6, 3q 2013-04-26 03:52:51
华丽的飘过: 金币+6, 3q 2013-04-26 03:52:51
|
Program main implicit none Real,External :: Gamma Write(*,*) "Gamma=",Gamma(X) !输出Gamma积分结果 End !子函数,计算Gamma积分 REAL Function Gamma(X) implicit none integer i Real X,Y,S,T,U Real :: A(11)=(/0.0000677106,-0.0003442342,0.0015397681,-0.0024467480,0.0109736958,-0.0002109075,0.0742379071,0.0815782188,0.4118402578,0.4227843370,1.0/) IF (X<=0) THEN Write(*,*) "ERR**X<0" Gamma=-1.0 End IF Y=X IF (Y<=1.0) THEN T=1.0/(Y*(Y+1.0)) Y=Y+2.0 ELSE IF (Y<=2.0) THEN T=1.0/Y Y=Y+1.0 ELSE IF (Y<=3.0) THEN T=1.0 ELSE T=1.0 10 IF (Y>3.0) THEN Y=Y-1.0 T=T*Y GOTO 10 END IF ENd IF S=A(1) U=Y-2.0 Do i=1,10 S=S*U+A(i+1) END do S=S*T Gamma=S Return END !子函数结束 !---------------------------------------------------------------- |
jyd3611
7楼2013-04-20 16:10:30
2楼2013-04-17 08:39:43
3楼2013-04-17 15:45:38
4楼2013-04-17 21:20:42