function v = helix(x,y,z)
        v = -y*e1 + x*e2 + (x^2+y^2)*z*e3;                                                
