MathGroup Archive 1999

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Real integrand->complex result.

  • To: mathgroup at
  • Subject: [mg20310] Re: Real integrand->complex result.
  • From: Jens-Peer Kuska <kuska at>
  • Date: Fri, 15 Oct 1999 20:20:41 -0400
  • Organization: Universitaet Leipzig
  • References: <7tuou7$>
  • Sender: owner-wri-mathgroup at

Hi William,

a) please use InputForm[] to paste formulas into e-mails or append a

b) you should avoid mixing floting point numbers and analytical
   terms. Introducing omega instead of 3.66558239083868908` and using
   you will see that all complex terms are removed

c) what do you think means 0.`\ \[ImaginaryI] ? something like
   0*I or 0? You can use Chop[] to remove this zero part of the complex
d) the evaluation is to slow ? There are several ways to avoid this. You
   buy a new faster computer. You may wait for a new faster computer and
a new
   faster Mathematica. You may put some ideas into you program. The 
   typical case is that you can find some kind of recursion for your 
   simpFphi[n], than you can use dynamic programming and store
   to reduce the computation time. You may also be able to find some
   for Integrate[simpFphi[n],{x,0,l}]= someterm[n,l]
   and avoid the integration at all. 

Hope that helps


William Golz wrote:
> Synopsis of the problem:
> I get a complex result of the form f(t)=(a+bi)g(t) when evaluating a
> real integrand F(x,t)*phi(x) over the interval {x:0<x<L). All of the
> variables have been properly declared as real and the functions have
> been declared as real for real arguments (using the package "ReIm").

> The evaluation of the integral is also quite slow, even for n=2, even
> when the expression is simplified as much as possible, and the final
> program will need to evaluate a larger number of terms.
> Thanks in advance, Bill

  • Prev by Date: Re: Re: Replacing Part of a Matrix
  • Next by Date: My brackets/Parentheses are getting translated to text. Help!
  • Previous by thread: Re: Real integrand->complex result.
  • Next by thread: Programming the front end