%x=10*rand(1,100)-5;
%y=10*rand(1,100)-5;
%z=10*rand(1,100)-5;
%u=x+2*y+3*z+4+randn(size(x));
%data=round([x; y; z; u]);
data=[1 4 2 1 3 4 5 -5 4 1 5 0 0 3 -3 0 4 1 3 2 1 -3 2 -4 1 2 2 4 5 3 1 4 1 -5 -4 4 0 3 -3 1 1 -5 1 -1 -5 0 -3 -4 -3 -4 -3 -5 1 -2 0 2 0 0 -1 -4 0 4 4 -2 -3 1 1 -1 -3 4 -4 -4 -4 -3 1 1 -4 4 2 2 -4 4 4 5 4 3 0 -3 -1 -4 -5 4 -2 -2 -2 0 1 -5 3 1;
	4 -2 -1 -4 -3 2 -2 4 -4 5 0 2 5 -2 -1 0 3 3 -4 -3 -1 -4 0 -2 -3 -3 4 2 0 4 -4 2 2 1 -3 1 -2 -4 -3 4 -4 -3 -4 -1 -5 4 -3 -4 -2 0 -4 5 -2 -2 -4 -2 -5 0 3 1 -4 -4 3 4 0 -4 3 -2 -2 2 -5 -5 2 1 0 2 2 3 -2 2 1 -1 -4 3 -2 1 2 -4 -4 0 0 4 3 2 -4 -4 -4 3 4 2;
	-4 2 -4 -4 1 -2 2 2 1 2 -3 2 5 4 -4 -1 -1 2 1 3 -1 -3 -4 3 -3 -1 1 -3 1 0 -3 3 -4 -2 -3 0 -4 -1 -4 -4 3 -2 1 5 -1 2 3 -1 2 -4 4 -3 -2 3 0 3 -1 -2 -5 2 -1 0 1 -4 -2 3 2 -4 -4 -4 -5 -1 2 2 0 -4 1 -4 -4 -4 -4 -3 -3 -2 -2 -3 -2 4 2 1 -3 -3 -4 4 2 1 -2 -3 1 5;
	-1 11 -7 -18 6 5 12 14 3 21 1 14 28 15 -13 0 9 17 4 9 -3 -17 -6 5 -11 -3 17 3 12 12 -14 20 -3 -5 -12 12 -13 -3 -16 1 6 -11 -2 16 -12 17 2 -9 1 -13 6 2 -5 7 -6 12 -8 -2 -6 9 -6 -3 17 -4 -4 5 18 -12 -14 2 -23 -12 11 11 7 -3 7 4 -10 -2 -9 -4 -9 9 -2 0 2 4 3 4 -10 8 -4 23 -1 -5 -6 -4 19 26];
coef=hyperplaneFitViaTls(data)