In this paper, the unsteady squeezing flow and heat transfer of MHD nanofluid between the infinite parallel plates with thermal radiation effect is investigated. A similarity transformation is used to convert the governing momentum and energy equations into non-linear ordinary differential equations with the appropriate boundary conditions. These non-linear ordinary differential equations are solved analytically by Duan–Rach Approach (DRA). This method allows us to find a solution without using numerical methods to evaluate the undetermined coefficients. This method modifies the standard Adomian Decomposition Method (ADM) by evaluating the inverse operators at the boundary conditions directly. The effects of various parameters such as the squeeze number, the magnetic parameter, the volume fraction of nanofluid, the Eckert number and the radiation parameter are investigated on the velocity and temperature. Also, the values of skin friction coefficient and the Nusselt number are calculated and presented through figures. The results show that the temperature profile and Nusselt number increase with the increase of radiation parameter. Furthermore, the limiting cases are obtained and are found to be in good agreement with the previously published results.