Automatic calculation of two-loop ELWK corrections to the muon ($g$-2)

Two-loop electroweak corrections to the muon anomalous magnetic moment are automatically calculated by using GRACE-FORM system, as a trial to extend our system for two-loop calculation. We adopt the non-linear gauge (NLG) to check the reliability of our calculation. In total 1780 two-loop diagrams consisting of 14 different topological types and 70 one-loop diagrams composed of counter terms are calculated. We check UV- and IR-divergences cancellation and the independence of the results from NLG parameters. As for the numerical calculation, we adopt trapezoidal rule with Double Exponential method (DE). Linear extrapolation method (LE) is introduced to regularize UV- and IR- divergence and to get finite values.