I refer here to the Customization Manual where DEFINE_SOURCE is explained.
Another good reference is to go through the basics with reading the simple introduction of Patankar "Numerical Heat Transfer and Fluid Flow".
Linearization will help on stability.