This paper introduces a new non-stationary calculation of characteristic earthquake probability, motivated and derived from the Mohr-Coulomb failure criterion. A total of six parameters are present in the model, including earthquake return period, fault-plane strength parameters, rock unit weight, focal depth, and the level of uncertainty of annual stress increment. In addition to the detailed derivations in the algorithm, a model demonstration is also presented in the paper.