Improve code indentation
Remove unused differetiate() function
Integrate with tau, not eps; use a less complicated integration formula
Add caching support, improves speed
Better precision: eps=1e-7, tau=1e-4