This study investigates the effects of viscous dissipation and a heat source or sink on the magneto-hydrodynamic laminar boundary layer flow of a Jeffrey fluid past a vertical plate. The governing boundary layer non-linear partial differential equations are reduced to non-linear ordinary differential equations using suitable similarity transformations. The resulting system of dimensionless differential equations is then solved numerically using the bivariate spectral quasi-linearisation method. The effects of some physical parameters that include the Schmidt number, Eckert number, radiation parameter, magnetic field parameter, heat generation parameter, and the ratio of relaxation to retardation times on the velocity, temperature, and concentration profiles are presented graphically. Additionally, the influence of some physical parameters on the skin friction coefficient, local Nusselt number, and the local Sherwood number are displayed in tabular form.