Am getting terms with expressions like:
7.10000000000000undef + 3.98967171417774e277x
That is 7.10000000000000 times undef
what does undef mean in sage math and how can I simplify the expression or remove the undef?
EXAMPLE:
beta= 7.10000000000000
y = 2.66778786276411e493*t^3*x - 1.93258971534823e247*x
result=beta*(integral(y, t,0,1))
print('value is:',result)
integral(y, t,0,1)
This is a bug in SageMath. I've reported it here: https://github.com/sagemath/sage/issu...
Here is a workaround, specific to polynomials:
var('x,t')
beta = 7.10000000000000
y = 2.66778786276411e493*t^3*x - 1.93258971534823e247*x
Y = SR(y.polynomial(RR).integral(t))
result = beta*(Y.subs(t=1) - Y.subs(t=0))
Asked: 2024-05-05 10:08:47 +0100
Seen: 221 times
Last updated: May 06
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
I haven't seen this personally so please post a (minimal) reproducible example of SageMath code producing such a result, so it can be properly investigated.
Interesting, the huge coefficients break every algorithm used by Sage:
integral(y, t, 0, 1, algorithm=a)breaks for eachain['maxima', 'sympy', 'mathematica_free', 'fricas', 'giac', 'libgiac']. Theundefcomes fromlibgiac. The error from choosing'sympy'looks like it should be fixable (on the developer side, not the user side).Thanks for investigating. I see the point. As an alternative is there a way to round the expression? round(y,10) <- something like that