AXIOM SESSIONscale=2.5]
FermatNumber(n:PositiveInteger):PositiveInteger == 2**2**n + 1
Type: Void
Time: 0 sec
l := [FermatNumber(i) for i in 1..6]
(2) [5,17,257,65537,4294967297,18446744073709551617]
Type: List PositiveInteger
Time: 0.01 (IN) + 0.01 (OT) = 0.02 sec
[prime?(n) for n in l]
(3) [true,true,true,true,false,false]
Type: List Boolean
Time: 0.01 (EV) + 0.02 (OT) = 0.03 sec
[factor(n) for n in l]
(4) [5,17,257,65537,641 6700417,274177 67280421310721]
Type: List Factored Integer
Time: 0.04 (EV) + 0.01 (OT) = 0.05 sec
-- factor FermatNumber(7)
[factorial(n) for n in 99..101]
(5)
[
93326215443944152681699238856266700490715968264381621468592963895217599993_
2299156089414639761565182862536979208272237582511852109168640000000000000_
000000000
,
93326215443944152681699238856266700490715968264381621468592963895217599993_
2299156089414639761565182862536979208272237582511852109168640000000000000_
00000000000
,
94259477598383594208516231244829367495623127947025437683278893534169775993_
1622147650308786159180834691162349000354959958336970630260326400000000000_
0000000000000
]
Type: List Integer
Time: 0.01 (EV) + 0.01 (OT) = 0.02 sec
S(n:PositiveInteger):Fraction Integer == reduce(+,[1/k**2 for k in 1..n])
Type: Void
Time: 0 sec
[S(n) for n in 1..20]
(7)
5 49 205 5269 5369 266681 1077749 9778141 1968329 239437889
[1, -, --, ---, ----, ----, ------, -------, -------, -------, ---------,
4 36 144 3600 3600 176400 705600 6350400 1270080 153679680
240505109 40799043101 40931552621 205234915681 822968714749
---------, -----------, -----------, ------------, ------------,
153679680 25971865920 25971865920 129859329600 519437318400
238357395880861 238820721143261 86364397717734821 17299975731542641
---------------, ---------------, -----------------, -----------------]
150117385017600 150117385017600 54192375991353600 10838475198270720
Type: List Fraction Integer
Time: 0.03 (OT) = 0.03 sec
Marc Moreno Maza