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

x2 = -4;
y2 = 5;
z2 = 6;

dprod = x1 * x2 + y1 * y2 + z1 * z2
dot([x1 y1 z1], [x2 y2 z2])
