simplify a bit

Author: ChrisRackauckas

benchmarks/AstroChem/astrochem.jmd

@@ -381,9 +381,7 @@ refsol = solve(oprob, Rodas5P(), abstol=1e-14, reltol=1e-14)
 ```
 
 ```julia
-refsol = solve(oprob, Rodas5P(), abstol=1e-13, reltol=1e-13)
-
-solve(oprob, RadauIIA9(), abstol=1e-9, reltol=1e-9)
+refsol = solve(oprob, Rodas5P(), abstol=1e-13, reltol=1e-13))
 
 # Run Benchmark
 
@@ -392,17 +390,16 @@ reltols = 1.0 ./ 10.0 .^ (9:10)
 
 setups = [
           Dict(:alg=>FBDF()),
-          #Dict(:alg=>QNDF()),
+          Dict(:alg=>QNDF()),
           Dict(:alg=>CVODE_BDF()),
           #Dict(:alg=>ddebdf()),
-          #Dict(:alg=>Rodas4()),
           Dict(:alg=>Rodas5P()),
-          #Dict(:alg=>KenCarp4()),
-          #Dict(:alg=>KenCarp47()),
-          #Dict(:alg=>RadauIIA9()),
+          Dict(:alg=>KenCarp4()),
+          Dict(:alg=>KenCarp47()),
+          Dict(:alg=>RadauIIA9()),
           #Dict(:alg=>rodas()),
           #Dict(:alg=>radau()),
-          #Dict(:alg=>lsoda()),
+          Dict(:alg=>lsoda()),
           #Dict(:alg=>ImplicitEulerExtrapolation(min_order = 5, init_order = 3,threading = OrdinaryDiffEqCore.PolyesterThreads())),
           #Dict(:alg=>ImplicitEulerExtrapolation(min_order = 5, init_order = 3,threading = false)),
           #Dict(:alg=>ImplicitEulerBarycentricExtrapolation(min_order = 5, threading = OrdinaryDiffEqCore.PolyesterThreads())),

benchmarks/AstroChem/nelson.jmd

@@ -174,10 +174,10 @@ sol4 = solve(prob, lsoda(), saveat = 1e10)
 ```julia
 using Plots
 colors = palette(:acton, 5)
-p1 = plot(sol1, vars = (0,11), lc=colors[1], legend = false, titlefontsize = 12, lw = 3, title = "HCO+ solved using Rodas5 (correct solution)")
-p2 = plot(sol2, vars = (0,11), lc=colors[2], legend = false, titlefontsize = 12, lw = 3, title = "HCO+ solved using FBDF")
-p3 = plot(sol3, vars = (0,11), lc=colors[3], legend = false, titlefontsize = 12, lw = 3, title = "HCO+ solved using lsoda")
-p4 = plot(sol4, vars = (0,11), lc=colors[4], legend = false, titlefontsize = 12, lw = 3, title = "HCO+ solved using lsoda with saveat")
+p1 = plot(sol1, vars = (0,11), lc=colors[1], legend = false, titlefontsize = 12, lw = 3, title = "Rodas5")
+p2 = plot(sol2, vars = (0,11), lc=colors[2], legend = false, titlefontsize = 12, lw = 3, title = "FBDF")
+p3 = plot(sol3, vars = (0,11), lc=colors[3], legend = false, titlefontsize = 12, lw = 3, title = "lsoda")
+p4 = plot(sol4, vars = (0,11), lc=colors[4], legend = false, titlefontsize = 12, lw = 3, title = "lsoda with saveat")
 ```
 
 ```julia