The system of partial differential equations for the C0 function

We present an approach to analyze the scalar integrals of any Feynman diagrams in detail here. This method not only completely recovers some well-known results in the literature, but also produces some new results on the C0 function. The approach can be employed to evaluate the coefficient of arbitrary power of ε in the expansion of a scalar integral, where D=4−2ε denotes the time–space dimension.