Ein Programm für Molecular-Dynamics-Simulation wird entwickelt, um Gleichgewicht der Staubkörner in Plasmen zu studieren. Von der Ausgangsphaseraum-Verteilung der Staubkörner lässt sich das Programm nach Newton-Gesetz der Staubewegung mit einem Modell für Staub-Staub-Interaktions-Potential entwickelt. Zwei Modelle, die abstoßendes sowie attraktives Regime aufweisen, werden benutzt. Die Simulation beginnt mit konstanten-Temperaturbegrenzung, bis ein stationärer Zustand erreicht wird. Sie wird dann ohne der Begrenzung fortgesetzt. Die Fluktuations der Schlüsselparameter werden überwacht. Die Abhängigkeit der Endgleichgewichtzustände von der Temperatur und Dichte in zwei- und dreidimensionalen (3D), endlichen und unendlichen, Systemen werden betrachtet. Gas und flüssig, sowie Clusterstrukturen werden gefunden und geforscht. Eine Monte-Carlo-Simulation des Staubgleichgewichts in externen Feld wird durchgeführt. Transportkoeffizienten für 3D- Staubsysteme werden erhalten.