The matrix form of magnetohydrodynamics (MHD) is presented. The approach follows the recent development of matrix methods in optics. Maxwell's equations, Ohm's law, the equation of continuity, the equation of motion with the JxB body force and the energy equation were formulated in matrix form without using vector or tensor analysis. These equations were applied to describe the propagation of plane MHD waves in conducting fluids. Two transverse modes were obtained due to z, coupling between viscous and magnetic diffution and Alfven wave. The ordinary acoustic waves was found to split into fast, slow and intermediate magnetoacoustic waves.