  Efficient simulation of probabilistic memristors and their networks requires
novel modeling approaches. One major departure from the conventional memristor
modeling is based on a master equation for the occupation probabilities of
network states [arXiv:2003.11011 (2020)]. In the present article, we show how
to implement such master equations in SPICE - a general-purpose circuit
simulation program. In the case studies, we simulate the dynamics of ac-driven
probabilistic binary and multi-state memristors, and dc-driven networks of
probabilistic binary and multi-state memristors. Our SPICE results are in
perfect agreement with known analytical solutions. Examples of LTspice codes
are included.
