MathGroup Archive 2007

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

Search the Archive

Re: Integration using mathematica

  • To: mathgroup at smc.vnet.net
  • Subject: [mg72640] Re: Integration using mathematica
  • From: "dimitris" <dimmechan at yahoo.com>
  • Date: Sat, 13 Jan 2007 04:30:02 -0500 (EST)
  • References: <eo7mt2$1ps$1@smc.vnet.net>

First things first!

In order to use Mathematica you should first learn the Basic Syntax of
it!

Here is a good place to start:
Copy/paste the following command in a notebook and execute it!

FrontEndExecute[{HelpBrowserLookup["MainBook", "T.0"]}]

You are right for example about the first letter being capitalized; but
on the other hand arguments of functions must be enclosed on square
brackets.

Apart from this your input misses some parentheses and I was not able
to figure out
exactly your expression!

Suppose I am right and I guess well, let f be the following expression

f = 0.013*Exp[15.16*ArcTan[0.002*x]^2]*ArcTan[0.002*x];

One thing you should learn for your very begining in using Mathematica
is not use inexact numbers for symbolic calculations like integration.
You can consult the archives of this forum if you want more details
about why.

A quick way to convert your real numbers to exact Rational numbers is
by using Rationalize.
Then

ff=Rationalize[f]
(13*E^((379/25)*ArcTan[x/500]^2)*ArcTan[x/500])/1000

Let's try to integrate (indefinite integral) this exression w.r.t. x.
We get

Integrate[ff, x]
(13*Integrate[E^((379/25)*ArcTan[x/500]^2)*ArcTan[x/500], x])/1000

Mathematica fails to return a analytical result.
There are three alternative explanations:

1) There is not a closed form expression. Especially for indefinite
integrals usually this is the
case within Mathematica.
2) Mathematica does not know the integral. Seldom especially for
indefinte integrals.
3) Mathematica "needs" a little help in order to determine the
integral.

Let know try definite integral of this expression. There are many
occasions that there is not a closed form expression for an indefinite
integral but there are closed form expressions for definite integrals.
For example

Integrate[BesselJ[0, x]/(1 + x^2), x]
Integrate[BesselJ[0, x]/(1 + x^2), x]

However,

Integrate[BesselJ[0, x]/(1 + x^2), {x, 0, Infinity}]
(1/2)*Pi*(BesselI[0, 1] - StruveL[0, 1])

Moreover, in definite integration you have also the numerical "way".

A quick check of the last result:

N[%]
0.8730842426508671

NIntegrate[BesselJ[0, x]/(1 + x^2), {x, 0, Infinity}, Method ->
Oscillatory]
0.8730842426508675

So, let integrate definite ff, say between 0 and Pi.
Before performing a definite integration (symbolic or numerical) always
plot the integrand for possible singularities in the intagrtion range.

Plot[(13*E^((379/25)*ArcTan[x/500]^2)*ArcTan[x/500])/1000, {x, 0, Pi},
PlotPoints -> 1000]

The graph, show normal behavior of the integrand. Warning, howver, that
the graph does
not guarranty that!

Anyway,

Integrate[ff, {x, 0, Pi}]
Integrate[(13*E^((379/25)*ArcTan[x/500]^2)*ArcTan[x/500])/1000, {x, 0,
Pi}]

Again no closed form result!

Finally,

NIntegrate[ff, {x, 0, Pi}]
0.00012834241443264268

I hope you will find my reply helfpful.
Try to repost your message with correct syntax if my guess was not
right!

Regards
Dimitris


Negede wrote:
> I have difficulty in integrating the following equation on Mathematica.
> Does any one of you know how to do it on Mathematica? Any trick to
> aply? or another software to use? Please help me. The function is as
> follows:
>
> 0.013 *Exp [15.16(ArcTan(0.002 x)])^2 *(ArcTan(0.002 x)
> 
> Thanks in advance


  • Prev by Date: Re: Integration using mathematica
  • Next by Date: Re: An Introduction to Programming with Mathematica, Third Edition
  • Previous by thread: On Matematica notebook options
  • Next by thread: On AspectRatio