function f_mel = Hz2mel( f_Hz )
  f_mel =  2595 * log10( 1 + f_Hz / 700 );
  