Math Is Fun Forum

  Discussion about math, puzzles, games and fun.   Useful symbols: ÷ × ½ √ ∞ ≠ ≤ ≥ ≈ ⇒ ± ∈ Δ θ ∴ ∑ ∫ • π ƒ -¹ ² ³ °

You are not logged in.

#1 2013-02-02 21:21:23

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Kummer and a tough integral.

Hi;

The evaluation numerically of this integral came up in a discussion with anonimnystefy:

This one can resist normal methods quite well. When we say evaluate numerically we typically mean 10 - 12 digits with a programmable calculator and at least 100 when using a CAS.

To start, although any CAS known will gag on this problem they are able to do pieces of so it is to our advantage to split the integral like this.

The first integral is not too difficult and any CAS can do it to 100 digits.

That is about 105 digits or so. This is the hard part:

we notice that asymptotically as x grows large the denominator will act like x^2. The x^2 term will completely drown out the cos(x). That suggests the following chain of moves,

Looking at the RHS we see the same idea works again. As x grows large the denominator acts like x^4 completely drowning out the x^2 cos(x) term. So let's add a cos(x) / x^4 to it, we get

As x grows larger the denominator acts like x^6 completely drowning out the x^4 cos(x) term. So let's subtract a cos(x)^2 / x^6 to it, we get

Two things should be apparent,

1) we can continue the series on the left indefinitely and the RHS is getting smaller and smaller for each term on the left. The RHS will eventually be 0 as we get an infinite series on the LHS.

This becomes:

What have we accomplished? We have replaced a tough integral with a sum of infinite number of hopefully easier to evaluate numerically integrals. The beauty is that although we have an infinite number of integrals we will only need a finite number of them for the accuracy required.
Actually we will only need the first 26 or so terms.

Now I am not proposing this method of evaluating that integral because in practice those 26 integrals can also pose problems and further work is needed. But you wanted to see how the acceleration works.


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

#2 2013-02-03 00:03:02

anonimnystefy
Real Member
From: Harlan's World
Registered: 2011-05-23
Posts: 16,049

Re: Kummer and a tough integral.

How many terms would you need for 1000 digits?


“Here lies the reader who will never open this book. He is forever dead.
“Taking a new step, uttering a new word, is what people fear most.” ― Fyodor Dostoyevsky, Crime and Punishment
The knowledge of some things as a function of age is a delta function.

Offline

#3 2013-02-03 01:41:13

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Re: Kummer and a tough integral.

Hi;

I estimate about 255 terms. As I said this method is too weak for that many digits. I just showed it because the transformation is one of the weapons used in numerical work.


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

#4 2013-02-03 01:46:42

anonimnystefy
Real Member
From: Harlan's World
Registered: 2011-05-23
Posts: 16,049

Re: Kummer and a tough integral.

Do sum accelerators work on it (like RRA and stuff...) ?


“Here lies the reader who will never open this book. He is forever dead.
“Taking a new step, uttering a new word, is what people fear most.” ― Fyodor Dostoyevsky, Crime and Punishment
The knowledge of some things as a function of age is a delta function.

Offline

#5 2013-02-03 01:53:24

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Re: Kummer and a tough integral.

I do not know the answer to that. I gave up on it when the 26 integrals were more difficult than they should have been.

Want to see another solution?


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

#6 2013-02-03 01:58:35

anonimnystefy
Real Member
From: Harlan's World
Registered: 2011-05-23
Posts: 16,049

Re: Kummer and a tough integral.

Of course!


“Here lies the reader who will never open this book. He is forever dead.
“Taking a new step, uttering a new word, is what people fear most.” ― Fyodor Dostoyevsky, Crime and Punishment
The knowledge of some things as a function of age is a delta function.

Offline

#7 2013-02-03 02:11:26

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Re: Kummer and a tough integral.

This method was showed to me.

Divide the integration region into lengths of
2 π and use cos(x) = cos(x + n*2π) having n as an integer. Then

now according to him we can interchange the order of the sum and the integral to get:

Mathematica can get the sum.

Now you just numerically integrate


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

#8 2013-02-03 02:28:40

anonimnystefy
Real Member
From: Harlan's World
Registered: 2011-05-23
Posts: 16,049

Re: Kummer and a tough integral.

Have you been able to get 1000 digits out of that?


“Here lies the reader who will never open this book. He is forever dead.
“Taking a new step, uttering a new word, is what people fear most.” ― Fyodor Dostoyevsky, Crime and Punishment
The knowledge of some things as a function of age is a delta function.

Offline

#9 2013-02-03 02:32:08

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Re: Kummer and a tough integral.

I have not tried!


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

#10 2013-02-03 02:33:20

anonimnystefy
Real Member
From: Harlan's World
Registered: 2011-05-23
Posts: 16,049

Re: Kummer and a tough integral.

Why not?

And, it seems to me that the point is to get the integral into a form with finite integration limits...


“Here lies the reader who will never open this book. He is forever dead.
“Taking a new step, uttering a new word, is what people fear most.” ― Fyodor Dostoyevsky, Crime and Punishment
The knowledge of some things as a function of age is a delta function.

Offline

#11 2013-02-03 02:35:43

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Re: Kummer and a tough integral.

For several reasons I do not like that solution.

1) I do not understand it!

2) Different versions of Mathematica were unable to do the sum.

