function lambda = evcyl(x)
lambda=x./sqrt(1-1./(x.*besselj(0,x)./besselj(1,x)-1).^2);
end