function [A] = constrpoly(T)
    n = size(T,1)
    A = zeros(n,3)
    A(:,1) = ones(n,1)
    A(:,2) = T
    A(:,3) = T.*T
endfunction