Merge pull request #1251 from SciML/ChrisRackauckas-patch-3

Fix docstring

Author: ChrisRackauckas

src/systems/diffeqs/odesystem.jl

@@ -19,7 +19,7 @@ eqs = [D(x) ~ σ*(y-x),
        D(y) ~ x*(ρ-z)-y,
        D(z) ~ x*y - β*z]
 
-de = ODESystem(eqs,t,[x,y,z],[σ,ρ,β])
+@named de = ODESystem(eqs,t,[x,y,z],[σ,ρ,β])
 ```
 """
 struct ODESystem <: AbstractODESystem

src/systems/diffeqs/sdesystem.jl

@@ -23,7 +23,7 @@ noiseeqs = [0.1*x,
             0.1*y,
             0.1*z]
 
-de = SDESystem(eqs,noiseeqs,t,[x,y,z],[σ,ρ,β])
+@named de = SDESystem(eqs,noiseeqs,t,[x,y,z],[σ,ρ,β])
 ```
 """
 struct SDESystem <: AbstractODESystem

src/systems/discrete_system/discrete_system.jl

@@ -18,7 +18,7 @@ eqs = [next_x ~ σ*(y-x),
        next_y ~ x*(ρ-z)-y,
        next_z ~ x*y - β*z]
 
-de = DiscreteSystem(eqs,t,[x,y,z],[σ,ρ,β])
+@named de = DiscreteSystem(eqs,t,[x,y,z],[σ,ρ,β])
 ```
 """
 struct DiscreteSystem <: AbstractTimeDependentSystem

src/systems/jumps/jumpsystem.jl

@@ -22,7 +22,7 @@ affect₂ = [I ~ I - 1, R ~ R + 1]
 j₁      = ConstantRateJump(rate₁,affect₁)
 j₂      = ConstantRateJump(rate₂,affect₂)
 j₃      = MassActionJump(2*β+γ, [R => 1], [S => 1, R => -1])
-js      = JumpSystem([j₁,j₂,j₃], t, [S,I,R], [β,γ])
+@named js      = JumpSystem([j₁,j₂,j₃], t, [S,I,R], [β,γ])
 ```
 """
 struct JumpSystem{U <: ArrayPartition} <: AbstractTimeDependentSystem

src/systems/nonlinear/nonlinearsystem.jl

@@ -15,7 +15,7 @@ $(FIELDS)
 eqs = [0 ~ σ*(y-x),
        0 ~ x*(ρ-z)-y,
        0 ~ x*y - β*z]
-ns = NonlinearSystem(eqs, [x,y,z],[σ,ρ,β])
+@named ns = NonlinearSystem(eqs, [x,y,z],[σ,ρ,β])
 ```
 """
 struct NonlinearSystem <: AbstractTimeIndependentSystem

src/systems/optimization/optimizationsystem.jl

@@ -13,7 +13,7 @@ $(FIELDS)
 @parameters σ ρ β
 
 op = σ*(y-x) + x*(ρ-z)-y + x*y - β*z
-os = OptimizationSystem(eqs, [x,y,z],[σ,ρ,β])
+@named os = OptimizationSystem(eqs, [x,y,z],[σ,ρ,β])
 ```
 """
 struct OptimizationSystem <: AbstractTimeIndependentSystem

src/systems/pde/pdesystem.jl

@@ -31,7 +31,7 @@ bcs = [u(t,0) ~ 0.,# for all t > 0
 domains = [t ∈ (0.0,1.0),
            x ∈ (0.0,1.0)]
 
-pde_system = PDESystem(eq,bcs,domains,[t,x],[u])
+@named pde_system = PDESystem(eq,bcs,domains,[t,x],[u])
 ```
 """
 struct PDESystem <: ModelingToolkit.AbstractMultivariateSystem