function a = quaternion_1_matrix ( )

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

  return
end
