function y = least_square_bpp(x0,spec,unc,wl);

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

