In this work, we solve nonlinear systems of ordinary differential equations coupled to noisy forcing, commonly used for models of neurons such as the Hodgkin-Huxley equation, over a memristor crossbar-based computing system. We demonstrate stability and faithfulness of the distributions even under the effects of nonidealities of the memristors and the system itself. We investigate the properties of the dynamical systems under quantization faithfulness, varying the level of precision of the fixed point integer representation and concluding that 24 bits is enough for the solution of the Hodgkin-Huxley equations, demon- strating that our solver can operate with both high precision and achieve speedups with low-precision approximate computation.