Split the code into functions because Q3 using the same code as
Q1 and Q2.

Created the functions :
*FindBasisFunctions which gets as an input an spd matrix and returns an
orderd list of basis functions best representing the SpdMatrix
*GetRGBValues which return a normalized rgb values for the input spd matrix
*PrintMeanSquareError which:
 	1)Print the mean squared error between the true illuminants
	  and their approximation (sum MSE over all illuminants)
	  using 1,2,3,4 basis functions
	2)Returns the approximations matrix which is a container
	  for all aproximations
Added Q3 which find basis for macbeth surfaces
Adde the helper function display_colors.m from hagit's color tutorials page