F=(1,x). On va ici dessiner le champ, et la ligne de courant demandée.
var('x,y,t')P(x,y)=1Q(x,y)=xChamp=plot_vector_field([P(x,y),Q(x,y)],(x,-3,3),(y,-3,3),color="blue")Courbe=plot(x^2/2,(x,-3,3),color="red")# La solutuion, calculée manuellement.show(Champ+Courbe,ymin=-3,ymax=3,aspect_ratio=1,figsize=7)