The numerical study of EMHD laminar flow of Bingham fluid flowing among two Riga plates with thermal radiation has been established. Both the Riga plates, upper and lower, are taken as stationary. Constant heat transfer is considered for both Riga plates. The governing equations for the problem have been transformed into dimensionless non-linear PDEs by using usual transformations. The obtained non-dimensional PDEs have been solved numerically by applying the explicit FDM. MATLAB has been conducted for the numerical simulation. The stability analysis stated that the governing equations will be converged for 0.09, r P 1.00 , 40, 40 mn = and = 4.00 respectively.Finally, a comparison with several published results has been discussed.