I am trying to compute a limit with sage, and I get incorrect answers.
Here is an example which should produce the result of 0 but does not.
reset()
n = var('n')
assume(n>0)
series = -(3*n^2 + 1)*(-1)^n/sqrt(n^5 + 8*n^3 + 8)
working_series = -(7*n^2 + 8)*(-1)^n/sqrt(n^5 + 4*n^2 + 2)
print "Limit of series = "+str(limit(series, n=infinity))
print "Limit of working_series = "+str(limit(working_series, n=infinity))
This yields these results:
Limit of series = 38/25*pi^2*und
Limit of working_series = 0
Both Limits should be equal to zero. Removing the assumption of n>0 or
adding assume(n,'integer') does not change the outcome.
If someone is aware of why this is happening and/or knows a method to
solving this; please let me know.
Thank You,
Pavel Y.
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-support.
For more options, visit https://groups.google.com/groups/opt_out.
