function ode(f, initcondition, timespan, alg=Tsit5())
    prob = ODEProblem(f,initcondition,timespan)
    return solve(prob,alg,reltol=1e-8,abstol=1e-8)
end


function plotode(sol)
    plot(sol,title="Solution to ODE")
end

function solve_plot_print(f, initcondition, timespan, alg=Tsit5(), isprint=true)
   sol =  ode(f, initcondition, timespan, alg)
   plotode(sol)
   
end

export ode,plotode, solve_plot_print