Sharp PC-E500/E500S Technikersoftware

GAMMA.161

Gamma-Funktion

GAMMA.161 (969 Bytes)

10 ' GAMMA.161 (gamma function)
15 CLEAR
20 MERGE "G:COMN.00$"
30 ON ERROR GOTO *ERR
40 GOSUB *CG
50 ERASE P1:DIM P1(8)
60 P1(0)=1:P1(1)=-0.577191652:P1(2)=0.988205891:P1(3)=-0.897056937:P1(4)=0.918206857
65 P1(5)=-0.756704078:P1(6)=0.482199394:P1(7)=-0.193527818:P1(8)=0.035868343
70 GOSUB 900:QX=0:PRINT "x (0<x";CHR$ &E5;"70)= ";:INPUT QX:IF QX>70 OR QX<=0 THEN 70
80 QX$=MID$ (STR$ QX,2,16)
90 GOSUB 900:PRINT "x=";QX$:GOSUB *DPC
100 IF QX=INT QX THEN R2=FACT (QX-1):GOTO 170
110 IF QX<1 THEN QY=QX:GOSUB 400:R2=R2/QY:GOTO 170
120 IF QX<=2 THEN QY=QX-1:GOSUB 400:GOTO 170
140 QY=QX-INT QX:GOSUB 400:FOR QI=1 TO QX-1:R2=R2*(QY+QI):NEXT
170 QXQ=R2:QYQ=-7:GOSUB *RNDN:R2=QXQ
180 GOSUB 900:PRINT CHR$ &83;"(";QX$;")=";R2:GOSUB *REK:GOTO 70
400 R2=P1(8):FOR QI=7 TO 0 STEP -1:R2=R2*QY+P1(QI):NEXT :RETURN
900 CLS :PRINT "        *  GAMMA FUNCTION  *":LOCATE 1,1:RETURN
1000 *ERR IF ERL =70 RESUME 70
1010 GOSUB 900:GOSUB *DPA:GOSUB *REK:RESUME 70

Generated by (c) Casio FX-850P/FX-880P Notepad 2.02

(c) by mani.becker@web.de - http://blog.manib.bplaced.net