% Beispiel für Kap 5.2.2
%% Auswertung größeres Schraubenbeispiel
daten = load('daten/schrauben100.dat');
xi = daten(:,1);
yi = daten(:,2);
n = length(xi);
m = length(yi);
sx = std(xi)
sy = std(yi)
al = 0.01;
c1 = finv(al/2, n-1, m-1)
c2 = finv(1-al/2, n-1, m-1)
Tx = (sx/sy)^2
pW = 1 - fcdf(Tx, n-1, m-1)