function q=findqopt(time)
%
% find fastest prob. dist, given timing matrix 
%
lam=nlinfit(time,1,'findl',-10^-6);
q=opt(time,lam);
