Recent MRI technology and pulse sequences have become increasingly complex and difficult to understand. However, the Bloch equations, which are the basic principle of MRI, have remained unchanged since it was proposed in 1973. Here we have presented a method to numerically solve the Bloch equations to understand the principle of MRI. This method enables us to reproduce complex pulse sequences on a computer and to understand various phenomena related to pulse sequences.