function q=findqopt1(time)
%
% find fastest prob. dist, given timing vector 
%
lam=nlinfit(time,1,'findl1',-10^-6);
q=opt1(time,lam);
