function main
  mrow = 5;
  ncol = 5;
  A = lo_matrix ( );
  disp ( A(1:10,1:10) );
  x = zeros ( mrow * ncol, 1 );
  x(7) = 1;
  Ax = mod ( A * x, 2 );
  disp ( Ax );
  y = round ( A \ Ax )
  return
end




