[Maxima] plotting a parametric integral
l.butler at ed.ac.uk
Fri Oct 9 11:19:04 CDT 2009
On Fri, 9 Oct 2009, Leo Butler wrote:
< On Fri, 9 Oct 2009, Stanislav Maslovski wrote:
< < Hello,
< < I am trying to plot this simple Fourier integral in Maxima 5.17.1:
< (%i2) f(t) ::= buildq([k:t],
< (%o2) f(t)::=buildq([k:t],integrate(exp(%i*k*x)*sin(x)/x,x,minf,inf))
< This will avoid any calls to 'asksign' (which is asking about the sign
< of k) when t is a constant. I have also set ratprint to false to avoid
< the annoying warnings about floats being converted to rationals.
< (%i3) l : makelist(i/10,i,-20,20)$
< Note that you can also use 'assume':
< (%i35) [ (map(forget,facts()), assume(k<-1),f(k)) , (forget(k<-1),
< assume(k>1),f(k)) , (forget(k>1), assume(k>-1,k<1),f(k)), f(-1), f(1) ];
< (%o35) [0,0,%pi,0,%pi]
I should say that there are two errors in the output here. f(-1) should
equal f(1) and both should be
Interestingly, Maple does not give the correct answer for f(-1) or f(1)
either, and Wolfram Alpha times out in doing the integral; both do the
trig integral correctly.
I have filed a bug report about this.
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.
More information about the Maxima