Abstract: Using the series expansion and the iterative multinomial expansion of the characteristic values of the Mathieu functions,Mathieu functions and its addition theorem are evaluated numerically.Calculated results are in largely conformity with the reference 1,and the precision appears an accuracy of better than 10^-5 in the range r≤30 and 0≤ q≤ 50.