% Estimation of the condition number for the matrix obtained
% for non-orthogonal polynomial approximation
function spectr_01
echo off, clc
clear all
Nmin = 2;
Nmax = 500;
N = Nmin;
k = 1;
while ( N <= Nmax );
% asseembling the matrix
for i=1:N
for j=1:N
A(i,j) = 1.0 / (i + j + 1);
end;
end;
% estimating the condition number;
D(k,1) = N;
D(k,2) = cond(A);
k = k + 1;
N = N * 2;
end;
figure(1); loglog(D(:,1),D(:,2),'-x');