syms x y
n=6;
i=1;m=0;
q=sym(0);

A=sed.M{n};
m=size(A);
m=m(1);

z1=sym([1,zeros(1,m-1)]);
z2=sym([1,zeros(1,m-1)]');

for j=1:m-1
    %     if sed.pows(1,3)==0
    %     z1(1,j+1)=x^sed.pows(j,2)*y^sed.pows(j,3);
    %     z2(j+1,1)=x^sed.pows(j,2)*y^sed.pows(j,3);
    %     else,
    %     z1(1,j+1)=0;
    %     z2(j+1,1)=0;
    %     end;
    z1(1,j+1)=x^sed.pows(j,1)*y^sed.pows(j,2);
    z2(j+1,1)=x^sed.pows(j,1)*y^sed.pows(j,2);
end;

q=sym(z1*A*(z2));
i=i+1;
q