Sarjat HA

L to 24.1.02

p-harmoniset sarjat

> seq(sum(1./n,n=1..m),m=1..100);

1., 1.500000000, 1.833333333, 2.083333333, 2.283333...
1., 1.500000000, 1.833333333, 2.083333333, 2.283333...
1., 1.500000000, 1.833333333, 2.083333333, 2.283333...
1., 1.500000000, 1.833333333, 2.083333333, 2.283333...
1., 1.500000000, 1.833333333, 2.083333333, 2.283333...
1., 1.500000000, 1.833333333, 2.083333333, 2.283333...
1., 1.500000000, 1.833333333, 2.083333333, 2.283333...
1., 1.500000000, 1.833333333, 2.083333333, 2.283333...
1., 1.500000000, 1.833333333, 2.083333333, 2.283333...
1., 1.500000000, 1.833333333, 2.083333333, 2.283333...

> log(100.);

4.605170186

> seq(evalf(log(m)),m=1..100);

0., .6931471806, 1.098612289, 1.386294361, 1.609437...
0., .6931471806, 1.098612289, 1.386294361, 1.609437...
0., .6931471806, 1.098612289, 1.386294361, 1.609437...
0., .6931471806, 1.098612289, 1.386294361, 1.609437...
0., .6931471806, 1.098612289, 1.386294361, 1.609437...
0., .6931471806, 1.098612289, 1.386294361, 1.609437...
0., .6931471806, 1.098612289, 1.386294361, 1.609437...
0., .6931471806, 1.098612289, 1.386294361, 1.609437...
0., .6931471806, 1.098612289, 1.386294361, 1.609437...
0., .6931471806, 1.098612289, 1.386294361, 1.609437...

> plot([seq([m,sum(1./n,n=1..m)],m=1..100)],style=point);

[Maple Plot]

> sum(1./n,n=1..1000);

7.485470861

>

> ithprime(10);

29

> ?prime

> seq(sum(1./ithprime(k),k=1..N),N=[100,200,300]);

2.106342121, 2.227059658, 2.290945735

> seq(sum(1./ithprime(k),k=1..N),N=[100,500]);

2.106342121, 2.365329466

>

> sum(1/n^2,n=1..infinity);

1/6*Pi^2

> sum(1/n^2,n=1..N);

-Psi(1,N+1)+1/6*Pi^2

> sum(1/n^3,n=1..N);

1/2*Psi(2,N+1)+Zeta(3)

>

Vuorottelevat sarjat

> seq(Sum((-1)^(k-1)*1/k,k=1..n),n=1..10);

Sum((-1)^(k-1)/k,k = 1 .. 1), Sum((-1)^(k-1)/k,k = ...
Sum((-1)^(k-1)/k,k = 1 .. 1), Sum((-1)^(k-1)/k,k = ...

> map(value,[%]);

[1, 1/2, 5/6, 7/12, 47/60, 37/60, 319/420, 533/840,...

> evalf(%);

[1., .5000000000, .8333333333, .5833333333, .783333...

> evalf(seq(sum((-1)^(k-1)*1/k,k=1..n),n=1..30));

1., .5000000000, .8333333333, .5833333333, .7833333...
1., .5000000000, .8333333333, .5833333333, .7833333...
1., .5000000000, .8333333333, .5833333333, .7833333...

> evalf(seq(sum((-1)^(k-1)*1/k,k=1..n),n=[1,3,5,7,9]));

1., .8333333333, .7833333333, .7595238095, .7456349...

> evalf(seq(sum((-1)^(k-1)*1/k,k=1..n),n=[2,4,6,8,10]));

.5000000000, .5833333333, .6166666667, .6345238095,...

> evalf(sum((-1)^(k-1)*1/k,k=1..100));

.6881721793

> evalf(sum((-1)^(k-1)*1/k,k=1..101));

>

.6980731694

>

Havaitaan:

Parittomat osasummat muodostavat pienenevän jonon.

Parilliset muodostavat kasvavan jonon.

Jokainen pariton > mikä tahansa parillinen

Suurilla n:n arvoilla näyttävät lähestyvän toisiaan (joskin hitaasti).