Determination of the maximum value of u for D-Flandrin distribution
Y = umaxdfla(u)
real vector
value of the function (H(u)+u/2)/(H(u)-u/2)-FMAX/FMIN.
umaxdfla is the function Y(u)=(H(u)+u/2)/(H(u)-u/2)-fmax/fmin. Doing UMAX = fzero('umaxdfla',0); gives the maximum value for U in the computation of the D-Flandrin distribution. For this distribution, H(u) = 1+(u/4)^2.