La détermination classique des équilibres des plasmas de tokamaks repose sur la résolution de l’équation de Grad-Shafranov. On y suppose que le plasma, vu comme un fluide conducteur, a une vitesse nulle. Expérimentalement, il existe de nombreuses observations d’une rotation spontanée dans les tokamaks notamment dans la direction toroïdale. Cette rotation du plasma s’avère avoir des effets très importants sur l’amélioration du confinement et le passage dans le mode H de confinement amélioré. Pour éclaircir ces observations expérimentales, nous nous proposons de déterminer numériquement les états stationnaires axisymétriques des équations non-linéaires de la magnétohydrodynamique visco-résistive obtenues en réintroduisant le terme convectif. Cette étude doit prendre en compte le forçage dû à la présence d’un champ électrique extérieur dans la direction toroïdale servant à créer le courant toroïdal nécessaire dans un tokamak à la création de la composante poloïdale du champ magnétique.