Summary: In this paper algorithms for numerical evaluation of the Mittag-Leffler function
and its derivative for all values of the parameters , and all values of the argument C are presented. For different parts of the complex plane different numerical techniques are used. In every case we provide estimates for accuracy of the computation; numerous pictures showing the behaviour of the Mittag-Leffler function for different values of the parameters and on different lines in the complex plane are included. The ideas and techniques employed in the paper can be used for numerical evaluation of other functions of the hypergeometric type. In particular, the same method with some small modifications can be applied for the Wright function which plays a very important role in the theory of partial differential equations of fractional order.