Numerical analysis of the modified AVD method [1]. This function implements
  the the numerical method [1] that is a modification
  over absolute difference (Fujii method)  [2]

  $Y \approx E[\frac{|i-j|}{i+j}]$
  $Y2\approx E[\frac{(i-j)^2}{(i+j)^2}]$  

  References:
  [1] Renan O Reis; Roberto Braga; Hector J Rabal.
      Light intensity independence during dynamic laser speckle analysis
  [2] FUJII, H. et al. Evaluation of blood flow by laser speckle image sensing. 
      Applied Optics, New York, v. 26, n. 24, p. 5321-5325, 1987.

  After starting the main routine just type the following command at the
  prompt:
           [Y] = numad(COM);
        [Y Y2] = numad(COM,2);

  Input:
  COM  is a 2D matrix, with 256 lines and 256 columns, that represents the 
       Co-Occurrence Matrix of THSP matrix. The element COM(a,b), in the 
       co-occurrence matrix, represents the quantity of times that, in two 
       successive columns of a THSP matrix, the intensity values jump of 
       a-1 to b-1.
  TYPE [Optional] the function returns an additional result.
       If TYPE is equal to 2, the function also returns the AD second moment. 

  Output:
  Y    is the value of the modified AVD first moment [1].

  Y2   if TYPE is equal to 2, the function also returns the AVD second moment. 


  For help, bug reports and feature suggestions, please visit:
  http://nongnu.org/bsltl/

Package: bsltl