The system of 4 differential equations in the external invariant satisfied by the 4 master integrals of the general massive 2-loop sunrise self-mass diagram is solved by the Runge-Kutta method in the complex plane. The method offers a reliable and robust approach to the direct and precise numerical evaluation of Feynman graph integrals.The relevance of the higher order calculations for the comparison with nowadays precision measurements in high energy physics is well known and comprehensively presented by G. Passarino in this conference.Therefore can be of some interest the exploitation of an alternative method (but still in the context of the integration by part identities and master integrals (MI) [1]) to the more common direct integration method for the numerical evaluation of the MI.The method uses directly the differential equations. Starting from the integral representation of the MI, related to a certain Feynman graph, by derivation with respect to one of the internal masses [2] or one of the external invariants [3] and with the repeated use of the integration by part identities, a system of independent first order partial differential equations is obtained in a number equal to the number of the MI (master differential equations).Enlarging the number of loops and legs grows the number of parameters, MI and equations, but does not change or spoil the method.