The unsteady laminar magnetohydrodynamic flow over an unsteady stretching sheet in the presence of thermal radiation and variable heat flux is investigated. The governing time-dependent boundary layer equations are transformed into ordinary differential equations containing a radiation parameter, space index parameter, time index parameter, Prandtl number, magnetic parameter, and unsteadiness parameter. These equations are solved numerically by applying the Chebyshev spectral method. The velocity profiles, temperature profiles, the skin friction coefficient, and the dimensionless surface temperature are computed and discussed in detail for various values of the different parameters. It is found that increasing the unsteadiness parameter leads to a fall for both the velocity and temperature distribution. Moreover, the temperature increases with the magnetic parameter, but the reverse is true for the velocity distribution. Likewise, the temperature decreases for increasing both values of space index parameter and time index parameter.