home *** CD-ROM | disk | FTP | other *** search
- Path: sparky!uunet!cs.utexas.edu!sun-barr!ames!data.nas.nasa.gov!mustang.mst6.lanl.gov!nntp-server.caltech.edu!harpo.tapir.Caltech.EDU!rauch
- From: rauch@harpo.tapir.Caltech.EDU (Kevin Rauch)
- Newsgroups: sci.math
- Subject: Re: 1+1/2+1/3+1/4+...+1/n
- Date: 20 Nov 1992 01:18:56 GMT
- Organization: California Institute of Technology, Pasadena
- Lines: 30
- Message-ID: <1ehee0INN9i0@gap.caltech.edu>
- References: <92324.132329K3032E2@ALIJKU11.BITNET>
- NNTP-Posting-Host: harpo.tapir.caltech.edu
-
- Mutter Christoph Johannes <K3032E2@ALIJKU11.BITNET> writes:
-
- >Hallo.
- >I've a problem. I have to calculate the sum 1+1/2+1/3+1/4+...+1/n.
- >The result should be 100. But my computer is far too slow, to solve this
- >problem. I need the index n when the sum > 100.0
- >And that exactly.
- >Chris. ------------
-
-
- Use the identity:
-
- 1+1/2+ ... +1/n = gamma + psi(1+n) ,
-
- where gamma=0.577... is the Euler gamma constant and psi(z) is the
- digamma function (logarithmic derivative of the gamma function). Using
- Mathematica I get the sum to be >100 for n greater than or equal to:
-
- 15092688622113788323693563264538101449859497
-
- (about 1.5e43; the estimate ceiling(exp(100)) is off by 1, which is good enough
- all the same, I think! :)
-
-
- &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
- & Kevin Rauch, Theoretical Astrophysics, California Institute of Technology &
- & email: rauch@tapir.caltech.edu &
- & "Got my sights on the stars. Won't get that far, but I'll try anyway..." &
- &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& --N. P. &&&&&&&&&
- &&&&&&&&&&&&&
-