3) I suspect it would take too long.


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

#12 2013-02-03 02:40:39

anonimnystefy
Real Member
From: Harlan's World
Registered: 2011-05-23
Posts: 16,049

Re: Kummer and a tough integral.

What do you not understand?

What can Mathematica do with the integral

Last edited by anonimnystefy (2013-02-03 02:40:57)


“Here lies the reader who will never open this book. He is forever dead.
“Taking a new step, uttering a new word, is what people fear most.” ― Fyodor Dostoyevsky, Crime and Punishment
The knowledge of some things as a function of age is a delta function.

Offline

#13 2013-02-03 02:55:19

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Re: Kummer and a tough integral.

The first line, beats me.

That will not work.


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

#14 2013-02-03 02:58:41

anonimnystefy
Real Member
From: Harlan's World
Registered: 2011-05-23
Posts: 16,049

Re: Kummer and a tough integral.

Well, if you do not understand it, there must be a part which you do not understand...

Why not?


“Here lies the reader who will never open this book. He is forever dead.
“Taking a new step, uttering a new word, is what people fear most.” ― Fyodor Dostoyevsky, Crime and Punishment
The knowledge of some things as a function of age is a delta function.

Offline

#15 2013-02-03 03:02:38

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Re: Kummer and a tough integral.

The trig substitution that produced the series, I do not get.

I tried a million substitutions to get the interval from 0 to infinity down to 0 to 1. Each one failed. Examine your integrand around 1.


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

#16 2013-02-03 03:07:30

anonimnystefy
Real Member
From: Harlan's World
Registered: 2011-05-23
Posts: 16,049

Re: Kummer and a tough integral.

He just splits it up in intervals and uses cos(x+2n*pi)=cos(x). I do not see what the big deal is...

The integrand of my integral tends to 1 at 1.


“Here lies the reader who will never open this book. He is forever dead.
“Taking a new step, uttering a new word, is what people fear most.” ― Fyodor Dostoyevsky, Crime and Punishment
The knowledge of some things as a function of age is a delta function.

Offline

#17 2013-02-03 03:15:06

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Re: Kummer and a tough integral.

The cos(x) does not change, the x^2 does. Why?

There is a big fat indeterminate at x = 1 for your function. Numerical techniques will have problems with that.


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

#18 2013-02-03 03:38:04

anonimnystefy
Real Member
From: Harlan's World
Registered: 2011-05-23
Posts: 16,049

Re: Kummer and a tough integral.

Because cos(x+n*2pi)=cos(x).

Did you actually try entering the integral into Mathematica?


“Here lies the reader who will never open this book. He is forever dead.
“Taking a new step, uttering a new word, is what people fear most.” ― Fyodor Dostoyevsky, Crime and Punishment
The knowledge of some things as a function of age is a delta function.

Offline

#19 2013-02-03 03:41:49

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Re: Kummer and a tough integral.

Hi;

Because cos(x+n*2pi)=cos(x).

I am not getting that.

Yes, I entered it into the built in command and one I wrote. Both could not deal with it.


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

#20 2013-02-03 03:45:01

anonimnystefy
Real Member
From: Harlan's World
Registered: 2011-05-23
Posts: 16,049

Re: Kummer and a tough integral.

Cosine has a period of 2pi, that's why how ever many times you add 2pi to the argument of cosine you will get the same thing.

Interestingly enough, WA gives 5 digits for that integral...


“Here lies the reader who will never open this book. He is forever dead.
“Taking a new step, uttering a new word, is what people fear most.” ― Fyodor Dostoyevsky, Crime and Punishment
The knowledge of some things as a function of age is a delta function.

Offline

#21 2013-02-03 03:47:19

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Re: Kummer and a tough integral.

I know that one but I still do not undertsand what he did.

Probably that is all the digits could get. That integral is very difficult to compute.


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

#22 2013-02-03 03:54:28

anonimnystefy
Real Member
From: Harlan's World
Registered: 2011-05-23
Posts: 16,049

Re: Kummer and a tough integral.

Hi

What do you get when you substitute t=x+n*2pi into the integrals?


“Here lies the reader who will never open this book. He is forever dead.
“Taking a new step, uttering a new word, is what people fear most.” ― Fyodor Dostoyevsky, Crime and Punishment
The knowledge of some things as a function of age is a delta function.

Offline

#23 2013-02-03 04:12:24

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Re: Kummer and a tough integral.

I get a big mess:

I left out the intervals of integration because they remain the same.


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

#24 2013-02-03 04:15:48

anonimnystefy
Real Member
From: Harlan's World
Registered: 2011-05-23
Posts: 16,049

Re: Kummer and a tough integral.

What do you get as the differential of t, i.e. as dt?


“Here lies the reader who will never open this book. He is forever dead.
“Taking a new step, uttering a new word, is what people fear most.” ― Fyodor Dostoyevsky, Crime and Punishment
The knowledge of some things as a function of age is a delta function.

Offline

#25 2013-02-03 04:17:24

bobbym
bumpkin
From: Bumpkinland
Registered: 2009-04-12
Posts: 109,606

Re: Kummer and a tough integral.

I got dt = 2 π n dx


In mathematics, you don't understand things. You just get used to them.
If it ain't broke, fix it until it is.
Always satisfy the Prime Directive of getting the right answer above all else.

Offline

Board footer

Powered by FluxBB