crossprd
==============================================
Purpose
----------------
Computes the cross-products (vector products) of sets of 3x1 vectors.
Format
----------------
.. function:: z = crossprd(x, y)
:param x: each column is treated as a 3x1 vector.
:type x: 3xK matrix
:param y: each column is treated as a 3x1 vector.
:type y: 3xK matrix
:return z: each column is the cross-product
(sometimes called vector product) of the
corresponding columns of *x* and *y*.
:rtype z: 3xK matrix
Examples
----------------
::
// First matrix
x = { 10 4,
11 13,
14 13 };
// Second matrix
y = { 3 11,
5 12,
7 9 };
z = crossprd(x, y);
After the above code,
::
7 -39
z = -28 107
17 -95
Remarks
-------
The cross-product vector *z* is orthogonal to both *x* and *y*. :code:`sumc(x .* z)`
and :code:`sumc(y .* z)` will be Kx1 vectors, all of whose elements are 0
(except for rounding error).
Source
------------
crossprd.src