k = 18;
tmp = twones;
for deg = 1:k
	tmp = [tmp; train1.^deg];
end
Y = tmp;
tmp = twones;
for deg = 1:k
	tmp = [tmp; train2.^deg];
end
preY = [Y tmp]';
Y = [Y -1*tmp]';

Yi = pinv(Y);

tmp =Yi*b;

for r = 1:100
	ar = r/10;
	x = [1];
	for deg = 1:k
		x = [x ar.^deg];
	end
	x
	if((tmp'*x') > 0)
		disc(r) = 1;
	else
		disc(r) = -1;
	end
end

return;
