function a = quaternion_j_matrix ( )

%*****************************************************************************80
%
%% quaternion_j_matrix returns a 4 by 4 matrix that behaves like the quaternion unit J.
%
%  Formula:
%
%    0  0  1  0
%    0  0  0  1
%   -1  0  0  0
%    0 -1  0  0
%
%  Properties:
%
%    J * 1 =   J
%    J * I = - K
%    J * J = - 1
%    J * K =   I
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    19 October 2007
%
%  Author:
%
%    John Burkardt
%
%  Output:
%
%    real A(4,4), the matrix.
%
  a = [   0.0,   0.0,   1.0,   0.0; ...
          0.0,   0.0,   0.0,   1.0; ...
         -1.0,   0.0,   0.0,   0.0; ...
          0.0,  -1.0,   0.0,   0.0 ];

  return
end
