Harj. 4 loppuviikko, ratkaisuja

16.3.2000

Teht. 2

> A:=n->int(1/x^3,x=n..infinity);A(n);

> evalf(seq((1/2)*(A(n)-A(n+1)),n=21..30));

> va:=n->(1/2)*(A(n)-A(n+1)); # virhearvio(n)

> va(37.);

> S:=n->add(1/k^3,k=1..n); # Osasumma(n)

Warning, `k` in call to `add` is not local

> St:=n->S(n)+(1/2)*(A(n)+A(n+1)); # Stähti(n)

> St(37.);

> evalf(S(37.));

> s:=sum(1/n^3,n=1..infinity); # Maple osaa/tuntee

> s:=evalf(s);

> s-St(37.);

> s-S(37);

> s-S(230);

> s-S(1000);

>