add other README example

Author: ChrisRackauckas

README.md

@@ -103,3 +103,5 @@ sol = solve(prob,Rodas5())
 
 using Plots; plot(sol,vars=(:α,Symbol(lorenz1.x),Symbol(lorenz2.y)))
 ```
+
+![](https://user-images.githubusercontent.com/1814174/79122361-6fdaca80-7d65-11ea-87fd-0f6c4a85cd0d.png)

test/lowering_solving.jl

@@ -24,3 +24,41 @@ tspan = (0.0,100.0)
 prob = ODEProblem(sys,u0,tspan,p,jac=true)
 sol = solve(prob,Tsit5())
 #using Plots; plot(sol,vars=(:x,:y))
+
+@parameters t σ ρ β
+@variables x(t) y(t) z(t)
+@derivatives D'~t
+
+eqs = [D(x) ~ σ*(y-x),
+       D(y) ~ x*(ρ-z)-y,
+       D(z) ~ x*y - β*z]
+
+lorenz1 = ODESystem(eqs,name=:lorenz1)
+lorenz2 = ODESystem(eqs,name=:lorenz2)
+
+@variables α
+@parameters γ
+connections = [0 ~ lorenz1.x + lorenz2.y + α*γ]
+connected = ODESystem(connections,t,[α],[γ],systems=[lorenz1,lorenz2])
+
+u0 = [lorenz1.x => 1.0,
+      lorenz1.y => 0.0,
+      lorenz1.z => 0.0,
+	  lorenz2.x => 0.0,
+	  lorenz2.y => 1.0,
+	  lorenz2.z => 0.0,
+	  α => 2.0]
+
+p  = [lorenz1.σ => 10.0,
+      lorenz1.ρ => 28.0,
+      lorenz1.β => 8/3,
+	  lorenz2.σ => 10.0,
+	  lorenz2.ρ => 28.0,
+	  lorenz2.β => 8/3,
+	  γ => 2.0]
+
+tspan = (0.0,100.0)
+prob = ODEProblem(connected,u0,tspan,p)
+sol = solve(prob,Rodas5())
+
+#using Plots; plot(sol,vars=(:α,Symbol(lorenz1.x),Symbol(lorenz2.y)))