28 lines
789 B
Plaintext
28 lines
789 B
Plaintext
DIM U0%(2,2), U1%(2,2), U2%(2,2), seed%(2)
|
|
U0%() = 1, -2, 2, 2, -1, 2, 2, -2, 3
|
|
U1%() = 1, 2, 2, 2, 1, 2, 2, 2, 3
|
|
U2%() = -1, 2, 2, -2, 1, 2, -2, 2, 3
|
|
|
|
seed%() = 3, 4, 5
|
|
FOR power% = 1 TO 7
|
|
all% = 0 : prim% = 0
|
|
PROCtri(seed%(), 10^power%, all%, prim%)
|
|
PRINT "Up to 10^"; power%, ": " all% " triples" prim% " primitives"
|
|
NEXT
|
|
END
|
|
|
|
DEF PROCtri(i%(), mp%, RETURN all%, RETURN prim%)
|
|
LOCAL t%() : DIM t%(2)
|
|
|
|
IF SUM(i%()) > mp% ENDPROC
|
|
prim% += 1
|
|
all% += mp% DIV SUM(i%())
|
|
|
|
t%() = U0%() . i%()
|
|
PROCtri(t%(), mp%, all%, prim%)
|
|
t%() = U1%() . i%()
|
|
PROCtri(t%(), mp%, all%, prim%)
|
|
t%() = U2%() . i%()
|
|
PROCtri(t%(), mp%, all%, prim%)
|
|
ENDPROC
|