% Routine 13.2
N=5;
for i=1:N 
 for j=1:N
  G(i,j) = 1/(i+j+1)-2/(i+j+2) + 1/(i+j+3); 
  S(i,j) = i*j/(i+j-1)-(2*i*j+i+j)/(i+j)+(i+1)*(j+1)/(i+j+1);  
  end; 
  end; 
%
[X,D] = eig(S,-G);