This function implements the Inertia Moment (IM) method  [1]. 
  This method can be used with different normalizations over the co-occurrence 
  matrix. Thus, they can be:

  $Y\approx E[(i-j)^2]$ 

  TYPE 1: The function uses the normalized co-occurrence matrix (COM) proposed by 
          CARDOSO, R.R. et al. [2]: $Y$.
  TYPE 2: The function uses the normalized co-occurrence matrix (COM) proposed by 
          ARIZAGA, R. et al. [1]: $Y2$.

  References:
  [1]  ARIZAGA, R. et al. Speckle time evolution characterization by the 
       co-occurrence matrix analysis. Optics and Laser Technology, Amsterdam, 
       v. 31, n. 2, p. 163-169, 1999.
  [2]  BRAGA R.A. CARDOSO, R.R. Enhancement of the robustness on dynamic speckle 
       laser numerical analysis. Optics and Lasers in Engineering, 
       63(Complete):19-24, 2014.

  After starting the main routine just type the following command at the
  prompt:
           [Y] = inertiamoment(COM);
        [Y Y2] = inertiamoment(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 in the same position in the output.
       If TYPE is equal to 2, the function also returns the inertia moment
       with ARIZAGA [2] co-occurrence normalization.

  Output:
  Y     is the value of inertia moment [1] with CARDOSO normalization [2].

  Ytype if TYPE is equal to 2, the function also returns the inertia moment. 
        ARIZAGA [2] co-occurrence normalization.


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

Package: bsltl