Adds a threshold controller.

Introduces a simple thresholding controller, which provides full power
if the temperature is below the threshold and off when it's above.
Adds odebdf function for stiff systems.

Adds a function which uses the backward differential formula as a
solver. Odeint ran into trouble with stiff systems such as the threshold
controller. The odebdf function allows simulation of such systems
without problems.