vet_1=vector(QQ,(1,2,3))
vet_2=vector(QQ,(-1,2,3))
def vetorial(u,v):
w1=u[1]*v[2]-u[2]*v[1]
w2=u[2]*v[0]-u[0]*v[2]
w3=u[0]*v[1]-u[1]*v[0]
return vector(QQ,(w1,w2,w3))
vetorial(vet_1,vet_2)
P=plot(vet_1,color='black')
P+=plot(vet_2,color='red')
P+= line3d([(-3,0,0),(3,0,0)])
P += line3d([(0,-3,0),(0,3,0)])
P += line3d([(0,0,-3),(0,0,3)])
P+=plot(vetorial(vet_1,vet_2),color='green')
P.show()