x1 = 1;
y1 = 2;
z1 = -3;

mag = sqrt(x1 * x1 + y1 * y1 + z1 * z1);

ux = x1 / mag
uy = y1 / mag
uz = z1 / mag
