Non standard finite difference scheme is proposed. The method is applied to solve fractional point kinetics equations. Comparisons are made with the classical solutions. Comparisons are made for different values of the anomalous order α. Stability conditions are investigated to control the step size.