Hello,
I want to do some calculations over an extension of the polynomial ring
Z[q] which many quantum integers have square roots.
Initially I thought that I was going to be easy and that I could do the
following:
sage: R.q=PolynomialRing(QQ)
sage: s3=(1+q+q^2).sqrt(name='s3')
sage:
On 2014-10-09 10:31, Andrew wrote:
sage:R.q=PolynomialRing(QQ)
sage:s3=(1+q+q^2).sqrt(name='s3')
sage:s4=(1+q+q^2+q^3).sqrt(name='s4')
sage:s3*s4 # this blows up:(
If you don't mind approximations, you could consider using
PowerSeriesRing(QQ) instead. There, every polynomial with constant