Magnetic resonance imaging (MRI) is used in physics, chemistry, engineering and medicine to study complex materials. In this paper, numerical solution of fractional Bloch equations in MRI is obtained using fractional variation iteration method (FVIM) and fractional homotopy perturbation transform method (FHPTM). Sufficient conditions for the convergence of FVIM and its error estimate are established. The obtained results are compared with the existing as well as recently developed methods and with the exact solution. The obtained numerical results for different fractional values of time derivative are discussed with the help of figures and tables. Figures are drawn using the Maple package. Test examples are provided to illustrate the accuracy and competency of the proposed schemes.
