function y = least_square_cp(x0,spec,wl);

y=sum((spec-x0(1).*(550./wl).^x0(2)).^2);
