workaround v0.5.2 bug

Author: ChrisRackauckas

src/solutions/solution_interface.jl

@@ -41,7 +41,8 @@ function done(sol::AbstractTimeseriesSolution,state)
   state >= length(sol)
 end
 
-const DEFAULT_PLOT_FUNC = (x...) -> (x...)
+DEFAULT_PLOT_FUNC(x...) = (x...)
+DEFAULT_PLOT_FUNC(x,y,z) = (x,y,z) # For v0.5.2 bug
 
 @recipe function f(sol::AbstractTimeseriesSolution;
                    plot_analytic=false,
@@ -309,8 +310,24 @@ function solplot_vecs_and_labels(dims,vars,plot_timeseries,plott,sol,plot_analyt
     end
 
     f = x[1]
-    tmp = f.(tmp...)
 
+    if VERSION < v"0.6-" && length(tmp) == 3
+      # Work around bug
+      #=
+      f = (x...) -> x
+      test = Any[rand(5),rand(5)]
+      f.(test...)
+      test = Any[rand(5),rand(5),rand(5)]
+      f.(test...)
+      =#
+      res = Any[]
+      for i in eachindex(tmp[1])
+        push!(res,(tmp[1][i],tmp[2][i],tmp[3][i]))
+      end
+      tmp = res
+    else
+      tmp = f.(tmp...)
+    end
 
     tmp = tuple((getindex.(tmp,i) for i in eachindex(tmp[1]))...)
     for i in eachindex(tmp)