function [vec1] = som_par(l)

d = log(2)/log(3);

for n=1:50
    som_par(n) = l^(2*complex(0,1)*n*pi/log(3))*epsi(n)/(d+2*n*complex(0,1)*pi/log(3))+l^(-2*complex(0,1)*n*pi/log(3))*epsi(-n)/(d-2*n*complex(0,1)*pi/log(3));
end

vec1 = som_par

end


