Let's consider the following functions:

1/x,1/(exp(x)-1),1/sqrt(x)-all have a singularity at x=0.All got almost similar plot.Now,making the decision if the integration blows up or not at x=0,just by looking at the plot could be misleading.So after integration from 0 to a finite value,say 1, Maple provides the answer-the first two diverge,the last one converges to 2.Now how can I be certain that Maple just simply didn't make a mistake about the first two?Is there any...

In the following worksheet I want the plot, but Maple can't generate it . I tried two ways-the first one took almost 40 minutes and gave an empty plot, the second one took almost 10 minutes and generate an incomplete plot.Also the integration is very slow,but still,it can do it.

So,can anybody help me with a better way to do things-using other numerical schemes/advanced plotting options-whatever it takes?    Or at least explain what is going on?

By...

I'm running this numerical integration 14 hours continuously,without result.

The problem is the numerical integration

 evalf(int(F(q)*sin(q*Theta)/q)) assuming(0

This is taking forever.

Now,this is already done using Fortran (not time consuming at all),but the choice is not mandatory,so I'm using Maple 13.

Now,I think, to do this I have to do the integral using any other numerical method, evalf(int(,)) isn't good enough.

I need to send a Maple worksheet (.mw) to someone who is not a Maple user,doesn't have Maple installed.So what is the simplest possible solution now?How can I convert this into a pdf/doc/ps/eps  file so that any document viewer can open it?

Sorry to bother you again,fellow primers.A few days ago I asked a similar kind of q.,but here is some change, & I'm stuck again.

int(p^(D-2)/(exp(beta*p)-1), p = 0 .. infinity) assuming beta::positive, D::posint #D>0 & is an integer,this two are the only restriction on D,nothing else.

 I'm not getting any result (correct/incorrect).

The  result is  GAMMA(D-2)*Zeta(D-2)/beta^(D-2).  #Zeta is the Riemann zeta function.