[seqfan] Re: A079039(8), decreasing trunc(cosh(n))
zak seidov
zakseidov at yahoo.com
Wed Jul 28 13:17:42 CEST 2010
1. Mmca code with 2000-digit accuracy gives:
(assuming that fractional part of cos(0)=1
as in A079039)
1, 3, 22, 29, 45, 75, 135, 259, 863, 1786, 2483, 2538
2. Actually cos(0)=1 and fractional part of cos(0) is zero not one. Hene defn in A079039'd be clarified, right?
--- On Wed, 7/28/10, Richard Mathar <mathar at strw.leidenuniv.nl> wrote:
> From: Richard Mathar <mathar at strw.leidenuniv.nl>
> Subject: [seqfan] A079039(8), decreasing trunc(cosh(n))
> To: seqfan at seqfan.eu
> Cc: benoit7848c at orange.fr
> Date: Wednesday, July 28, 2010, 6:21 AM
>
> I'd argue that the a(8) in www.oeis.org/classic/A079039 is
> 275 (not 863)
> and that the next term a(9) > 33500, right?
> This seems to be supported by the PARI code in the
> database, and a by
> a Maple implementation like
>
> Digits :=120 ;
> A079039 := proc(n)
> option remember;
> if n = 0 then
>
> return 1 ;
> end if;
> k := 1;
> prevmin :=
> frac(cosh(procname(n-1))) ;
> k := n+ 1;
> while
> evalf(frac(cosh(k))-prevmin) >= 0 do
> k
> := k+1 ;
> if
> k mod 500 = 0 then
>
> print("...",k) ;
> end
> if;
> end do;
> k ;
> end proc:
> for n from 0 do
> print( A079039(n)) ;
> end do;
>
>
>
> _______________________________________________
>
> Seqfan Mailing list - http://list.seqfan.eu/
>
More information about the SeqFan
mailing list