Статистическая обработка результатов эксперимента

CLS
REM Stati4eskay Obrabotka danih
PRINT "Otnoshenie zarplati Y(i) ot 4asov X(i)"
INPUT "Vvedite n="; n
PRINT "Dlya n="; n - 2; "Vvedite tabli4niy koefficient studenta pri 0.80"; "pri f"; n - 1
INPUT "Tt="; Tt
DIM x(n), y(n), tx(n), ty(n)
'4asi
FOR i = 1 TO n
    PRINT "Vvedite x("; i; ")";
    INPUT x(i)
NEXT i
'Zarplata
FOR i = 1 TO n
    PRINT "Vvedite y("; i; ")";
    INPUT y(i)
NEXT i
'matemati4eskoe ogidanie
mx = 0: my = 0
FOR i = 1 TO n
    mx = mx + x(i)
    my = my + y(i)
NEXT i
mx = mx / n
my = my / n
'Dispersia
dx = 0: dy = 0: S = 0
FOR i = 1 TO n
    dx = dx + (x(i) - mx) ^ 2
    dy = dy + (y(i) - my) ^ 2
    S = S + (x(i) - mx) * (y(i) - my)
NEXT i
dx = dx / (n - 1): dy = dy / (n - 1)
'Koefficient variacii
bx = SQR(dx): by = SQR(dy)
'Otklonenie
FOR i = 1 TO n
    tx(i) = (x(i) - mx) / bx
    ty(i) = (y(i) - my) / by
NEXT i
PRINT "N   Y(i)   Zarplata   X(i)   4asi   Tx   Ty"
FOR i = 1 TO n
    PRINT USING "|##   ####.##   ####. ## ##.#### ####.#### "; i; x(i); y(i); tx(i); ty(i)
 
NEXT i
rxy = S / ((n - 1) * bx * by)
tr = (rxy * SQR(n - 1)) / SQR(1 - rxy ^ 2)
IF ABS(rxy) > 0.7 AND ABS(rxy) < 1 THEN
    PRINT "Lineinay 3avisimost' syschestvyet";
    IF rxy > 0 THEN
        PRINT "Lin/zav pologitel'nay";
    ELSE
        PRINT "Lin/3av otricatel'nay";
    END IF
ELSE
    PRINT "Lineynaya 3avisimost' ne syschestvyet";
END IF
PRINT
 
PRINT "mx="; mx;
PRINT "my="; my
PRINT "dx="; dx;
PRINT "dy="; dy
PRINT "bx="; bx;
PRINT "by="; by
PRINT "rxy="; rxy
PRINT "tr="; tr; "Tt="; Tt
IF tr > Tt THEN
    PRINT "Koefficient korrelyacii zna4im"
ELSE
    PRINT "Koefficient korrelyacii ne zna4im"
END IF
END

CLS
REM Stati4eskay Obrabotka danih
PRINT "Otnoshenie Koe`ffitcient aktivnosti Y(i) ot Ionnaia sila rastvora X(i)"
INPUT "Vvedite n="; n
PRINT "Dlya n="; n - 2; "Vvedite tabli4niy koefficient studenta pri 0.80"; "pri f"; n - 1
INPUT "Tt="; Tt
DIM x(n), y(n), tx(n), ty(n)
'Ionnaia sila rastvora 
DATA 0.001, 0.002, 0.005, 0.01, 0.04, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5
 FOR i = 1 TO n
 READ x(i)
 NEXT i
' Koe`ffitcient aktivnosti 
 DATA 0.98, 0.97, 0.95, 0.92, 0.9, 0.84, 0.81, 0.8, 0.8, 0.81, 0.82, 0.84
 FOR i = 1 TO n
 READ y(i)
 NEXT i
 
'matemati4eskoe ogidanie
mx = 0: my = 0
FOR i = 1 TO n
    mx = mx + x(i)
    my = my + y(i)
NEXT i
mx = mx / n
my = my / n
'Dispersia
dx = 0: dy = 0: S = 0
FOR i = 1 TO n
    dx = dx + (x(i) - mx) ^ 2
    dy = dy + (y(i) - my) ^ 2
    S = S + (x(i) - mx) * (y(i) - my)
NEXT i
dx = dx / (n - 1): dy = dy / (n - 1)
'Koefficient variacii
bx = SQR(dx): by = SQR(dy)
'Otklonenie
FOR i = 1 TO n
    tx(i) = (x(i) - mx) / bx
    ty(i) = (y(i) - my) / by
NEXT i
PRINT "N   Y(i)   Koe`ffitcient aktivnosti   X(i)   Ionnaia sila rastvora   Tx   Ty"
FOR i = 1 TO n
    PRINT USING "|##   ####.##   ####. ## ##.#### ####.#### "; i; x(i); y(i); tx(i); ty(i)
 
NEXT i
rxy = S / ((n - 1) * bx * by)
tr = (rxy * SQR(n - 1)) / SQR(1 - rxy ^ 2)
IF ABS(rxy) > 0.7 AND ABS(rxy) < 1 THEN
    PRINT "Lineinay 3avisimost' syschestvyet";
    IF rxy > 0 THEN
        PRINT "Lin/zav pologitel'nay";
    ELSE
        PRINT "Lin/3av otricatel'nay";
    END IF
ELSE
    PRINT "Lineynaya 3avisimost' ne syschestvyet";
END IF
PRINT
 
PRINT "mx="; mx;
PRINT "my="; my
PRINT "dx="; dx;
PRINT "dy="; dy
PRINT "bx="; bx;
PRINT "by="; by
PRINT "rxy="; rxy
PRINT "tr="; tr; "Tt="; Tt
IF tr > Tt THEN
    PRINT "Koefficient korrelyacii zna4im"
ELSE
    PRINT "Koefficient korrelyacii ne zna4im"
END IF
END