In this work, linear stability of an electrically conductive fluid experiencing Poiseuille flow for minimum Reynolds value under a normal magnetic field is analysed using the Chebyshev collocation method. The neutral curves of linear instability are derived by utilising Qualitat and Zuverlassigkeit (QZ) method. Instability of the magnetohydrodynamics for plane Poiseuille flow is introduced by solving the generalised Orr–Sommerfeld equation to determine the growth rates, wave number and spatial shapes of the eigenmodes. To solve linear problems, we use numerical methods which help us at each time step of the simulation, uncoupled by physical processes, which can improve the computational performance. This article provides the stability and error analysis, presents a concise study of the Poiseuille flow, and produces computational tests to support the given theory.