Mathematical models of neurons are formulated as nonlinear stochastic ordinary or partial differential equations. Progress in determining the activity of single neurons for models of varying degrees of mathematical complexity by analytical and simulation methods will be described. Exact results for the global activity in a network of such elements will also be obtained and in particular their approximations by means of diffusion processes. Analogous models for viral dynamics and epidemic networks may also be considered.