Abstract. We consider the equations of stationary incompressible magnetohydrodynamics posed in three dimensions, and treat the full coupled system of equations with inhomogeneous boundary conditions. We prove the existence of solutions without any conditions on the data. Also we discuss a finite element discretization and prove the existence of a discrete solution, again without any conditions on the data. Finally, we derive error estimates for the nonlinear case.