reduction

We should be able to attach names to rules, so that normalization steps (and error messages?) can be better understood by users. For example:

normalize(@term(diff(x + y, x)))
 - @term(diff(x,x) + diff(y,x)) by sum rule in differentiation
 - @term(one(x) + diff(y, x)) by linear rule of differentiation
 - @term(1 + diff(y, x)) by multiplicative identity of a number
 - @term(1 + zero(x)) 

Attempting to generate the sphinx documentation fails.

The may be a configuration issue and/or a code documentation style issue.

We need to pick a format of code documatation that works (preferably astropy, or maybe numpy), convert all doc strings to it, and get the documentation working.