10 REM *********** ΠΡΟΓΡΑΜΜΑ
Π114 ************
20 REM * *
30 REM * x²-Kατανομή *
40 REM
*
*
50 REM
******************************************
60 CLS
70
EPS=.0000005
80 DIM
X(24),W(24)
90 INPUT "Βαθμοί
ελευθερίας ν της x²=";NN
100 PRINT :INPUT "x²";X2:PRINT
110
AL=.5*NN
120
B=.5*X2
130
PI=2*SQR(ATN(1))
140 GOSUB 350
150 IF
INT(NN/2)*2=NN THEN GOTO 190
160 NC=(NN-1)/2
170 GOSUB 240
180 GOTO
220
190 FOR
I=2 TO NN/2-1
200 SUM=SUM/I
210 NEXT
I
220
PRINT "Κατανομή
x²= ";SUM
230 END
240 IF
NC=0 THEN GOTO 330
250 IF
NC=1 THEN GOTO 310
260 FOR
I=1 TO NC
270 SUM=2*SUM/(2*I-1)
280 NEXT
I
290
SUM=SUM/PI
300 RETURN
310
SUM=2*SUM/PI
320
RETURN
330
SUM=SUM/PI
340 RETURN
350 REM Το πρόγραμμα αυτό
υπολογίζει την
360 REM ελλειπτική Γάμμα συνάρτηση γ(α,x)
370 REM με τη σύνθετη μέθοδο
του Gauss
380 DEF
FNF(AL,T)=T^(AL-1)*EXP(-T)
390 REM Διάβασε δεδομένα
400 FOR
I=1 TO 24
410 READ
X(I),W(I)
420 NEXT
I
430
Q=.00001
440
K=0:A=0
450
TEMP=0
460
N=2^K
470
H=(B-A)/N
480
AH=A-.5*H
490
SUM=0
500 FOR
J=1 TO N
510
XJ=J*H
520 FOR
I=1 TO 24
530
R=.5*H*X(I)
540
Y1=AH+R+XJ
550 Y2=AH-R+XJ
560
SUM=SUM+W(I)*(FNF(AL,Y1)+FNF(AL,Y2))
570 NEXT
I
580 NEXT
J
590
SUM=.5*H*SUM
600 IF
ABS(SUM-TEMP)<EPS THEN RETURN
610 K=K+1
620 IF
K=6 THEN RETURN
630
TEMP=SUM
640 GOTO 460
650 PRINT
660 END
670 PRINT
680 PRINT "Δεν έχουμε την
ακρίβεια που θέλεις"
690 PRINT"Αλλαξε ακρίβεια
στην 150 ή αύξησε το κ στην 740"
700 END
710 DATA
.0323801709628694,.0647376968126839
720 DATA
.0970046992094627,.0644661644359501
730 DATA
.1612223560688917,.0639242385846482
740 DATA
.2247637903946891,.0631141922862540
750 DATA
.2873624873554556,.0620394231598927
760 DATA
.3487558862921607,.0607044391658939
770 DATA
.4086864819907167,.0591148396983956
780 DATA
.4669029047509584,.0572772921004032
790 DATA
.5231609747222330,.0551995036999842
800 DATA
.5772247260839727,.0528901894851937
810 DATA
.6288673967765136,.0503590355538545
820 DATA
.6778723796326639,.0476166584924905
830 DATA
.7240341309238147,.0446745608566943
840 DATA
.7671590325157403,.0415450829434647
850 DATA
.8070662040294426,.0382413510658307
860 DATA
.8435882616243935,.0347772225647704
870 DATA
.8765720202742479,.0311672278327981
880 DATA
.9058791367155697,.0274265097083569
890 DATA
.9313866907065543,.0235707608393244
900 DATA
.9529877031604309,.0196161604573555
910 DATA
.9705915925462472,.0155793157229438
920 DATA
.9841245837228269,.0114772345792346
930 DATA
.9935301722663508,.0073275539012763
940 DATA
.9987710072524261,.0031533460523058