Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #!/usr/bin/python3
- from sympy import *
- from fractions import Fraction
- a = var('a')
- p1 = -5/(6*a); q1 = Fraction(-2, 3) + (5/(6*a))**2
- p2 = a/4; q2 = Fraction(3, 2) + (a/4)**2
- f=((p1-p2)**2-(q1+q2))**2-4*q1*q2
- print(f.expand())
- g=Poly(f*a**2)
- print(roots(g))
Advertisement
