% FKL The function to be transformed by the Kanrorovich-Lebedev
%     transformation
%
function y=fKL(x)
global a lambda 
y=ones(size(x));
%y=x;
%y=x.^lambda;
%y=exp(-x*cos(a));
%y=exp(-x*cosh(a));
%y=exp(-x);
%y=x.^lambda.*exp(-x);
%y=exp(-x)./sqrt(x);
%y=exp(-x*cosh(a))./sqrt(x);
%y=exp(-x-a)./((x+a).*sqrt(x));
%y=exp(-x-.5*a^2./x)./(2*x);
%y=exp(-a^2*x.^2);
%y=cos(x);
%y=sin(x);
%y=cos(x*sinh(a/2));
%y=sin(x*sinh(a/2));
%y=zeros(size(x));
%low=find(abs(x)<1e-10); high=find(abs(x)>=1e-10);
%xl=x(low); xh=x(high);
%if isempty(low)==0
%  y(low)=sinh(a/2);
%   y(low)=zeros(size(xl));
%   y(low)=.5*(sinh(a/2))^2*ones(size(xl));
%end
%if isempty(high)==0
%  y(high)=sin(xh*sinh(a/2))./xh;
%   y(high)=(1-cos(xh*sinh(a/2)))./xh;
%   y(high)=2*(sin(.5*sinh(a/2)*xh)).^2./xh.^2;
%end
%y=exp(-x*cosh(a)*cos(lambda)).*cos(x*sinh(a)*sin(lambda));
%y=bessel(0,x*sinh(a/2));
%y=x.^lambda.*bessel(lambda,x*sinh(a/2));
%y=besselk(lambda,x);
%y=exp(-x-a)./sqrt(x+a);
%y=erfc(sqrt(x));
%y=besselk(0,x);