|
|
LET N=6
INPUT PROMPT "GAMMA(X) X=":U
SELECT CASE N
CASE 2
RESTORE 2
CASE 3
RESTORE 3
CASE 4
RESTORE 4
CASE 5
RESTORE 5
CASE 6
RESTORE 6
CASE 7
RESTORE 7
END SELECT
FOR I=1 TO N
READ X,W
LET S=S+F(X,U)*EXP(X)*W !'ガウス・ラゲール法(半無限区間)
NEXT I
PRINT "Γ(";STR$(U);")=";S
LET H=1/1024
FOR T=-5 TO 5 STEP H
LET SS=SS+F(EXP(PI/2*SINH(T)),U)*PI/2*COSH(T)*EXP(PI/2*SINH(T))*H !'二重指数関数法(半無限区間)
NEXT T
PRINT "Γ(";STR$(U);")=";SS
LET A=1
FOR I=1 TO 30000
LET A=A*((1+1/I)^U)/(1+U/I)
NEXT I
PRINT "Γ(";STR$(U);")=";A/U
2 DATA .5857864376269049511983113,8.5355339059327376220042218E-01
DATA 3.4142135623730950488016887,1.4644660940672623779957782E-01
3 DATA .4157745567834790833115339,7.1109300992917301544959019E-01
DATA 2.2942803602790417198220504,2.7851773356924084880144489E-01
DATA 6.2899450829374791968664158,1.0389256501586135748964920E-02
4 DATA .3225476896193923118003615,6.0315410434163360163596602E-01
DATA 1.7457611011583465756868167,3.5741869243779968664149202E-01
DATA 4.5366202969211279832792854,3.8887908515005384272438168E-02
DATA 9.3950709123011331292335364,5.3929470556132745010379057E-04
5 DATA .2635603197181409102030619,5.2175561058280865247586093E-01
DATA 1.4134030591065167922184080,3.9866681108317592745413335E-01
DATA 3.5964257710407220812231866,7.5942449681707595387653311E-02
DATA 7.0858100058588375569221242,3.6117586799220484544612626E-03
DATA 12.6408008442757826594332193,2.3369972385776227891149085E-05
6 DATA .2228466041792606894643548,4.5896467394996359356828488E-01
DATA 1.1889321016726230307431509,4.1700083077212099411337757E-01
DATA 2.9927363260593140776913253,1.1337338207404497573870619E-01
DATA 5.7751435691045105018398304,1.0399197453149074898913303E-02
DATA 9.8374674183825899177155470,2.6101720281493205947924286E-04
DATA 15.9828739806017017825457916,8.9854790642962123882529205E-07
7 DATA .1930436765603624138382479,4.0931895170127390213043288E-01
DATA 1.0266648953391919503451994,4.2183127786171977992928101E-01
DATA 2.5678767449507462069077862,1.4712634865750527839537418E-01
DATA 4.9003530845264845681017144,2.0633514468716939865705615E-02
DATA 8.1821534445628607910818276,1.0740101432807455221319596E-03
DATA 12.7341802917978137580126425,1.5865464348564201268732622E-05
DATA 19.3957278622625403117125821,3.1703154789955805622713222E-08
END
EXTERNAL FUNCTION F(U,X)
LET F=EXP(-U)*U^(X-1)
END FUNCTION
|
|
