Brief description of a FORTRAN 77 program for calculating energy values, refection and transmission matrices, and corresponding wave functions in a coupled-channel approximation of the adiabatic approach is presented. In this approach, a multidimensional Schro¨dinger equation is reduced to a system of the coupled second-order ordinary differential equations on a finite interval with the homogeneous boundary conditions of the third type at the leftand right-boundary points for continuous spectrum problem, or a set of first, second and third type boundary conditions for discrete spectrum problem. The resulting system of these equations containing the potential matrix elements and first-derivative coupling terms is solved using high-order accuracy approximations of the finite element method.