syms x y a1 a2equ1=1/120 - ((2*a2)/81 + (pi*a1)/12 + pi^2*(a2/405 + 1/12) + (pi^3*a1)/120)/pi^4equ2=-(a2/128 + (2*pi*a1)/81 + pi^2*(a2/960 + 2/81) + (pi^3*a1)/405 - pi^4/405)/pi^5[a1,a2]=solve(equ1, equ2, "a1", "a2")a1=vpa(a1)a2=vpa(a2)執行結果：a1 = -50761304626263552301236414542950/140129736313387243/(2814749767106560*pi+8727491006471547)a2 = -675/280259472626774486a1 = -.20616910499519232147350835793861e-1a2 = -.24084823741137452008216073646926e-14