week13.sagews

x,y,z,t = var('x y z t')

F(x,y,z) = (3*y+4*z,6*z,0)
P = plot_vector_field3d(F(x,y,z),(x,-1,3),(y,-1,3),(z,-1,3),plot_points=7)
Q = parametric_plot3d([1+7*t+9*t^2,1+6*t,1],(t,-1/3,2/9),color="black")
(P+Q).show()

3D rendering not yet implemented

G(x,y) = (x+y,y)
P = plot_vector_field(G(x,y),(x,-1,1),(y,-1,1),aspect_ratio=1.0)
Q = parametric_plot([-.25*exp(t)+.25*exp(t)*t,.25*exp(t)],(t,-5,1.45))
(P+Q).show()