A question about A079490
Pfoertner, Hugo
Hugo.Pfoertner at muc.mtu.de
Wed Jan 22 14:49:12 CET 2003
Mark,
my first guess for the Pi^n near integers sequence is
C PI^n near to Integers
INTEGER ARG
REAL*16 QARG, EX, DIFF, DIFMIN, PI
PI = QATAN2 (0.0Q0,-1.0Q0)
WRITE (*,*) ' PI=', PI
DIFMIN = 100.0Q0
DO 10 ARG = 1, 65
QARG = QEXT ( ARG )
EX = PI**QARG
DIFF = EX - QNINT(EX)
IF ( QABS(DIFF) .LT. DIFMIN ) THEN
DIFMIN = QABS(DIFF)
WRITE (*,*) ARG, DIFF
ENDIF
10 CONTINUE
END
with result:
PI= 3.14159265358979323846264338327950
1 0.141592653589793238462643383279503
2 -0.130395598910641381165509000123849
3 6.276680299820175476315067101392597E-0003
58 -5.378723144531250000000000000000000E-0003
Good luck with more terms...
Hugo
-----Ursprüngliche Nachricht-----
Von: Mark Hudson [mailto:mrmarkhudson at hotmail.com]
Gesendet am: 22 January, 2003 14:30
An: Pfoertner, Hugo
Betreff: Re: AW: A question about A079490
Thanks Hugo.
I believed that I had done it correctly, but I feel that it is best to come
at things from the angle of "What have I done wrong".
I'm currently going through the same procedure but using Pi (as in the
3.14159... and not the number of primes function) instead of exp().
I don't even know if it is in the database already!
Ho hum... my real work is so boring by comparison.
Thanks,
Mark.
More information about the SeqFan
mailing list