24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jyd3611

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.8
帖子: 33
在线: 10小时
虫号: 2059571
注册: 2012-10-14
性别: GG
专业: 光学

[求助] δ函数或者comb怎么编程

如何在matlab上编写二维梳状函数

回复此楼

» 猜你喜欢

存款400万可以在学校里躺平吗已经有22人回复
最失望的一年已经有4人回复
拟解决的关键科学问题还要不要写已经有7人回复
国自然申请面上模板最新2026版出了吗？已经有19人回复
请教限项目规定已经有3人回复
基金委咋了？2026年的指南还没有出来？已经有10人回复
基金申报已经有6人回复
推荐一本书已经有13人回复
疑惑？已经有5人回复
溴的反应液脱色已经有7人回复

1楼 2013-04-16 20:46:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Hades_h

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 962.8
红花: 1
帖子: 235
在线: 68.5小时
虫号: 2212832
注册: 2012-12-29
性别: GG
专业: 光学

★ ★ ★ ★ ★ ★
华丽的飘过: 金币+6, 3q 2013-04-26 03:52:51

引用回帖:

6楼: Originally posted by jyd3611 at 2013-04-18 23:15:25
能把你变过的代码贴上吗，我参考下，谢谢了！！！...

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

!子函数结束

!----------------------------------------------------------------