This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
soddy(amax) = {
my(L=List(), abc, t, u);
for(a=1, amax,
for(b=1, a-1,
for(c=1, b-1,
abc=a*b*c;
if(issquare(abc*(a+b+c), &t),
u=a*b+a*c+b*c;
if(abc%(u+2*t) == 0,
if(u-2*t != 0,
if(abc%(u-2*t) == 0,
listput(L, [a,b,c,abc\(u+2*t),-abc\(u-2*t)])
)
,
listput(L, [a,b,c,abc\(u+2*t),0])
)
)
)
)
)
);
Vec(L)
}
soddy(253)
soddy(k) = {
s=[];
for(a=1, k,
for(b=1, a-1,
for(c=1, b-1,
if(issquare(a*b*c*(a+b+c), &t),
if(a*b*c % (a*b+a*c+b*c+2*t) == 0,
s=concat(s, a)
)
)
)
)
);
s
}
soddy(500)
Comments