Single particle calculations for a Woods-Saxon potential with triaxial deformations, and large Cartesian oscillator basis

We present a computer program which solves the Schrodinger equation of the stationary states for an average nuclear potential of Woods-Saxon type. In this work, we take specifically into account triaxial (i.e. ellipsoidal) nuclear surfaces. The deformation is specified by the usual Bohr parameters. The calculations are carried out in two stages. In the first, one calculates the representative matrix of the Hamiltonian in the cartesian oscillator basis. In the second stage one diagonalizes this matrix with the help of subroutines of the EISPACK library. If it is wished, one can calculate all eigenvalues, or only the part of the eigenvalues that are contained in a fixed interval defined in advance. In this latter case the eigenvectors are given conjointly. The program is very rapid, and the run-time is mainly used for the diagonalization. Thus, it is possible to use a significant number of the basis states in order to insure a best convergence of the results