syms x y
n=size(sed.K.s);
n=n(2);
i=1;m=0;
p=sym(0);
while i<=n
    k=sed.K.s(i);
    A=reshape(sed.x(m+1:m+k^2),k,k);
    j=1;
    z1=sym([1,zeros(1,k-1)]);
    z2=sym([1,zeros(1,k-1)]');
    while(j<k)
        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);
        j=j+1;
    end;
    m=m+k^2;
    ezsurf(p,[-1.2,1.2,-1.2,1.2]);
    set(gca,'NextPlot','add');
    p=p+sym(z1*A*(z2));
    i=i+1;
end
p