test(docs): update Multinomial doc tests

src/distribution/multinomial.rs

@@ -19,10 +19,10 @@ use rand::Rng;
 /// ```
 /// use statrs::distribution::Multinomial;
 /// use statrs::statistics::MeanN;
-/// use nalgebra::DVector;
+/// use nalgebra::vector;
 ///
-/// let n = Multinomial::new(&[0.3, 0.7], 5).unwrap();
-/// assert_eq!(n.mean().unwrap(), DVector::from_vec(vec![1.5, 3.5]));
+/// let n = Multinomial::new_from_nalgebra(vector![0.3, 0.7], 5).unwrap();
+/// assert_eq!(n.mean().unwrap(), (vector![1.5, 3.5]));
 /// ```
 #[derive(Debug, Clone, PartialEq)]
 pub struct Multinomial<D>
@@ -87,9 +87,10 @@ where
     ///
     /// ```
     /// use statrs::distribution::Multinomial;
+    /// use nalgebra::dvector;
     ///
-    /// let n = Multinomial::new(&[0.0, 1.0, 2.0], 3).unwrap();
-    /// assert_eq!(n.p(), [0.0, 1.0, 2.0]);
+    /// let n = Multinomial::new(vec![0.0, 1.0, 2.0], 3).unwrap();
+    /// assert_eq!(*n.p(), dvector![0.0, 1.0/3.0, 2.0/3.0]);
     /// ```
     pub fn p(&self) -> &OVector<f64, D> {
         &self.p
@@ -103,7 +104,7 @@ where
     /// ```
     /// use statrs::distribution::Multinomial;
     ///
-    /// let n = Multinomial::new(&[0.0, 1.0, 2.0], 3).unwrap();
+    /// let n = Multinomial::new(vec![0.0, 1.0, 2.0], 3).unwrap();
     /// assert_eq!(n.n(), 3);
     /// ```
     pub fn n(&self) -> u64 {
@@ -300,6 +301,7 @@ mod tests {
         distribution::{Discrete, DiscreteCDF, Multinomial},
         statistics::{Max, MeanN, Min, Mode, VarianceN},
     };
+    use approx::UlpsEq;
     use nalgebra::{
         dmatrix, dvector, matrix, vector, Const, DimMin, Dyn, Matrix, OMatrix, OVector,
         VecStorage,