Hi community.
I'm interested in manipulating a formal expression
$$\int \mathrm{d}t \; e^{- 2 \int \mathrm{d}t \; h(t)}.$$
My notebook contains the following code
var('t')
h = function('h')(t)
integrate(exp(integrate(-2*h,t)),t)
but the result is +Infinity.
QUESTION(S):
h is not explicit?+Infinity come from?Thanks to the comment by @fredericc I found the answer! The solution is to use the algorithm='maxima' flag, e.g.
var('t')
h = function('h')(t)
integrate(exp(integrate(-2*h,t)),t, algorithm='maxima')
With the above, I can also manipulate the expression... differentiate or other:
var('t')
h = function('h')(t)
HH = function('HH')(t)
HH = exp(integrate(2*h,t))*(Cg + k*integrate(exp(integrate(-2*h,t)),t, algorithm='maxima'))
diff(HH,t)
The above result is compatible with the by-hand calculation and the result by Mathematica
Asked: 2022-10-09 10:34:55 +0100
Seen: 241 times
Last updated: Oct 10 '22
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.
The culprit is giac:
Interestingly :
Enhanceable ?
Indeed, Giac misbehaves.
And :
Maxima is correct.
Would you mind filing a ticket ?