function v = radial(x,y,z)
     v = x*e1+y*e2+z*e3;
