function CM(ALPHA,EL)
      A=[0.2050    0.0810   -0.0460   -0.1740   -0.2590
         0.1680    0.0770   -0.0200   -0.1450   -0.2020
         0.1860    0.1070   -0.0090   -0.1210   -0.1840
         0.1960    0.1100   -0.0050   -0.1270   -0.1930
         0.2130    0.1100   -0.0060   -0.1290   -0.1990
         0.2510    0.1410    0.0100   -0.1020   -0.1500
         0.2450    0.1270    0.0060   -0.0970   -0.1600
         0.2380    0.1190   -0.0010   -0.1130   -0.1670
         0.2520    0.1330    0.0140   -0.0870   -0.1040
         0.2310    0.1080         0   -0.0840   -0.0760
         0.1980    0.0810   -0.0130   -0.0690   -0.0410
         0.1920    0.0930    0.0320   -0.0060   -0.0050]
      #
      S= 0.2 * ALPHA
      K= floor(S)+3
      if(K <= 1)
            K= 2
      end
      if(K >= 12)
            K=  11
      end
      DA= S - K +3
      L = K + fix( 1.1*sign(DA) )
      #
      S= EL/12.0
      M= floor(S)+3
      if(M <= 1)
            M= 2
      end
      if(M >=  5)
            M=  4
      end
      DE= S - M +3
      N= M + fix( 1.1*sign(DE) )
      #
      K=convert(Int64, K)
      L=convert(Int64, L)
      M=convert(Int64, M)
      N=convert(Int64, N)
      T= A[K,M]
      U= A[K,N]
      V= T + abs(DA) * (A[L,M] - T)
      W= U + abs(DA) * (A[L,N] - U)
      return cmv = V + (W-V)  * abs(DE)
end
