s13rout1.m

Industrial Mathematics

% Routine 13.1
clear
N= 5; 
% 
for i=1:N
 b(i) = 1/(i+1); 
 for j=1:N 
  A(i,j) = j/(i+j)+  1/(i+j+1);
    end; 
  end;
% 
c=A\b';