Uses a Newton Raphson method where the Jacobian is computed via hyperdual numbers.
pub use crate::md::Variable;
pub use crate::md::Vary;