around \(D=4\) , and the well-known result of \(C_{_{0}}^{(0)}\)  is presented in the
literature[1]}. However, evaluating \(C_{_0}^{(i)}\) 
precisely is also necessary to obtain high order radiative corrections
to the amplitude self-consistently when the virtual interactions
originate from the counter terms[2]}.
Considering Feynman integrals as the generalized hypergeometric functions [3]},
one finds that the \(D-\) module of a Feynman diagram [4]} is isomorphic to
Gel'fand-Kapranov-Zelevinsky (GKZ) \(D-\) module [5]}, [6]}, [7]}, [8]}, [9]}.
Correspondingly Feynman integrals satisfy indeed the systems of
holonomic linear partial differential equations (PDEs) [10]} whose singularities
are determined by the Landau singularities.
