Is there a way to find out the rate of decay of function f(x) when x-> +infinity? Hi all, I just wanted to use this real-world engineering problem to practice problem-solving using Maple. Suppose I have a function f(x) (probably complex valued), how to find out its rate of decay as x approaches +infinity? I used "limit" but it will give me a value, not a functional form. I already knew my function has a limit =0 at x=+infinity. And it is of fast decay. I would like to bound the limiting behavior of f(x) by exp(-alpha*x) for some suitable constant alpha. Is there a way to handle this in Maple?

Hi Folks, Here is a numerical integration program: http://people.scs.fsu.edu/~burkardt/f_src/quadpack/quadpack.f90 In evaluating the integrand functions, it is good to do two function evaluations at the same time simultaneously. I just found that I have two CPUs Pentium Xeon 2.4GHz, each supporing only up to SSE2. Although the CPUs are slow, but if I can utilize them to do parallel computing of the numerical integration. It would be great! But I am not a Fortran expert, nor am I a CS major... Could anybody give some easy-to-follow advice on how to make that code parallel on two CPUs?

Hi all, I am having lots of trouble using C codegen in generating C code from Maple expressions. My questions are: 1. How to assume each element of a whole array or vector >0? Of course I can do: Assume (a[1]>0, a[2]>0, a[3]>0)... But I have "m" elements, how can I do that for all m elements? 2. What to do if the output says "the output length exceeds maximum limit of 1000000"? This occurred when I tried to "simplify" some complicated expression... obviously Maple did get some results, but it refused to print them out due to length limit, ... 3. There is no way to translate the result of "JJ16:=optimize(JJ15, 'tryhard')"?

Hi all, I am asking this for my engineering friend. We've posted on Maple newsgroup and we just heard that here is another place where Maple experts appear, so we posted it again here to seek more advice -- because we have a struggling problem which annoyed us for several days already! Thank you! ---------------- Here is what I want to compute: t22 := Re((1/2*(a^2+2*sigma2*(c+g*(ee+i*v)))^(1/2)*exp(.5*a*t)/ (1/2*(a^2+2*sigma2*(c+g*(ee+i*v)))^(1/2)*cosh(1/2*(a^2+2*sigma2*(c +g*(ee+i*v)))^(1/2)*t)+.5*a*sinh(1/2*(a^2+2*sigma2*(c+g*(ee +i*v)))^(1/2)*t)))^(2*a*xbar/sigma2)) Assuming all variables are real and positive...

Let's say you want to compute Re(Gamma(a+b*I)) where a and b are real. Is there a way to obtain a good form of the real part and imaginary part?