n:= READ; m:= READ; s:= 0;
	begin	real array  X[1:n, -2:m];
		for i:= 1 step 1 until n do
		for j:=-1 step 1 until m do X[i,j]:= READ;
		for i:= 1 step 1 until n do
		begin	a:= b:= 0;
			for j:= 1 step 1 until m do a:= a+X[i,j];
			X[i,-1]:= a:= a/m;
			for j:= 1 step 1 until m do b:= b+(X[i,j]-a)*(X[i,j]-a);
			s:= s+b; X[i,-2]:= b/(m-1);
		end;
		s:= s/(n*m-n);

		PRINTTEXT(" result program apv 300969  5238
berekening constanten voor vergelijking van twee steekproeven

aantal steekproeven                     "); ABSFIXT(4,0,n); PRINTTEXT("
omvang steekproeven                     "); ABSFIXT(4,0,m);PRINTTEXT ("
geschatte spreiding                     "); ABSFIXT(4,4,sqrt(s)); b:= sqrt(s);
		PRINTTEXT("

tweezijdig toetsen:

onbetrouwbaarheid        constante        constante * spreiding"); NLCR;
		for p:= .1, .05, .01 do
		begin	a:= proct(n*m-n, p/2) * sqrt(2/m);
			SPACE(6); ABSFIXT(2,2,p);
			SPACE(11); ABSFIXT(4,4,a);
			SPACE(11); ABSFIXT(4,4,a*b);NLCR
		end;
		PRINTTEXT("

eenzijdig toetsen:

onbetrouwbaarheid        constante        constante * spreiding"); NLCR;
		for p:= .1, .05, .01 do
		begin	a:= proct(n*m-n, p) * sqrt(2/m);
			SPACE(6); ABSFIXT(2,2,p);
			SPACE(11); ABSFIXT(4,4,a);
			SPACE(11); ABSFIXT(4,4,a*b);NLCR
		end;
		NEWPAGE; PRINTTEXT("
  code      gemiddelde      spreiding    kans");
		for i:= 1 step 1 until n do
		begin	NLCR; ABSFIXT(8,0,X[i,0]);
			ABSFIXT(4,3,X[i,-1]); ABSFIXT(4,3,X[i,-2]);
			ABSFIXT(1,4,Chisq(m-1, (m-1)*X[i,-2]/s))
		end
	end
end
apv