function [ i, j, v ] = usa_st ( )

%*****************************************************************************80
%
%% usa_st() returns the USA matrix in sparse triplet (ST) format.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    05 November 2022
%
%  Author:
%
%    John Burkardt
%
%  Output:
%
%    integer i(210), j(210): the indices of nonzero entries. 
%
%    real v(210): the nonzero values.
%
  i = [ ...
    1, 1, 1, 1, 3, 3, 3, 3, 4, 4, ...
    4, 4, 4, 4, 5, 5, 5, 6, 6, 6, ...
    6, 6, 6, 7, 7, 7, 8, 8, 8, 9, ...
    9,10,10,10,10,10,12,12,12,12, ...
   12,12,13,13,13,13,13,14,14,14, ...
   14,15,15,15,15,15,15,16,16,16, ...
   16,17,17,17,17,17,17,17,18,18, ...
   18,19,20,20,20,20,21,21,21,21, ...
   21,22,22,22,23,23,23,23,24,24, ...
   24,24,25,25,25,25,25,25,25,25, ...
   26,26,26,26,27,27,27,27,27,27, ...
   28,28,28,28,28,29,29,29,30,30, ...
   30,31,31,31,31,32,32,32,32,32, ...
   33,33,33,33,34,34,34,35,35,35, ...
   35,35,36,36,36,36,36,36,37,37, ...
   37,37,38,38,38,38,38,38,39,39, ...
   40,40,41,41,41,41,41,41,42,42, ...
   42,42,42,42,42,42,43,43,43,43, ...
   44,44,44,44,44,45,45,45,46,46, ...
   46,46,46,47,47,48,48,48,48,48, ...
   49,49,49,49,50,50,50,50,50,50 ];

  j = [ ...
    9,10,42,24,31,44,28, 5,18,24, ...
   42,25,36,43, 3,28,37,31,36,16, ...
   27,50,44,39,21,32,30,38,20,10, ...
    1,40,33,42, 1, 9,47,37,28,44, ...
   50,26,49,15,25,17,14,13,17,35, ...
   22,23,41,27,25,13,49,36,25,27, ...
    6,42,46,48,35,14,13,25,24, 4, ...
   43,29, 8,38,48,46,29,45,32, 7, ...
   39,49,14,35,34,41,15,49, 1,42, ...
    4,18, 4,42,17,13,15,27,16,36, ...
   12,50,41,34,16,25,15,41,50, 6, ...
    3,44,12,37, 5,45,21,19,32,38, ...
    8,43,36, 6, 3,38,30, 7,21,45, ...
   46,42,10,40,26,41,23,22,14,17, ...
   48,38,43, 4,25,16, 6,31, 5,28, ...
   12,47,35,48,20, 8,30,32,21, 7, ...
   33,10,27,15,23,34,26,50, 1,10, ...
   33,46,17,25, 4,24,18, 4,36,31, ...
    3, 6,50,12,28,32,21,29,20,48, ...
   17,42,33,37,12,46,20,38,35,17, ...
   23,15,13,22, 6,27,41,26,12,44 ];

  v = ones ( 1, 210 );

  return
end
 
