The role of differential equations in the process of calculating Feynman integrals is shown. An example of a diagram is given, in the calculation of which the method of differential equations was introduced, the properties of the inverse-mass-expansion coefficients are shown, and modern methods based on differential equations are considered.