HP 39g-Grafenberechner - Handbuch (Seite 158)

Exakte Berechnungen und Mathematik mit HP40G
158 Arithmetische Programme
sonst
D+1 -
> D
/wenn/
/solange/
Ergebnis FACT
/Funktion/
Erste Verbesserung
Man testet nur den Teiler D zwischen 2 q E(
ÖN).
Wenn N = D1* D2 so hat man also:
D1
£ E(ÖN) oder D2 £ E(ÖN), sonst würde:
D1 * D2
³ (E(ÖN) + 1)² > N.
Funktion facprem(N)
lokal D FACT
2-> D
{} -> FACT
solange D*D
£ N durchführen
wenn N mod D = 0 so
FACT + D -> FACT
N/D -> N
sonst
D+1-> D
/wenn/
/solange/
FACT + N -> FACT
Ergebnis FACT (Ergebnis)
/Funktion/
Zweite Verbesserung