syms x y z;
F = [x*y^2, x*y, z];
vectline(F,[x,y,z],[-1,1,-2,3,1,3])