Hi, I have done this calculation using a very tedious way and have checked that it is correct. Can I possibly perform this using Sage only.
I have a polynomial : $$D=M^2-(A/(2p^4))M+(B/(16p^4))$$ where $$A=18p^{10} - 54p^9 + 59p^8 + 130p^7 - 209p^6 - 98p^5 + 407p^4 + 362p^3 + 49p^2 - 16p + 8$$ and $$B=9( p + 1 )^2(p^4 - 2p^3 + 2p^2 + 2p + 1)(4p^8 - 52p^7 + 373p^6 + 68p^5 - 445p^4 + 72p^3 + 163p^2 - 48*p + 9)$$.
I have checked using multiple software that the factorization of D using the quartic in $p$ $$v^2= p^4-2p^3+5p^2+8p+4$$ gives $$[M-((A+2Fv)/(4p^4))][M-((A-2Fv)/(4p^4))]$$ where $$F=9p^8-18p^7-7p^6+45p^5-21p^4-74p^3-18p^2+6p-2$$.
Can someone help me obtain the same result by using Sage only. Thank you.