OFFSET
1,2
LINKS
David Broadhurst, Table of n for n = 1..1000
Kevin Acres and David Broadhurst, Rational points on y^2 = x^3 + 10*x^2 + 5*x
EXAMPLE
[1, 0.6931471805599453094172321215]
[2, 0.5994738181995926360154858945]
[3, 0.5898075219334671955890478209]
[5, 0.5866039232314788114510488867]
[10, 0.5865963134453746145925657102]
[13, 0.5864356249925781873672553233]
[18, 0.5864232678262123920533476064]
PROG
(PARI)
T(k)={
local(P=ellpow(ellinit([0, 10, 0, 5, 0]), [-1, 2], k), s, t);
s=P[1]^2;
t=abs(numerator(P[2]^4/s-80*s));
while(t%2==0, t=t/2; );
t;
}
M(k)=1+k%2;
V(k)=log(M(k)*T(k))/k^2;
{
lowest=1;
for(i=1, 20000,
l=V(i);
if(l<lowest,
lowest=l;
print([i, l]);
);
);
}
CROSSREFS
KEYWORD
nonn
AUTHOR
Kevin T. Acres, Feb 12 2011
STATUS
approved