Auch grosse Primzahlen koennen leicht mit Mathematica aufgestellt
werden; dabei wird aber vorausgesetzt, dass QPrime[] 100 % verlaesslich ist
(was ich nicht ueberpruefen kann).
Man startet von einer Zufallszahl der gewuenschten Groesse
(npp = Zahl der Stellen), incremtiert um 1, solange bis QPrime[]
True gibt.
Auf meinem PowerPC mit 300 MHz dauert das einige Sekunden.
In[1]:=
liprime[npp_] := Module[{i, pn },
For[i = Random[Integer, 10^npp {1, 9}], PrimeQ[i] == False, pn = i + 1;
i++]; Return[pn]]
In[2]:=
liprime[150] // Timing
**************************************************************************
ao. Prof. Dr. B. Schnizer Institut fuer Theoretische Physik
Tel: +43 (316) 873-8173, -8171 (Secr.) Technische Universitaet Graz
Fax: +43 (316) 873 8678 Petersgasse 16
e-mail: schnizer@XXXXXXX.at A-8010 Graz
http://www.itp.tu-graz.ac.at/schnizer Austria/Europe
https://online.tu-graz.ac.at/tug_online/tgo_org.display?CORG=2382
***************************************************************************