In this work, we present a numerical scheme based on the operational matrix of fractional Caputo-Fabrizio (CF) integration for handling fractional Bloch equation (FBE) in nuclear magnetic resonance (NMR). The understanding of Bloch equation provides us a fundamental framework for describing magnetic resonance phenomena, facilitating breakthrough in diverse fields such as medical diagnostics, quantum computing and materials characterization. The non-integer order derivative and integration are presented in the Caputo-Fabrizio sense. To construct the operational matrix, Jacobi polynomial is used as a basis. The fractional Bloch equation is transformed into a set of algebraic equations by using the operational matrix. In order to examine the fractional order problem, we obtain an approximate solution for FBE and present the numerical results in graphical and tabular forms.