Hallo Peter,das folgt der üblichen Mathematicaprogrammierung, machen wir zunächst einen Arrayplot der bei Wolfram bekannten Isotope
In[45]:= Clear[isotopenPlot]
isotopenPlot[] := With[{zMax = 120, aMax = 300},
ArrayPlot[
Table[If[
Head[Quiet[
Check[IsotopeData[{o, oo}, "StandardName"],
0, {IsotopeData::notent}]]] === String,
1, (* else *)
0 ], {o, zMax}, {oo, aMax}],
ColorRules -> {1 -> Cyan, 0 -> Yellow}, Mesh -> {zMax, aMax}/5,
Frame -> True, FrameTicks -> Automatic, FrameLabel -> {"Z", "A"},
PlotLabel -> "Wolfram::Isotope"]
]
auf einschlägigen Seiten findet man, dass Kernladungszahlen > 120 und
Massenzahlen > 300 nicht gemischelt werden konnten, demnach nicht
existieren. Das gibt das Bild im Anhang.
Jetzt müssten Sie wg.
X^{A}_{Z} -> Y^{A-4}_{Z-2} + He^{4}_{2} + dE
was der Alphazerfall in LaTeX-Notation ist, die Kette herunterlaufen, mit
DaughterNuclides kontrollieren, dass Y ein daughter nuclide von X ist, die
Halbwertszeit abholen und dieses rekursiv ausführen, bis kein Zweig mehr
Töchter meldet. Am Schluss einen schönen Graphen fertigen.
Go ahead! Grüsse Udo.On Sun, 08 Sep 2019 20:30:31 +0200, Peter Klamser via demug <demug@XXXXXXX.ch> wrote:
Hallo,Ich würde gerne mit MMA12 über IsotopeData["Pu241", "HalfLife"] eine Liste der Zerfallsprodukte in der Zerfallskette mit der jeweiligen Halbwertszeitnur für Alphazerfälle erstellen. Die Beispiele bei wri.com sind sehr komplex. Weiß jemand wie das geht? Danke sagt Peter _______________________________________________ DMUG Deutschsprachiges Mathematica-Forum demug@XXXXXXX.ch http://www.mathematica.ch/mailman/listinfo/demug Archiv: http://www.mathematica.ch/archiv.html
klamser-Isotope.jpg
Description: JPEG image
_______________________________________________ DMUG Deutschsprachiges Mathematica-Forum demug@XXXXXXX.ch http://www.mathematica.ch/mailman/listinfo/demug Archiv: http://www.mathematica.ch/archiv.html
