%!	View coefficients of row indexed sparce matrix, A, stored as ija(k) & sa(k)
%!	ija(1)=N+2
%!	ija(N+1)= one greater than the index in sa of the last off-diagonal element 
%!		of last row
%!	ija(k), k=1,2,...,N contain the index of sa that contain the 
%!		first off-diagonal element of the corresponding row of the matrix
%!	ija(k), k=N+2,N+3,.... contains the column number of the 
%!		corresponding element in sa
%!	sa(k), k=1,2,...,N stores the diagonal elements of A
%!	sa(N+1)is not used
%!	sa(k), k=N+2,N=3,... contains the off-diagonal components of A
clear all
clf reset
load('matrix.dat');
ija=matrix(:,1);
sa=matrix(:,2);
N=ija(1)-2
A=zeros(N);
for k=1:N
   A(k,k)=sa(k);	% diagonal elements
end
for k=1:N			% for each row of matrix A
   for kk=ija(k):(ija(k+1)-1)
      A(k,ija(kk))=sa(kk);
   end
end
spy(A)
   