10 REM *********** ΠΡΟΓΡΑΜΜΑ Π116 ************

20 REM * *

30 REM * F - κατανομή *

40 REM * *

50 REM ******************************************

60 CLS

70 EPS=.0000005

80 DIM X(24),W(24)

90 INPUT "υ=";B:PRINT

100 PRINT "Δώσε τους βαθμούς ελευθερίας":PRINT

110 INPUT "ν1 ";N1:PRINT

120 INPUT "ν2 ";N2

130 PI=4*ATN(1)

140 S=1

150 IF N1=1 AND N2=1 THEN GOTO 270

160 IF INT(N1/2)*2=N1 AND INT(N2/2)*2=N2 THEN GOTO 440

170 IF INT(N1/2)*2=N1 OR INT(N2/2)*2=N2 THEN GOTO 290

180 K1=(N1-1)/2

190 K2=(N2-1)/2

200 R=0

210 S=1

220 FOR I=1 TO K1+K2

230 R=R+1

240 S=2*S*I/(2*R-1)

250 IF I=K1 THEN R=0

260 NEXT I

270 S=S/PI

280 GOTO 540

290 IF INT(N1/2)*2=N1 THEN GOTO 340

300 K1=(N1-1)/2

310 K2=N2/2

320 L=K2

330 GOTO 370

340 K1=N1/2

350 K2=(N2-1)/2

360 L=K1

370 S=1

380 FOR J=1 TO L

390 S=S*(2*(K1+K2)+1-2*J)/2

400 IF J=1 THEN GOTO 420

410 S=S/(J-1)

420 NEXT J

430 GOTO 540

440 K1=N1/2

450 K2=N2/2

460 IF K1<K2 THEN L=K1

470 IF K1>K2 THEN L=K2

480 S=1

490 FOR J=1 TO L

500 S=S*(K1+K2-J)

510 IF J=1 THEN GOTO 530

520 S=S/(J-1)

530 NEXT J

540 GOSUB 580

550 FD=N1^(.5*N1)*N2^(.5*N2)*S*SUM

560 PRINT "Η F - κατανομή P(";B;"│";N1;",";N2;")= ";FD

570 END

580 REM Σύνθετη μέθοδο του Gauss

590 DEF FNF(N1,N2,T)=T^(.5*N1-1)*(N1*T+N2)^(-.5*(N1+N2))

600 REM Διάβασε δεδομένα

610 FOR I=1 TO 24

620 READ X(I),W(I)

630 NEXT I

640 A=0

650 K=0

660 TEMP=0

670 N=2^K

680 H=(B-A)/N

690 AH=A-.5*H

700 SUM=0

710 FOR J=1 TO N

720 XJ=J*H

730 FOR I=1 TO 24

740 R=.5*H*X(I)

750 Y1=AH+R+XJ

760 Y2=AH-R+XJ

770 SUM=SUM+W(I)*(FNF(N1,N2,Y1)+FNF(N1,N2,Y2))

780 NEXT I

790 NEXT J

800 SUM=.5*H*SUM

810 IF ABS(SUM-TEMP)<EPS THEN RETURN

820 K=K+1

830 IF K=6 THEN RETURN

840 TEMP=SUM

850 GOTO 670

860 DATA .0323801709628694,.0647376968126839

870 DATA .0970046992094627,.0644661644359501

880 DATA .1612223560688917,.0639242385846482

890 DATA .2247637903946891,.0631141922862540

900 DATA .2873624873554556,.0620394231598927

910 DATA .3487558862921607,.0607044391658939

920 DATA .4086864819907167,.0591148396983956

930 DATA .4669029047509584,.0572772921004032

940 DATA .5231609747222330,.0551995036999842

950 DATA .5772247260839727,.0528901894851937

960 DATA .6288673967765136,.0503590355538545

970 DATA .6778723796326639,.0476166584924905

980 DATA .7240341309238147,.0446745608566943

990 DATA .7671590325157403,.0415450829434647

1000 DATA .8070662040294426,.0382413510658307

1010 DATA .8435882616243935,.0347772225647704

1020 DATA .8765720202742479,.0311672278327981

1030 DATA .9058791367155697,.0274265097083569

1040 DATA .9313866907065543,.0235707608393244

1050 DATA .9529877031604309,.0196161604573555

1060 DATA .9705915925462472,.0155793157229438

1070 DATA .9841245837228269,.0114772345792346

1080 DATA .9935301722663508,.0073275539012763

1090 DATA .9987710072524261,.0031533460523058