function lsf = DavilaGeislerLSFMinutes(distance)
%DavilaGeislerLSFMinutes  Compute Davila/Geisler estimate of 3mm pupil human LSF
%   lsf = DavilaGeislerLSFMinutes(distance)
%
%   Compute the 3mm pupil LSF from Davila and Geisler (1991, Vision
%   Research, 31, 1369-1380), Figure 1. They give the parameters of a sum
%   of Gaussians LSF as a function of passed distance.  Distance passed in
%   minutes of arc.
%
%   Return is normalized to a maximum of 1.
%
%   See also GeislerLSFMinutes, WestPSFMinutes, WestLSFMinutes, LsfToPsf, PsfToLsf

weight1 = 0.409;
sigma1 = 0.417;
sigma2 = 1.42;
lsf = weight1*normpdf(distance,0,sigma1) + ...
       (1-weight1)*normpdf(distance,0,sigma2);
  
lsf = lsf/max(lsf(:));