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