function [A, b] = getproblem(n)
	A = gallery('moler', n);
	b = A*ones(n, 1);
