function [mel]=f2mel(f)
  mel=2595 * log10(1+f/700);

  
  