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

mag = sqrt(x1 * x1 + y1 * y1 + z1 * z1)
sqrt(dot([x1 y1 z1], [x1 y1 z1]))
norm([x1 y1 z1])
