Saso@1: function t_mat = poly_mat(tt, dgr, N)
Saso@1: t_mat = ones(dgr,N);
Saso@1: if dgr > 1
Saso@1:   t_mat(2,:) = tt;
Saso@1:   for k=3:dgr
Saso@1:     t_mat(k,:) = t_mat(k-1,:) .* tt;
Saso@1:   end
Saso@1: end