Merge pull request #1 from mlubin/master

use Calculus for derivative rules
Renamed imag to epsilon