function z=bern_sum(x,b)
     [r,c]=size(b);
     r=max([r c]);
     z=sym(0);
     for k=2:r
     z=z+b(k)*(-1)^k/(k*(k-1)*x^(k-1));
end