function P = tridikarpol(A)
n=max(size(A));
P=zeros(n+1);
P(1,1)=1;
P(2,1)=1;
P(2,2)=-A(1,1);
for k=3:(n+1), 
    s1=zeros(1,k);
    s2=zeros(1,k);
    s3=zeros(1,k);
    s1(1:k-1)=P(k-1,1:k-1);
    s2(2:k)=-A(k-1,k-1)*P(k-1,1:k-1);
    s3(3:k)=-A(k-1,k-2)^2*P(k-2,1:k-2);
    s=s1+s2+s3;
    P(k,1:k)=s;
end;
