[Python] 연립방정식 풀기

연립방정식을 파이썬을 이용해서 풀어보자. Pade apporoximation을 통해 f'j를 구하기 위해 아래의 테일러 테이블의 방정식을 풀어야 했다. 테일러 테이블 for Pade approximation 연립방정식은 아래와 같다. 코드를 짜보면, from sympy import symbols, Eq, solve # Pade scheme # 미지수 정의 a0, a1, b0, b1, b2, h = symbols('a0 a1 b0 b1 b2 h') # 연립방정식 정의 eq1 = Eq(b0 + b1 + b2, 0) eq2 = Eq(1 + a0 + a1 + b0*(-1*h) + b2*(h), 0) eq3 = Eq(a0*(-1*h) + a1*h + (1/2)*b0*(-1*h)**2 + (1/2)*b2*(h)**2 , 0) eq4 = Eq((1/2)*a0*(-1*h)**2 + (1/2)*a1*(h)**2 + (1/6)*b0*(-1*h)**3 + (1/6)*b2*(h)**3, 0) eq5 = ...

원문링크 : [Python] 연립방정식 풀기