Several quantum models were widely used
    for simplified interpretation of quantum phenomena.
    Some of them are solved in terms of
    classical special functions. However the others,
    more sophisticated - two Coulomb centers problem,
    Stark effect in hydrogen, isotonic oscillator ets.
    are solved in terms of Heun functions. Typical for
    these is the presence of an additional parameter
    which can play the role of an adiabatic variable.
    The problem posesses the hamiltonian formulation.
    The procedure of antiquantization leads to dynamical
    model which is integrable since it is exposed in
    terms of Painleve equations. The physical understanding
    of these models should be clarified in the future.
    We give results of symbolic and numeric computations
    developed with the help of the Computer Algebra
    package Mathematica.