在Julia中，在同一个图形上绘制多个ODE解决方案的最佳方法是什么？

我有一个常微分方程，想求解它，并绘制一些初始条件的解。开始时，我决定做一个“for”来增加初始条件的值，然后一个接一个地绘制解决方案

using DifferentialEquations, Plots

for i = -2.00:0.25:2.00
    f(x,p,t) = x-x^3
    x0 = i
    tspan = (0.00,2.00)
    prob = ODEProblem(f,x0,tspan)
    sol = solve(prob)
    if i==-2.00
        plot(sol,linewidth=4, title="Solution",
             xaxis="t",yaxis="x(t)",legend=false)
    else
        plot!(sol,linewidth=4, title="Solution",
            xaxis="t",yaxis="x(t)",legend=false)
    end
end
savefig("Graphic.png")

---

有更好的方法吗？

您可以创建一个只有标题、轴标签等的空绘图，然后用解决方案填充绘图：

# main plot settings
plot(linewidth=4, title="Solution",
     xaxis="t",yaxis="x(t)",legend=false)

for i = -2.00:0.25:2.00
    f(x,p,t) = x-x^3
    x0 = i
    tspan = (0.00,2.00)
    prob = ODEProblem(f,x0,tspan)
    sol = solve(prob)

    # add to existing plot
    plot!(sol)
end

---

好主意！非常感谢你！不客气！如果答案解决了你原来的问题，请别忘了标记。这样人们将来就可以过滤搜索了。