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

more »

... 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.