In this paper we consider the confluent Heun equation, which is a linear differential equation of second order with three singular points -two of them are regular and the third one is irregular of rank 1. The purpose of the work is to propose a procedure for numerical evaluation of the equation's solutions (confluent Heun functions). A scheme based on power series, asymptotic expansions and analytic continuation is described. Results of numerical tests are given. arXiv:1804.01007v1 [math.NA] 2 Apr 2018 2 Statement and basic notationsWe use the following form of the confluent Heun equation:This second order linear differential equation has regular singularities at z = 0 and 1, and an irregular singularity of rank 1 at z = ∞ (see e.g. [20]). The parameter q ∈ C is usually referred to as an accessory or auxiliary parameter and γ, δ, ε, α (also belonging to C) are exponent-related parameters.It is important to note that in this paper the parameters ε and α are assumed to be independent. Below, we will use notation c H (q, α, γ, δ, ε; z) or c H (z) for brevity. There are 6 local solutions of equation (1). The Frobenius method can be used to derive local power-series solutions to (1) near z = 0 and z = 1 (two per a singular point), while two solutions at z = ∞ can be obtained in the form of asymptotic series. In § 3 we will present the local solutions near the point z = 0. One of them is analytic in a vicinity of zero and if γ is not a nonpositive integer, we normalize this solution to unity at z = 0 and call it the local confluent Heun function. It is denoted by c Hl (q, α, γ, δ, ε; z). For the second Frobenius local solution, we will use the notation c Hs(q, α, γ, δ, ε; z). When γ is a nonpositive integer, one solution of (1) is analytic in a vicinity of z = 0 but it is equal to zero at z = 0, whereas the second solution can be normalized to unity at zero but generally it is not analytic. Following [16], the normalized solution will be denoted by c Hl (z) and another one by c Hs(z).