function y = white_reference(l)

% y = white_reference(lambda)
%
% y = Illumination(lambda)*Spectralon_Reflectance(lambda)*Detector_Response(lambda)
% IN THEORY!
% In practise, use it is some approximate hack.
% lambda is wavelength /[nm].

% Liam Pedersen, NASA ARC, August 1999.


%y = (((l-200)/10).^5).*exp(-((l+3000))/180);
maxx = 650;
k = 1.2;
m = 1.05;
s = 500;
u = maxx - s*((k/m)^(1/m));
%u = 0;
x = (l-u)/s;
y = (x.^k).*exp(-(x.^m)) + 0.5;

my = max(y);
y = y/my;