change coupled sde example to use trajectory

Author: oameye

docs/src/CoupledSDEs.md

@@ -35,7 +35,8 @@ Lyapunov spectrum of a `CoupledSDEs` in the absence of noise, here exemplified b
 FitzHugh-Nagumo model, can be computed by typing:
 
 ```julia
-using DynamicalSystems
+using DynamicalSystemsBase, StaticArrays
+using DynamicalSystemsBase: idfunc
 
 function fitzhugh_nagumo(u, p, t)
     x, y = u
@@ -47,8 +48,8 @@ function fitzhugh_nagumo(u, p, t)
     return SA[dx, dy]
 end
 
-sys = CoupledSDEs(fitzhugh_nagumo, idfunc, zeros(2), [1.,3.,1.,1.,1.,0.], 0.1)
-ls = lyapunovspectrum(CoupledODEs(sys), 10000)
+sde = CoupledSDEs(fitzhugh_nagumo, idfunc, zeros(2), [1.,3.,1.,1.,1.,0.], 0.1)
+tr = trajectory(sde, 10)
 ```
 
 ## Available noise processes

src/core_systems/continuous_time_sde.jl

@@ -183,7 +183,7 @@ end
 
 Converts a [`CoupledSDEs`](@ref) into [`CoupledODEs`](@ref).
 """
-function CoupledODEs(sys::CoupledSDEs; diffeq=DEFAULT_STOCH_DIFFEQ, t0=0.0)
+function CoupledODEs(sys::CoupledSDEs; diffeq=DEFAULT_DIFFEQ, t0=0.0)
     return CoupledODEs(
         sys.integ.f, SVector{length(sys.integ.u)}(sys.integ.u), sys.p0; diffeq=diffeq, t0=t0
     )
@@ -262,7 +262,7 @@ Identity function for a diffusion function `g` of `CoupledSDEs` (out-of-place).
 Equivalent to `(u, p, t) -> ones(length(u))`,
 """
 function idfunc(u, p, t)
-    return typeof(u)(ones(eltype(u), length(u)))
+    return  SVector{length(u)}(ones(eltype(u), length(u)))
 end;
 
 """