An algorithm to obtain global solutions of the double confluent Heun   equation
  A procedure is proposed to construct solutions of the double confluent Heun equation with a determinate behaviour at the singular points. The connection factors are expressed as quotients of Wronskians of the involved solutions. Asymptotic expansions are used in the computation of those Wronskians. The feasibility of the method is shown in an example, namely, the Schroedinger equation with a quasi-exactly-solvable potential. 
J. Abad, F. J. Gomez, J. Sesma