function B = ssmat2(nl,pk)
% Function to generate Simon and Sturm matrices.
%
% by D. Holmgren, Mar 13 94. holmgren@phobos.astro.uwo.ca
% This version for use with SSNORMAL.
B = zeros(nl,nl);
for i = 1:nl
   for j = 1:nl
    if i == j + pk
    B(i,j) = 1.;
    end
   end
end