Use dual number for automatic jacobian

[huand]

Specifying jacobian is complex and error prone,

Currently successfully using dual numbers and https://docs.rs/num-dual/0.11.0/num_dual/fn.jacobian.html to automatically get the jacobian from residual computation.

Maybe a feature could be added to this crate to allow residual as dual numbers and the jacobian could be automatically computed?