| 查看: 1344 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
δ函数或者comb怎么编程
|
||
| 如何在matlab上编写二维梳状函数 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有25人回复
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有30人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有10人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
-
filecode
已经有16人回复
★ ★ ★ ★ ★ ★
华丽的飘过: 金币+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
