You are not logged in.

- Topics: Active | Unanswered

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

Any other way? This one seems pretty impractical.

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

**bobbym****Administrator**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 103,828

Hmmm, impractical? In 135 years of working with that rule it failed only once! Now that may not be good enough for a load of mathematicians like A that never calculated anything but it is a workhorse for us.

Sometimes we can do a tail analysis which can help us know how close we are. You are aware of the methods of numerical integration? Know some of the theory?

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**

**Online**

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

I'm not sure. I have heard of a few formulas, including the trapezoidal, simpsons, romberg,... Of these, I know only what the trapezoidal is. I am willing to learn it though.

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

**bobbym****Administrator**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 103,828

You know of a simpler one than those. You know about stuffing rectangles to get an area?

By the way, when you say things like you said in the previous post I am reminded of what pappym told me.

bobbym, everything you know is wrong! - pappym

http://www.youtube.com/watch?v=z4jeREy7Pbc

Everything they taught you is wrong. Think for yourself, do not opine the words of others. The double method is highly practical and we can improve it.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**

**Online**

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

Oh, yes, Riemann integration. I also remember Euler-Mclaurin...

So, what kind of tail analysis do you have in mind?

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

**bobbym****Administrator**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 103,828

The simplest one, If areas ( integrals ) can be approximated by rectangles then rectangles can be approximated by integrals. Each 1 x n rectangle can be thought of as a term in the sum.

Testing this we see that we are off by

by actual computation. There are some conclusions, notice that the double rule was predicting 8 and we really had about 10. Showing that it is on the conservative side and can be trusted.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**

**Online**

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

Ah, that one's nice.

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

**bobbym****Administrator**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 103,828

With it you can answer your question of how many terms you would need for 100 digits of accuracy.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**

**Online**

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

10^50 with the original series! That is really a lot!!

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

**bobbym****Administrator**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 103,828

For that many digits none of the methods looked at so far will suffice.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**

**Online**

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

I know. It needs exponential convergence or some trick, I guess...

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

**bobbym****Administrator**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 103,828

We will get into a few sequence converters that might help tomorrow.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**

**Online**

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

Can't wait. Tomorrow, then!

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

**bobbym****Administrator**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 103,828

I did a little preliminary work. It was not very encouraging.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**

**Online**

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

So, those sequence converters. What do they look like?

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

**bobbym****Administrator**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 103,828

I use a couple of them. Shanks, Romberg, Aitken, Euler, RRA, and a new one.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**

**Online**

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

Well, how does the first one work?

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

**bobbym****Administrator**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 103,828

Romberg acceleration worked best so far. We will start there.

Romberg integrates by taking a sequence of values and using Richardson extrapolation on them. This means it can accelerate many sequences. Have you ever programmed a Romberg integration?

If you have not and do not know how to write it I will give you mine.

See you a little later, chores are calling.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**

**Online**

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

I haven't programmed it, but I remember you helping me find it among the Maxima packages and use it for the 1/(cos x+x^2) integral.

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

**bobbym****Administrator**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 103,828

Hi;

I can give you mine if you want and we can go on.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**

**Online**

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

That would be okay.

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

**bobbym****Administrator**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 103,828

Hi;

```
romberg[l_] := Module[{x, a, o, l1}, a = Dimensions[l][[1]]; x = 1/2;
o = Table[0, {a}, {a - 1}];
l1 = augment[Transpose[{l}], o];
For[i = 2, i <= a, ++i,
For[j = 2, j <= a, ++j,
l1[[i, j]] = (l1[[i, j - 1]] -
x^(j - 1)*l1[[i - 1, j - 1]])/(1 - x^(j - 1));];];
l1[[a, a]]];
augment[a_, b_] := Transpose[Join[Transpose[a], Transpose[b]]];
```

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**

**Online**

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

Does it work?

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

**bobbym****Administrator**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 103,828

Yes, try it on this simple one:

romberg[{1, 5/4, 205/144, 1077749/705600,

822968714749/519437318400}]

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**

**Online**

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

Isn't it better to use recursion?

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