A dust grain in a plasma acquires an electric charge by collecting electron and ion currents. These currents consist of discrete charges, causing the charge to fluctuate around an equilibrium value /spl lang/Q/spl rang/. Electrons and ions are collected at random intervals and in a random sequence, with probabilities that depend on the grain's potential. We developed a model for these probabilities and implemented it in a numerical simulation of the collection of individual ions and electrons, yielding a time series Q(t) for the grain's charge. Electron emission from the grain is not included, although it could be added easily to our method. We obtained the power spectrum and the RMS fluctuation level, as well as the distribution function of the charge. Most of the power in the spectrum lies at frequencies much lower than 1//spl tau/, the inverse charging time. The RMS fractional fluctuation level varies as 0.5 /spl verbar//spl lang/N/spl rang//spl verbar//sup /spl minus/1/2/, where /spl lang/N/spl rang/=/spl lang/Q/spl rang//e is the average number of electron charges on the grain. This inverse square-root scaling means that fluctuations are most important for small grains. We also show that very small grains can experience fluctuations to neutral and positive polarities, even in the absence of electron emission.< <ETX xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">></ETX>