from sympy import *

expand((1+sqrt(5))**5)