Compute with Properties of Isotopes

Mathematica 6 has a full database of nuclear isotopes. This computes a decay network for isotopes with A=100 -102.

With[{is = 

   Select[IsotopeData[], 

    100 <= IsotopeData[#, "MassNumber"] <= 102 &]}, 

 LayeredGraphPlot[

  Table[(IsotopeData[#, "Symbol"] & /@ (i -> #)) & /@ 

     IsotopeData[i, "DaughterNuclides"], {i, is}] // Flatten, 

  DirectedEdges -> True, VertexLabeling -> True]]

Out[1]=