| 查看: 1111 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
δ函数或者comb怎么编程
|
||
| 如何在matlab上编写二维梳状函数 |
回复此楼» 猜你喜欢
三无产品还有机会吗
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有7人回复
-
压汞仪和BET测气凝胶孔隙率
已经有4人回复
-
博士申请都是内定的吗？
已经有14人回复
-
谈谈两天一夜的“延安行”
已经有13人回复
-
氨基封端PDMS和HDI反应快速固化
已经有11人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有11人回复
-
论文投稿求助
已经有4人回复
-
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛
已经有3人回复
-
投稿精细化工
已经有6人回复
★ ★ ★ ★ ★ ★
华丽的飘过: 金币+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