function res = phi(x)
	res = x^2 * (x-2)^2;