Le schéma itératif le plus simple est la méthode du pas de temps
imaginaire. Elle consiste à poser dans l'équation de Schrödinger
dépendant du temps 
, d'où son nom. On voit alors qu'en résolvant
cette équation par la méthode d'Euler sur un intervalle assez grand pour
, en orthonormalisant les fonctions d'onde à chaque pas,
la solution tend vers l'état fondamental du
système [56]. Nous avons appliqué cette méthode dans la
version à symétrie axiale de notre programme : dans des cas où la
déformation de l'agrégat, elle donne des résultats acceptables, mais en
général il faut entre 200 et 300 itérations avec un pas de temps
imaginaire très petit pour aboutir à une solution correcte.