Update scipy.md

Author: ChrisRackauckas

docs/src/optimization_packages/scipy.md

@@ -97,7 +97,7 @@ sol = solve(prob, ScipyCOBYLA())
 ### Differential evolution (global) with custom options
 
 ```@example SciPy3
-using Optimization, OptimizationSciPy, Random
+using Optimization, OptimizationSciPy, Random, Statistics
 Random.seed!(123)
 
 ackley(x, p) = -20exp(-0.2*sqrt(mean(x .^ 2))) - exp(mean(cos.(2π .* x))) + 20 + ℯ
@@ -129,4 +129,4 @@ If SciPy raises an error it is re-thrown as a Julia `ErrorException` carrying th
 
 ## Contributing
 
-Bug reports and feature requests are welcome in the [Optimization.jl](https://github.com/SciML/Optimization.jl) issue tracker.  Pull requests that improve either the Julia wrapper or the documentation are highly appreciated. 
+Bug reports and feature requests are welcome in the [Optimization.jl](https://github.com/SciML/Optimization.jl) issue tracker.  Pull requests that improve either the Julia wrapper or the documentation are highly appreciated.