General

There are a lot of works[1] on this subject, but, the present situation is far from complete because of the complexity of the whole problem. We have to overcome many technical problems such as calculating tree Feynman diagrams, one-loop diagrams, further two-loop,..., removing infrared and ultraviolet divergences, utilizing gauge invariance, imposing renormalization conditions, and so on. Thus if we want to solve all of these problems in an automatic way, we have to do a number of things. Listing up the main issues required to complete the program for automatic calculation of amplitudes including loop-diagrams, we have the following items:

1. Build the physics model, i.e. construct the Lagrangian and quantize it.

2. Generation of all Feynman diagrams and amplitudes for a given process.

3. Analytic calculation of amplitudes of these diagrams and generation of the expression of the total squared amplitude

4. Feynman parameter integration

5. Phase space integration.

The paper is organized as follows. In Section 2, we give a description about the package FDC1 for Step 1. Section 3 is devoted to the package FDC2 which generate Feynman diagrams for a given process. In section 4, we present FDC3 which is used to perform analytic calculation of each diagram and generate total squared amplitude. Finally, the summary and conclusions are presented in Section 5.