The following code calculates the inverse quark propagator of a Wilson-type action with improved Gauge action at one loop in a general covariant gauge. It generates an optimized C++ program to numerically evaluate the amplitude for an arbitrary momentum-point.