pub fn merge_dual(y: Vec<f64>, dy: Vec<f64>) -> Vec<Dual>