'  Programma che calcola chi-quadro per 20 lanci di un dado per cui si sono
' ottenuti 3 uno, 4 due, 1 tre, 5 quattro, 3 cinque, 4 sei (confronto con
' distribuz. unif.)
'  Modifcando NC (n. classi), N (n. prove), il DATA e la definiz. delle Pr(i)
' si puo` calcolare il chi-quadro relativo ad altri esperimenti e ad altre
' distribuzioni
'  Al posto di Data/Read si puo` usare un file (se hai avviato il programma
' aprendo il QBasic da MaCoSa, indica il percorso a partire da Macosa; ad es.
' se il file e' XXX e sta nel sottodirectory Sosta di MaCoSa occorre battere
' sosta\XXXX come nome$). Le prime 2 righe del file conterranno NC e N, le
' successive conterranno le coppie  FreqAssolutaOsservata, Probabilita`.
PRINT
10 INPUT "1: usi Data/read o 2: usi un file "; r%
IF r% < 1 OR r% > 2 THEN GOTO 10

IF r% = 2 THEN
 INPUT "nome file"; nome$: OPEN nome$ FOR INPUT AS #1
 INPUT #1, nc: INPUT #1, n
 DIM FrOss(nc), pr(nc)  ' freq. assoluta osservata, probabilita`
 FOR i = 1 TO nc: INPUT #1, FrOss(i), pr(i): NEXT
ELSE
 nc = 6: n = 20
 DIM FrOss(nc), pr(nc)  ' freq. assoluta osservata, probabilita`
 DATA 3, 4, 1, 5, 3, 4
 FOR i = 1 TO nc: READ FrOss(i): NEXT
 FOR i = 1 TO nc: pr(i) = 1 / 6: NEXT
END IF
s = 0: FOR i = 1 TO nc: s = s + FrOss(i): NEXT
IF s <> n THEN PRINT "controlla: FrOss(1)+FrOss(2)+...="; s; "<> N": END
FOR i = 1 TO nc
    FrAtt = n * pr(i)
    chi2 = chi2 + (FrOss(i) - FrAtt) ^ 2 / FrAtt
NEXT
PRINT "chi quadro ="; chi2