include("odehelper.jl")

function lorenz!(du,u,p,t)
    du[1] = 10.0*(u[2]-u[1])
    du[2] = u[1]*(28.0-u[3]) - u[2]
    du[3] = u[1]*u[2] - (8/3)*u[3]
end
solve_plot_print(lorenz!,[1.0;0.0;0.0],(0.0,100.0))
#http://www.math.utah.edu/~gustafso/2250systems-de.pdf
function brinetankvascade!(du,u,p,t)
    du[1] = -(1/2)*u[1]
    du[2] = (1/2) * u[1] - (1/4)*u[2]
    du[3] = (1/4)*u[2] - (1/6)*u[3]
end

solve_plot_print(lorenz!,[20,40,60],(0.0,10.0))


