Risposta:
# (X- (1 + sqrt (5)) a / 2) (x- (1-sqrt (5)) a / 2) #
# (X + (3 + sqrt (5)) a / 2) (x- (sqrt (5) -3) y / 2) = 0 #
Spiegazione:
# "Risolvi l'equazione quartica caratteristica senza prima la y:" #
# x ^ 4 + 2 x ^ 3 - 3 x ^ 2 - 4x - 1 = 0 #
# => (x ^ 2-x-1) (x ^ 2 + 3x + 1) = 0 "(*)" #
# "1)" x ^ 2 + 3x + 1 = 0 => x = (-3 pm sqrt (5)) / 2 #
# "2)" x ^ 2-x-1 = 0 => x = (1 pm sqrt (5)) / 2 #
# "Se applichiamo questo sul polinomio dato otteniamo" #
# (x ^ 2 - x y - y ^ 2) (x ^ 2 + 3 xy + y ^ 2) = 0 #
# => (x- (1 + sqrt (5)) y / 2) (x- (1-sqrt (5)) y / 2) #
# (X + (3 + sqrt (5)) a / 2) (x- (sqrt (5) -3) y / 2) = 0 #
# "(*) Con la sostituzione" x = y-1/2 "otteniamo:" #
# y ^ 4 - (9/2) y ^ 2 + 1/16 = 0 #
# "Ora inserisci" z = y ^ 2 "e moltiplica con 16:" #
# 16 z ^ 2 - 72 z + 1 = 0 #
# "disco:" 72 ^ 2 - 4 * 16 = 5120 = 32 ^ 2 * 5 #
# => z = (72 pm 32 sqrt (5)) / 32 = 9/4 pm sqrt (5) #
# => y = pm sqrt (9/4 pm sqrt (5)) #
# => y = pm sqrt (9 pm 4 sqrt (5)) / 2 #
# => y = pm sqrt ((2 pm sqrt (5)) ^ 2) / 2 #
# => y = pm (1 pm sqrt (5) / 2) #
# => x = (1 pm sqrt (5)) / 2 "or" (-3 pm sqrt (5)) / 2 #