var('p')
w=vector([0,0,p])
print w
var('x')
var('y')
var('z')
r=vector([x , y , z])
print r
w.cross_product(r)

/* code_ends