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