function F = lmaxpdf(x,l,w,gbar,Etaext,Rext)
F1 = repmat(expcdf(l/w*((Etaext(3:end)-Etaext(2:end-1))./(Rext(3:end)-Rext(2:end-1))),gbar)',1,length(x));
F2 = repmat(expcdf(l/w*((Etaext(2:end-1)-Etaext(1:end-2))./(Rext(2:end-1)-Rext(1:end-2))),gbar)',1,length(x));
F3 = (repmat(w*Rext(2:end-1)',1,length(x))>repmat(x,length(Rext(2:end-1)),1)).*...
  (exppdf(repmat(l*Etaext(2:end-1)',1,length(x))./(repmat(w*Rext(2:end-1)',1,length(x))...
  -repmat(x,length(Rext(2:end-1)),1)),gbar)).*...
  (repmat(l*Etaext(2:end-1)',1,length(x))./((repmat(w*Rext(2:end-1)',1,length(x))...
  -repmat(x,length(Rext(2:end-1)),1)).^2));
F3(isnan(F3)) = 0;
F3(repmat(l*Etaext(2:end-1)',1,length(x))./max(repmat(w*Rext(2:end-1)',1,length(x))...
  -repmat(x,length(Rext(2:end-1)),1),0)<repmat(l/w*((Etaext(2:end-1)-Etaext(1:end-2))./(Rext(2:end-1)-Rext(1:end-2)))',1,length(x)))=0;
F3(repmat(l*Etaext(2:end-1)',1,length(x))./max(repmat(w*Rext(2:end-1)',1,length(x))...
  -repmat(x,length(Rext(2:end-1)),1),0)>repmat(expcdf(l/w*((Etaext(3:end)-Etaext(2:end-1))./(Rext(3:end)-Rext(2:end-1))),gbar)',1,length(x)))=0;
% F3(w*Rext(2:end-1)-x<0)=1;
F = sum(F3.*(F1-F2));
end