function analyse

[S, lambdai] = CIE_daylight([5000 6000, 7000, 10000, 15000], 300:5:830);

load dxc_930p.dat
figure, plot(dxc_930p(:,1), dxc_930p(:,2:end))

load CV-M7_tristimulus
figure, plot(300:5:830, CV_M7)