Errors of (a) q2(a1) and (b) q1(a3) computed by discrete variation method (solid line) and Runge–Kutta method (dotted line).