Computer subroutines in C++ for computing Mathieu functions of integer orders are described. The routines can handle a large range of the order n and the parameter ft. Sample test results and graphs are given.