% GTD4
clc
randn('seed',0)
A = randn(20,6);
b = randn(20,1);
for r = [0 10 20]
    for p = [ 1 3 6]
        disp(' ' )
        disp(sprintf('r = %2d  p = %2d',r,p))
        format long
        x = pTLS(p,r,A,b)
    end
end
