t = LinRange(0,2,10) # define the times to solve at
y0 = 1
y = ode( frhs, y0, t) # solve the problem (find y(t))
soln = solve(ODEProblem((y,p,t) -> frhs(y,t), y0, (t[1],t[end]), ["dummy"]))
#-------------- plot the results
plt.resetfontsizes(); plt.scalefontsizes(1.5)
plt.scatter(t,y, label="Euler", markershape=:square, ms=6)
plt.scatter!(sol.t, sol.u, label="Julia", color="black")
plt.plot!(t->y0*exp(-t), label="Exact", color="black")
plt.plot!(xlabel="time", ylabel="y")
plt.plot!(foreground_color_legend=nothing, background_color_legend=nothing)