To understand the nickel iron electrodeposition process, we have developed an one dimensional numerical method. This model addresses dissociation, diffusion, electromigration, convection and deposition of multiple ion species. We consider a system of steady-state transport equations of five dissolved species with simultaneous homogeneous reactions coupled with electroneutrality. To compute the numerical solution of the non linear system of reaction diffusion differential equation obtained we introduce an iterative domain decomposition method. This method manage the non-linear reaction and the anisotropiy of the solution, The convergence of the domain decomposition method is proved. We present some numerical examples.