function a = fitpoly(x, y, n)
	a = randermonde(x', n)\y';