How to use PyPlot.jl to draw the graph of DifferentialEquations

[Minghao_Du]

The sol of the “DifferentialEquations.jl” package can be nicely drawn with Plots, like the following

But if you use the “PyPlot.jl” package, how do you draw the sol curve?

The IDE I use is Pluto

using DifferentialEquations
import PyPlot as plt

let
	l = 1.0                            
	m = 1.0                           
	g = 9.81                           

	function pendulum!(du,u,p,t)
   		du[1] = u[2]                 
    	du[2] = -3g/(2l)*sin(u[1]) + 3/(m*l^2)*p(t) 
	end

	θ₀ = 0.01                          
	ω₀ = 0.0                          
	u₀ = [θ₀, ω₀]                     
	tspan = (0.0,10.0)                

	M = t->0.1sin(t)                  

	prob = ODEProblem(pendulum!,u₀,tspan,M)
	sol = solve(prob)
	
	plt.clf()
	plt.plot(sol.t, sol[1,:])
	plt.gcf()
end

[Minghao_Du]

This is the image I drew with PyPlot.jl and it is not a curve

[tomaklutfu]

sol.t has adaptive time steps. The solution object is callable so you can use a range for time steps and call solution object with that range. This explains manual plotting.

[Minghao_Du]

Thanks!