please corrected attachment file

I wanted to draw streamlines between two concentric cylibders.

please corrected attachment file .

run("stream_H.m")

Unrecognized function or variable 'talbot_inversion'.

Error in stream_H (line 44)
mH1=talbot_inversion(psi_H,1);

Error in run (line 100)
evalin('caller', strcat(scriptStem, ';'));

[Torsten]

We miss the function "talbot_inversion" to run your code.

[dpb]

type stream_H

close all; clear all; clc;
xi=0.5^(1/2);
U_x=1;
j=0.1;s_1=3.0;s_2=s_1;sl=1;c=2;rho_s=1;
LL=xi; a=1;
c =-a/LL;
b =a/LL;
m =a*20; % NUMBER OF INTERVALS
%[x,y]=meshgrid([c:(b-c)/m:b],[c:(b-c)/m:b]');
[x,y]=meshgrid([-1/LL:1.01:1/LL]);
[I J]=find(sqrt(x.^2+y.^2)<(a-.25));
if ~isempty(I);
 x(I,J) = 0;
y(I,J) = 0;
end
rho=sqrt(x.^2+y.^2);
theta=atan2(y,x);
 for i1=1:length(x);
   for k1=1:length(x);
      if sqrt(x(i1,k1).^2+y(i1,k1).^2)>1./LL;
         rho(i1,k1)=0;
      end
   end
end
warning off
 L1 =@(S)( 1./2.*(2.*s_1.*S.*j-8.*s_1.*c.^2+2.*sl.^2+8.*s_1.*c+2.*S+2.*(S.^2.*j.^2.*s_1.^2-8.*S.*c.^2.*j.*s_1.^2+16.*c.^4.*s_1.^2+8.*S.*c.*j.*s_1.^2-2.*S.*j.*s_1.*sl.^2-32.*c.^3.*s_1.^2-8.*c.^2.*s_1.*sl.^2-2.*S.^2.*j.*s_1-8.*S.*c.^2.*s_1+16.*c.^2.*s_1.^2-8.*c.*s_1.*sl.^2+sl.^4-8.*S.*c.*s_1+2.*S.*sl.^2+S.^2).^(1./2)).^(1./2));
L2 =@(S)( 1./2.*(2.*s_1.*S.*j-8.*s_1.*c.^2+2.*sl.^2+8.*s_1.*c+2.*S-2.*(S.^2.*j.^2.*s_1.^2-8.*S.*c.^2.*j.*s_1.^2+16.*c.^4.*s_1.^2+8.*S.*c.*j.*s_1.^2-2.*S.*j.*s_1.*sl.^2-32.*c.^3.*s_1.^2-8.*c.^2.*s_1.*sl.^2-2.*S.^2.*j.*s_1-8.*S.*c.^2.*s_1+16.*c.^2.*s_1.^2-8.*c.*s_1.*sl.^2+sl.^4-8.*S.*c.*s_1+2.*S.*sl.^2+S.^2).^(1./2)).^(1./2));
AH =@(S)( -U_x.*(((L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(L1(S)+L2(S)).*(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L2(S))-besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*xi.*besseli(1,L2(S)./xi)-(L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L1(S))+besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)-L2(S).*((-2.*xi.^2.*(c-1).*L1(S).^2+((4.*c-4).*xi.^2-sl.^2-S).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L2(S))-besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)-(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*((-2.*c+2).*xi.^2+sl.^2+S).*L1(S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*besselk(0,L1(S)./xi)+(-2.*xi.^2.*(c-1).*L1(S).^2+((4.*c-4).*xi.^2-sl.^2-S).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*besselk(0,L2(S)./xi).*(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L1(S))+besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)))).*besseli(1,L1(S)./xi)+((L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(L1(S)+L2(S)).*(((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*xi.*besselk(1,L1(S)./xi)+((((4.*c-4).*xi.^2-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L2(S))-besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besseli(0,L1(S)./xi)-(((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*besselk(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*L2(S).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*besselk(0,L2(S)./xi).*((-2.*c+2).*xi.^2+sl.^2+S)).*L1(S)).*besseli(1,L2(S)./xi)+((((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*(L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)-L2(S).*(-(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*((-2.*c+2).*xi.^2+sl.^2+S).*L1(S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*besseli(0,L1(S)./xi)+((((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)-(((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)).*(-2.*xi.^2.*(c-1).*L1(S).^2+((4.*c-4).*xi.^2-sl.^2-S).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)))).*besselk(1,L1(S)./xi)+L1(S).*((-(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L1(S))+besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*besseli(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*L2(S).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*((-2.*c+2).*xi.^2+sl.^2+S).*besseli(0,L2(S)./xi)-(((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)./xi).*(((4.*c-4).*xi.^2-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S))).*besselk(1,L2(S)./xi)-2.*(((((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L2(S))-besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+besselk(0,L2(S)./xi).*(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L1(S))+besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S))).*besseli(0,L1(S)./xi)-((((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)-(((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)).*besselk(0,L1(S)./xi)).*(L1(S)-L2(S)).*L2(S).*(c-1).*(L1(S)+L2(S)).*xi))./(((L1(S)-L2(S)).*(-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)+4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L1(S)).*L1(S)).*(L1(S)+L2(S)).*xi.*besseli(1,L2(S)./xi)-(L1(S)-L2(S)).*(-4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)-4.*besseli(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)+(2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besseli(0,L2(S)).*L2(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)-((-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)+(2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L1(S)).*L1(S)).*besseli(0,L2(S)./xi)-(sl.^2-L2(S).^2+S).*L1(S).*(4.*besseli(1,L2(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(0,L2(S)./xi)+(sl.^2-L1(S).^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S)))).*besselk(0,L1(S)./xi)+((2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besseli(0,L2(S)).*L2(S)).*besselk(0,L2(S)./xi)-4.*besselk(1,L1(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*L2(S)).*besseli(1,L1(S)./xi)+((L1(S)-L2(S)).*(-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(0,L1(S)./xi)-4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(-2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*L1(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L1(S)./xi)+((-4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(2.*((-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^2+sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L1(S)).*L1(S)).*besseli(0,L1(S)./xi)+(-4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(2.*((-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^2+sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L1(S)).*besselk(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*((sl.^2-L1(S).^2+S).*L2(S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*besselk(0,L2(S)./xi)-4.*besselk(1,L2(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1)).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S)))).*L1(S)).*besseli(1,L2(S)./xi)+((L1(S)-L2(S)).*(4.*besseli(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(0,L1(S)./xi)-4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+(2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*L1(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)-L2(S).*(-(sl.^2-L2(S).^2+S).*(4.*besselk(1,L2(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besseli(0,L2(S)./xi)+4.*besseli(1,L2(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(0,L2(S)./xi)+(sl.^2-L1(S).^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S)))).*L1(S).*besseli(0,L1(S)./xi)+((-2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L1(S)).*besseli(0,L2(S)./xi)+((-2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L2(S))-(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L1(S)).*besselk(0,L2(S)./xi)+4.*besseli(1,L1(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S))).*besselk(1,L1(S)./xi)+(((-4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+(-2.*((-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^2+sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L1(S))-(sl.^2-L2(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S))-(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besseli(0,L2(S)).*L2(S)).*besseli(0,L1(S)./xi)+(sl.^2-L1(S).^2+S).*(-4.*besseli(1,L1(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L2(S).*besseli(0,L2(S)./xi)+((-2.*((-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^2+sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L1(S)).*besselk(0,L1(S)./xi)+4.*besseli(1,L2(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S)./xi)+4.*(L1(S)-L2(S)).*L2(S).*(c-1).*(L1(S)+L2(S)).*(((((-xi.^2.*L1(S).^2+xi.^2.*L2(S).^2).*besselk(1,L2(S))+1./2.*besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-1./2.*besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+(((-xi.^2.*L1(S).^2+xi.^2.*L2(S).^2).*besselk(1,L1(S))-1./2.*besselk(0,L1(S)).*L1(S).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))-1./2.*besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)./xi)-besselk(1,L1(S)).*xi.^2.*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)+((((-xi.^2.*L1(S).^2+xi.^2.*L2(S).^2).*besselk(1,L2(S))+1./2.*besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+1./2.*besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)./xi)+besselk(1,L2(S)).*xi.^2.*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+((((-xi.^2.*L1(S).^2+xi.^2.*L2(S).^2).*besseli(1,L2(S))-1./2.*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+1./2.*besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)-besseli(1,L1(S)).*xi.^2.*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)./xi)+besselk(0,L2(S)./xi).*besseli(1,L2(S)).*xi.^2.*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*xi).*L1(S)));
BH =@(S)( (((L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(-besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besselk(1,L1(S))+besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*(L1(S)+L2(S)).*xi.*besseli(1,L2(S)./xi)-(L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)-L2(S).*((-2.*xi.^2.*(c-1).*L1(S).^2+((4.*c-4).*xi.^2-sl.^2-S).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*(-besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besselk(1,L1(S))+besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)-(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*((-2.*c+2).*xi.^2+sl.^2+S).*L1(S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*besselk(0,L1(S)./xi)+(-2.*xi.^2.*(c-1).*L1(S).^2+((4.*c-4).*xi.^2-sl.^2-S).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*(besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)./xi)+4.*besselk(1,L1(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S))).*besseli(1,L1(S)./xi)+((besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*(L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(L1(S)+L2(S)).*xi.*besselk(1,L1(S)./xi)+((-besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besselk(1,L1(S))+besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*besseli(0,L1(S)./xi)-(besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*besselk(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*((sl.^2-L1(S).^2+S).*L2(S).*((-2.*c+2).*xi.^2+sl.^2+S).*besselk(0,L2(S)./xi)+4.*besselk(1,L2(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1)).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S)))).*L1(S)).*besseli(1,L2(S)./xi)+((L1(S)-L2(S)).*(L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)-L2(S).*(-(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*((-2.*c+2).*xi.^2+sl.^2+S).*L1(S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*besseli(0,L1(S)./xi)+(besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*(-2.*xi.^2.*(c-1).*L1(S).^2+((4.*c-4).*xi.^2-sl.^2-S).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*besseli(0,L2(S)./xi)-(-2.*xi.^2.*(c-1).*L1(S).^2+((4.*c-4).*xi.^2-sl.^2-S).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*(L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)-4.*besseli(1,L1(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S))).*besselk(1,L1(S)./xi)+((-(besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*besseli(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*L2(S).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*((-2.*c+2).*xi.^2+sl.^2+S).*besseli(0,L2(S)./xi)-(L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((-2.*c+2).*xi.^2+sl.^2+S)).*besselk(0,L1(S)./xi)-4.*besseli(1,L2(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S)./xi)-4.*(L1(S)-L2(S)).*L2(S).*(((-1./2.*besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besselk(1,L1(S))+1./2.*besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+(1./2.*besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S))+1./2.*besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)./xi)-besselk(1,L1(S)).*xi.^2.*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)+((-1./2.*besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))-1./2.*besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)./xi)+besselk(1,L2(S)).*xi.^2.*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+((1./2.*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))-1./2.*besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)-besseli(1,L1(S)).*xi.^2.*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)./xi)+besselk(0,L2(S)./xi).*besseli(1,L2(S)).*xi.^2.*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*(c-1).*(L1(S)+L2(S)).*xi).*L1(S)).*U_x./(((L1(S)-L2(S)).*(-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)+4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L1(S)).*L1(S)).*(L1(S)+L2(S)).*xi.*besseli(1,L2(S)./xi)-(L1(S)-L2(S)).*(-4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)-4.*besseli(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)+(2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besseli(0,L2(S)).*L2(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)-((-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)+(2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L1(S)).*L1(S)).*besseli(0,L2(S)./xi)-(sl.^2-L2(S).^2+S).*L1(S).*(4.*besseli(1,L2(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(0,L2(S)./xi)+(sl.^2-L1(S).^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S)))).*besselk(0,L1(S)./xi)+((2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besseli(0,L2(S)).*L2(S)).*besselk(0,L2(S)./xi)-4.*besselk(1,L1(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*L2(S)).*besseli(1,L1(S)./xi)+((L1(S)-L2(S)).*(-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(0,L1(S)./xi)-4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(-2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*L1(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L1(S)./xi)+((-4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(2.*((-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^2+sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L1(S)).*L1(S)).*besseli(0,L1(S)./xi)+(-4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(2.*((-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^2+sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L1(S)).*besselk(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*((sl.^2-L1(S).^2+S).*L2(S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*besselk(0,L2(S)./xi)-4.*besselk(1,L2(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1)).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S)))).*L1(S)).*besseli(1,L2(S)./xi)+((L1(S)-L2(S)).*(4.*besseli(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(0,L1(S)./xi)-4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+(2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*L1(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)-L2(S).*(-(sl.^2-L2(S).^2+S).*(4.*besselk(1,L2(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besseli(0,L2(S)./xi)+4.*besseli(1,L2(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(0,L2(S)./xi)+(sl.^2-L1(S).^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S)))).*L1(S).*besseli(0,L1(S)./xi)+((-2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L1(S)).*besseli(0,L2(S)./xi)+((-2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L2(S))-(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L1(S)).*besselk(0,L2(S)./xi)+4.*besseli(1,L1(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S))).*besselk(1,L1(S)./xi)+(((-4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+(-2.*((-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^2+sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L1(S))-(sl.^2-L2(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S))-(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besseli(0,L2(S)).*L2(S)).*besseli(0,L1(S)./xi)+(sl.^2-L1(S).^2+S).*(-4.*besseli(1,L1(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L2(S).*besseli(0,L2(S)./xi)+((-2.*((-sl.^2-S).*L1(S).^2-2.*xi.^2.*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^2+sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+(sl.^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L1(S)).*besselk(0,L1(S)./xi)+4.*besseli(1,L2(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S)./xi)+4.*(L1(S)-L2(S)).*L2(S).*(c-1).*(L1(S)+L2(S)).*(((((-xi.^2.*L1(S).^2+xi.^2.*L2(S).^2).*besselk(1,L2(S))+1./2.*besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-1./2.*besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+(((-xi.^2.*L1(S).^2+xi.^2.*L2(S).^2).*besselk(1,L1(S))-1./2.*besselk(0,L1(S)).*L1(S).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))-1./2.*besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)./xi)-besselk(1,L1(S)).*xi.^2.*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)+((((-xi.^2.*L1(S).^2+xi.^2.*L2(S).^2).*besselk(1,L2(S))+1./2.*besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+1./2.*besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)./xi)+besselk(1,L2(S)).*xi.^2.*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+((((-xi.^2.*L1(S).^2+xi.^2.*L2(S).^2).*besseli(1,L2(S))-1./2.*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+1./2.*besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)-besseli(1,L1(S)).*xi.^2.*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)./xi)+besselk(0,L2(S)./xi).*besseli(1,L2(S)).*xi.^2.*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*xi).*L1(S)));
CH =@(S)( 2.*(sl.^2-L2(S).^2+S).*((besselk(1,L2(S)).*(L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(L1(S)+L2(S)).*xi.*besseli(1,L2(S)./xi)-besseli(1,L2(S)).*(L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)-(besselk(1,L2(S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L2(S).^2-2.*(c-1).*(sl.^2+L1(S).^2+S).*xi.^2+(sl.^2+S).^2).*besseli(0,L2(S)./xi)+besseli(1,L2(S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L2(S).^2-2.*(c-1).*(sl.^2+L1(S).^2+S).*xi.^2+(sl.^2+S).^2).*besselk(0,L2(S)./xi)-2.*xi.^4.*(c-1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*L2(S)).*besseli(1,L1(S)./xi)+(besselk(1,L2(S)).*L1(S).*(-2.*xi.^2.*(c-1).*L2(S).^2-2.*(c-1).*(sl.^2-2.*L1(S).^2+S).*xi.^2+(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)+(sl.^2-L1(S).^2+S).*L2(S).*besseli(1,L1(S)).*((-2.*c+2).*xi.^2+sl.^2+S).*besselk(0,L2(S)./xi)-2.*(c-1).*(((-2.*L1(S).^2+2.*L2(S).^2).*besseli(1,L1(S))+L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))).*xi.^4).*besseli(1,L2(S)./xi)+(-besseli(1,L2(S)).*L1(S).*(-2.*xi.^2.*(c-1).*L2(S).^2-2.*(c-1).*(sl.^2-2.*L1(S).^2+S).*xi.^2+(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)+(sl.^2-L1(S).^2+S).*L2(S).*besseli(1,L1(S)).*((-2.*c+2).*xi.^2+sl.^2+S).*besseli(0,L2(S)./xi)-2.*(((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L1(S))-L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))).*(c-1).*xi.^4).*besselk(1,L2(S)./xi)+2.*(-(besselk(1,L2(S)).*(L1(S).^2-L2(S).^2).*besseli(0,L2(S)./xi)+besseli(1,L2(S)).*(L1(S).^2-L2(S).^2).*besselk(0,L2(S)./xi)+xi.^2.*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*L1(S).*besseli(0,L1(S)./xi)+((((-2.*L1(S).^2+2.*L2(S).^2).*besseli(1,L1(S))+L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))).*besseli(0,L2(S)./xi)-(((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L1(S))-L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))).*besselk(0,L2(S)./xi)).*xi.^2).*L2(S).*(c-1).*xi).*U_x./(((L1(S)-L2(S)).*(4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)+(2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L1(S)).*L1(S)).*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L2(S)).*L2(S)).*(L1(S)+L2(S)).*xi.*besseli(1,L2(S)./xi)-(L1(S)-L2(S)).*(-4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)-4.*besseli(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)+(2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besseli(0,L2(S)).*L2(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)-((-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)+(2.*((-sl.^2-S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L1(S)).*L1(S)).*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L2(S)).*L2(S)).*besseli(0,L2(S)./xi)+(-4.*besseli(1,L2(S)).*xi.^2.*L1(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)+(2.*((-sl.^2-S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besseli(0,L2(S)).*L2(S)).*besselk(0,L2(S)./xi)-(sl.^2-L1(S).^2+S).*((sl.^2-L2(S).^2+S).*L1(S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*besselk(0,L1(S)./xi)+4.*besselk(1,L1(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1)).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S)))).*L2(S)).*besseli(1,L1(S)./xi)+(L1(S).*(-4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)-4.*besselk(1,L2(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(1,L1(S)./xi)+(2.*(-2.*xi.^2.*(c-1).*L2(S).^2+((2.*c-2).*xi.^2+sl.^2-L1(S).^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)).*L1(S)).*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)).*L2(S)).*besseli(0,L1(S)./xi)+(sl.^2-L1(S).^2+S).*L2(S).*(-4.*besseli(1,L1(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(1,L1(S)./xi)+(-4.*besseli(1,L1(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L1(S)).*besselk(0,L2(S)./xi)+((-2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*L1(S)).*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L2(S)).*L2(S).*besseli(1,L1(S))).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L1(S)./xi)-(((-2.*(-2.*xi.^2.*(c-1).*L2(S).^2+((2.*c-2).*xi.^2+sl.^2-L1(S).^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*L1(S)).*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)).*L2(S).*besseli(1,L1(S))).*besselk(0,L1(S)./xi)+4.*besselk(1,L2(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S)).*besseli(1,L2(S)./xi)+(-(4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)-4.*besseli(1,L2(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(1,L1(S)./xi)+(2.*(-2.*xi.^2.*(c-1).*L2(S).^2+((2.*c-2).*xi.^2+sl.^2-L1(S).^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besseli(0,L2(S)).*L2(S)).*L1(S).*besseli(0,L1(S)./xi)+(sl.^2-L1(S).^2+S).*L2(S).*(-4.*besseli(1,L1(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(1,L1(S)./xi)+(-4.*besseli(1,L1(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L1(S)).*besseli(0,L2(S)./xi)+(L1(S)-L2(S)).*((2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besseli(0,L2(S)).*L2(S).*besseli(1,L1(S))).*(L1(S)+L2(S)).*xi.*besselk(1,L1(S)./xi)-(((2.*(-2.*xi.^2.*(c-1).*L2(S).^2+((2.*c-2).*xi.^2+sl.^2-L1(S).^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besseli(0,L2(S)).*L2(S).*besseli(1,L1(S))).*besselk(0,L1(S)./xi)-4.*besseli(1,L2(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S)).*besselk(1,L2(S)./xi)-L2(S).*(-L1(S).*(-2.*(L1(S)-L2(S)).*(-2.*besselk(1,L2(S)).*xi.*(sl.^2-L2(S).^2+S).*besselk(1,L1(S)./xi)+((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))+besselk(0,L1(S)).*L1(S).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))-besselk(0,L2(S)).*besselk(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*(c-1).*(L1(S)+L2(S)).*xi.*besseli(0,L2(S)./xi)-2.*(L1(S)-L2(S)).*(c-1).*(L1(S)+L2(S)).*(-2.*besseli(1,L2(S)).*xi.*(sl.^2-L2(S).^2+S).*besselk(1,L1(S)./xi)+((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))+besselk(0,L1(S)).*L1(S).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*xi.*besselk(0,L2(S)./xi)+(sl.^2-L1(S).^2+S).*((sl.^2-L2(S).^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*besselk(1,L1(S)./xi)-4.*besselk(1,L1(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1)).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S)))).*besseli(0,L1(S)./xi)+(((-2.*((-sl.^2-S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*L1(S)).*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L2(S)).*L2(S).*besseli(1,L1(S))).*besselk(1,L1(S)./xi)-4.*(L1(S)-L2(S)).*((((-xi.^2.*L1(S).^2+xi.^2.*L2(S).^2).*besseli(1,L1(S))+1./2.*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+1./2.*besselk(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)./xi)+besselk(1,L2(S)).*xi.^2.*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S).*(c-1).*(L1(S)+L2(S)).*xi).*besseli(0,L2(S)./xi)+(((-2.*((-sl.^2-S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*L1(S)).*besseli(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besseli(0,L2(S)).*L2(S).*besseli(1,L1(S))).*besselk(1,L1(S)./xi)-4.*(L1(S)-L2(S)).*((((-xi.^2.*L1(S).^2+xi.^2.*L2(S).^2).*besseli(1,L1(S))+1./2.*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))-1./2.*besseli(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)./xi)+besseli(1,L2(S)).*xi.^2.*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S).*(c-1).*(L1(S)+L2(S)).*xi).*besselk(0,L2(S)./xi)+4.*besseli(1,L1(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L1(S)./xi).*L1(S)+xi.*besselk(1,L1(S)./xi)).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S))));
DH =@(S)( -2.*(sl.^2-L1(S).^2+S).*U_x.*((besselk(1,L1(S)).*(L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(L1(S)+L2(S)).*xi.*besseli(1,L1(S)./xi)-(L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*besseli(1,L1(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L1(S)./xi)+L1(S).*(besselk(1,L1(S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L1(S).^2-2.*(c-1).*(sl.^2+L2(S).^2+S).*xi.^2+(sl.^2+S).^2).*besseli(0,L1(S)./xi)+besseli(1,L1(S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L1(S).^2-2.*(c-1).*(sl.^2+L2(S).^2+S).*xi.^2+(sl.^2+S).^2).*besselk(0,L1(S)./xi)-2.*xi.^4.*(c-1).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S))).*besseli(1,L2(S)./xi)+(-besselk(1,L1(S)).*L2(S).*(-2.*xi.^2.*(c-1).*L1(S).^2-2.*(c-1).*(sl.^2-2.*L2(S).^2+S).*xi.^2+(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)-(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*((-2.*c+2).*xi.^2+sl.^2+S).*L1(S).*besselk(0,L1(S)./xi)+2.*(((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S))).*(c-1).*xi.^4).*besseli(1,L1(S)./xi)+(L2(S).*(-2.*xi.^2.*(c-1).*L1(S).^2-2.*(c-1).*(sl.^2-2.*L2(S).^2+S).*xi.^2+(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besseli(1,L1(S)).*besseli(0,L2(S)./xi)-(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*((-2.*c+2).*xi.^2+sl.^2+S).*L1(S).*besseli(0,L1(S)./xi)-2.*(((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*(c-1).*xi.^4).*besselk(1,L1(S)./xi)+2.*L1(S).*(c-1).*xi.*(L2(S).*((-L1(S).^2+L2(S).^2).*besselk(1,L1(S)).*besseli(0,L1(S)./xi)+(-L1(S).^2+L2(S).^2).*besseli(1,L1(S)).*besselk(0,L1(S)./xi)+xi.^2.*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)-((((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S))).*besseli(0,L1(S)./xi)+(((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)./xi)).*xi.^2))./(((L1(S)-L2(S)).*(-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)+4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L1(S)).*L1(S)).*(L1(S)+L2(S)).*xi.*besseli(1,L1(S)./xi)+(L1(S)-L2(S)).*(-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(0,L1(S)./xi)-4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(-2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*L1(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L1(S)./xi)+((-4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(2.*((-sl.^2-S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L1(S)).*L1(S)).*besseli(0,L1(S)./xi)+(-4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(2.*((-sl.^2-S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*L1(S)).*besselk(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*((sl.^2-L1(S).^2+S).*L2(S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*besselk(0,L2(S)./xi)-4.*besselk(1,L2(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1)).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S)))).*L1(S)).*besseli(1,L2(S)./xi)+(-L2(S).*(-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)-4.*besselk(1,L1(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(1,L2(S)./xi)+(2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)).*L1(S)).*besseli(0,L2(S)./xi)+(sl.^2-L2(S).^2+S).*(4.*besseli(1,L2(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(1,L2(S)./xi)+L2(S).*(4.*besseli(1,L2(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(0,L2(S)./xi)+(sl.^2-L1(S).^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))))).*L1(S).*besselk(0,L1(S)./xi)-(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.*((2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besseli(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L1(S)).*L1(S)).*besselk(1,L2(S)./xi)-(((2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)).*L1(S)).*besselk(0,L2(S)./xi)-4.*besselk(1,L1(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*L2(S)).*besseli(1,L1(S)./xi)+(-(-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(0,L1(S)./xi)+4.*besseli(1,L1(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(1,L2(S)./xi)+(-2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*L1(S)).*L2(S).*besseli(0,L2(S)./xi)+(sl.^2-L2(S).^2+S).*(4.*besseli(1,L2(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(1,L2(S)./xi)+L2(S).*(4.*besseli(1,L2(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besselk(0,L2(S)./xi)+(sl.^2-L1(S).^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))))).*L1(S).*besseli(0,L1(S)./xi)+(L1(S)-L2(S)).*((2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*L1(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)+(((2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*L1(S)).*besselk(0,L2(S)./xi)-4.*besseli(1,L1(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*L2(S)).*besselk(1,L1(S)./xi)+(L2(S).*(-2.*(2.*besselk(1,L1(S)).*xi.*(sl.^2-L1(S).^2+S).*besselk(1,L2(S)./xi)+((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))-besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*(L1(S)-L2(S)).*(c-1).*(L1(S)+L2(S)).*xi.*besseli(0,L1(S)./xi)+2.*(-2.*besseli(1,L1(S)).*xi.*(sl.^2-L1(S).^2+S).*besselk(1,L2(S)./xi)+((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*(L1(S)-L2(S)).*(c-1).*(L1(S)+L2(S)).*xi.*besselk(0,L1(S)./xi)+((sl.^2-L1(S).^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*besselk(1,L2(S)./xi)+4.*besselk(1,L2(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1)).*(sl.^2-L2(S).^2+S).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S)))).*besseli(0,L2(S)./xi)+(((-2.*((-sl.^2-S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besseli(0,L2(S)).*L2(S)).*besselk(1,L1(S))-(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L1(S)).*L1(S)).*besselk(1,L2(S)./xi)-2.*(L1(S)-L2(S)).*L2(S).*((((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besselk(0,L1(S)).*besseli(1,L2(S)).*L1(S).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)+2.*besselk(1,L1(S)).*xi.^2.*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*(c-1).*(L1(S)+L2(S)).*xi).*besseli(0,L1(S)./xi)+(((-2.*((-sl.^2-S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*L1(S)).*besselk(1,L2(S)./xi)-4.*(L1(S)-L2(S)).*L2(S).*((((xi.^2.*L1(S).^2-xi.^2.*L2(S).^2).*besseli(1,L2(S))+1./2.*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-1./2.*besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)+besseli(1,L1(S)).*xi.^2.*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*(c-1).*(L1(S)+L2(S)).*xi).*besselk(0,L1(S)./xi)+4.*besseli(1,L2(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L2(S)./xi).*L2(S)+xi.*besselk(1,L2(S)./xi)).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S)));
EH =@(S)( -2.*(sl.^2-L2(S).^2+S).*((besselk(1,L2(S)).*(L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(L1(S)+L2(S)).*xi.*besseli(1,L2(S)./xi)-besseli(1,L2(S)).*(L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)-(besselk(1,L2(S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L2(S).^2-2.*(c-1).*(sl.^2+L1(S).^2+S).*xi.^2+(sl.^2+S).^2).*besseli(0,L2(S)./xi)+besseli(1,L2(S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L2(S).^2-2.*(c-1).*(sl.^2+L1(S).^2+S).*xi.^2+(sl.^2+S).^2).*besselk(0,L2(S)./xi)-2.*xi.^4.*(c-1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*L2(S)).*besselk(1,L1(S)./xi)+(-besselk(1,L2(S)).*L1(S).*(-2.*xi.^2.*(c-1).*L2(S).^2-2.*(c-1).*(sl.^2-2.*L1(S).^2+S).*xi.^2+(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)./xi)+(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*L2(S).*((-2.*c+2).*xi.^2+sl.^2+S).*besselk(0,L2(S)./xi)+2.*(c-1).*(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L1(S))+besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))-besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besselk(1,L1(S))).*xi.^4).*besseli(1,L2(S)./xi)+(besseli(1,L2(S)).*L1(S).*(-2.*xi.^2.*(c-1).*L2(S).^2-2.*(c-1).*(sl.^2-2.*L1(S).^2+S).*xi.^2+(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)./xi)+(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*L2(S).*((-2.*c+2).*xi.^2+sl.^2+S).*besseli(0,L2(S)./xi)-2.*(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L1(S))+besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*(c-1).*xi.^4).*besselk(1,L2(S)./xi)-2.*(-(besselk(1,L2(S)).*(L1(S).^2-L2(S).^2).*besseli(0,L2(S)./xi)+besseli(1,L2(S)).*(L1(S).^2-L2(S).^2).*besselk(0,L2(S)./xi)+xi.^2.*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*L1(S).*besselk(0,L1(S)./xi)+((((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L1(S))+besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))-besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besselk(1,L1(S))).*besseli(0,L2(S)./xi)+besselk(0,L2(S)./xi).*(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L1(S))+besselk(0,L1(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S))).*xi.^2).*L2(S).*(c-1).*xi).*U_x./(((L1(S)-L2(S)).*(-4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(0,L1(S)./xi)+(-2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*L1(S)).*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L2(S)).*L2(S).*besseli(1,L1(S))).*(L1(S)+L2(S)).*xi.*besseli(1,L2(S)./xi)+(L1(S)-L2(S)).*(-4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+4.*besseli(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(0,L1(S)./xi)+(2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besseli(0,L2(S)).*L2(S).*besseli(1,L1(S))).*(L1(S)+L2(S)).*xi.*besselk(1,L2(S)./xi)-L2(S).*((-4.*besselk(1,L2(S)).*xi.^2.*L1(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(0,L1(S)./xi)+(-2.*((-sl.^2-S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*L1(S)).*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L2(S)).*L2(S).*besseli(1,L1(S))).*besseli(0,L2(S)./xi)+(-4.*besseli(1,L2(S)).*xi.^2.*L1(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(0,L1(S)./xi)+(-2.*((-sl.^2-S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*L1(S)).*besseli(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besseli(0,L2(S)).*L2(S).*besseli(1,L1(S))).*besselk(0,L2(S)./xi)-(sl.^2-L1(S).^2+S).*((sl.^2-L2(S).^2+S).*L1(S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*besseli(0,L1(S)./xi)-4.*besseli(1,L1(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1)).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))))).*besselk(1,L1(S)./xi)+(-(4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+4.*besselk(1,L2(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(1,L1(S)./xi)+(-2.*(-2.*xi.^2.*(c-1).*L2(S).^2+((2.*c-2).*xi.^2+sl.^2-L1(S).^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*L1(S)).*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)).*L2(S).*besseli(1,L1(S))).*L1(S).*besselk(0,L1(S)./xi)+(sl.^2-L1(S).^2+S).*L2(S).*(4.*besselk(1,L1(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besseli(1,L1(S)./xi)+(-4.*besselk(1,L1(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besseli(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L1(S)).*besselk(0,L2(S)./xi)+((2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L1(S)).*L1(S)).*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L2(S)).*L2(S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.*besseli(1,L1(S)./xi)+(((2.*(-2.*xi.^2.*(c-1).*L2(S).^2+((2.*c-2).*xi.^2+sl.^2-L1(S).^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)).*L1(S)).*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)).*L2(S)).*besseli(0,L1(S)./xi)-4.*besselk(1,L2(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S)).*besseli(1,L2(S)./xi)+(-(4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)-4.*besseli(1,L2(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(1,L1(S)./xi)+(2.*(-2.*xi.^2.*(c-1).*L2(S).^2+((2.*c-2).*xi.^2+sl.^2-L1(S).^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besseli(0,L2(S)).*L2(S).*besseli(1,L1(S))).*L1(S).*besselk(0,L1(S)./xi)+(sl.^2-L1(S).^2+S).*L2(S).*(4.*besselk(1,L1(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besseli(1,L1(S)./xi)+(-4.*besselk(1,L1(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besseli(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S)).*L1(S)).*besseli(0,L2(S)./xi)-(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.*((2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S)./xi)-(((2.*(-2.*xi.^2.*(c-1).*L2(S).^2+((2.*c-2).*xi.^2+sl.^2-L1(S).^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L2(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L1(S).^2+S)).*besseli(0,L2(S)).*L2(S)).*besseli(0,L1(S)./xi)-4.*besseli(1,L2(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S)).*besselk(1,L2(S)./xi)-(-L1(S).*(2.*(L1(S)-L2(S)).*(2.*besselk(1,L2(S)).*xi.*(sl.^2-L2(S).^2+S).*besseli(1,L1(S)./xi)+((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))+L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*(c-1).*(L1(S)+L2(S)).*xi.*besseli(0,L2(S)./xi)-2.*(L1(S)-L2(S)).*(c-1).*(L1(S)+L2(S)).*(-2.*besseli(1,L2(S)).*xi.*(sl.^2-L2(S).^2+S).*besseli(1,L1(S)./xi)+((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))-L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besseli(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*xi.*besselk(0,L2(S)./xi)+((sl.^2-L2(S).^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*besseli(1,L1(S)./xi)-4.*besseli(1,L1(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1)).*(sl.^2-L1(S).^2+S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S)))).*besselk(0,L1(S)./xi)+(((2.*((-sl.^2-S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L1(S)).*L1(S)).*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S)./xi)-4.*(L1(S)-L2(S)).*L1(S).*((((-xi.^2.*L1(S).^2+xi.^2.*L2(S).^2).*besselk(1,L1(S))-1./2.*besselk(0,L1(S)).*L1(S).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+1./2.*besselk(0,L2(S)).*besselk(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)+besselk(1,L2(S)).*xi.^2.*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*(c-1).*(L1(S)+L2(S)).*xi).*besseli(0,L2(S)./xi)+(((2.*((-sl.^2-S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*besselk(1,L1(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L2(S).^2+2.*(c-1).*(sl.^2-L1(S).^2+S).*xi.^4+((2.*c-2).*L1(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S)./xi)-4.*((((-xi.^2.*L1(S).^2+xi.^2.*L2(S).^2).*besselk(1,L1(S))-1./2.*besselk(0,L1(S)).*L1(S).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))-1./2.*besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)+besseli(1,L2(S)).*xi.^2.*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*(L1(S)-L2(S)).*L1(S).*(c-1).*(L1(S)+L2(S)).*xi).*besselk(0,L2(S)./xi)+4.*besselk(1,L1(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besseli(0,L1(S)./xi).*L1(S)-xi.*besseli(1,L1(S)./xi)).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*L2(S)));
FH =@(S)( 2.*(sl.^2-L1(S).^2+S).*U_x.*((besselk(1,L1(S)).*(L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*(L1(S)+L2(S)).*xi.*besseli(1,L1(S)./xi)-(L1(S)-L2(S)).*((-4.*c+4).*xi.^2+sl.^2+S).*besseli(1,L1(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L1(S)./xi)+L1(S).*(besselk(1,L1(S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L1(S).^2-2.*(c-1).*(sl.^2+L2(S).^2+S).*xi.^2+(sl.^2+S).^2).*besseli(0,L1(S)./xi)+besseli(1,L1(S)).*(((4.*c-4).*xi.^2-sl.^2-S).*L1(S).^2-2.*(c-1).*(sl.^2+L2(S).^2+S).*xi.^2+(sl.^2+S).^2).*besselk(0,L1(S)./xi)-2.*xi.^4.*(c-1).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S))).*besselk(1,L2(S)./xi)+(besselk(1,L1(S)).*L2(S).*(-2.*xi.^2.*(c-1).*L1(S).^2-2.*(c-1).*(sl.^2-2.*L2(S).^2+S).*xi.^2+(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)-(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*((-2.*c+2).*xi.^2+sl.^2+S).*L1(S).*besselk(0,L1(S)./xi)+2.*(c-1).*(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L2(S))-besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*xi.^4).*besseli(1,L1(S)./xi)+(-L2(S).*(-2.*xi.^2.*(c-1).*L1(S).^2-2.*(c-1).*(sl.^2-2.*L2(S).^2+S).*xi.^2+(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besseli(1,L1(S)).*besselk(0,L2(S)./xi)-(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*((-2.*c+2).*xi.^2+sl.^2+S).*L1(S).*besseli(0,L1(S)./xi)+2.*(c-1).*(((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*xi.^4).*besselk(1,L1(S)./xi)-2.*L1(S).*(c-1).*xi.*(L2(S).*((-L1(S).^2+L2(S).^2).*besselk(1,L1(S)).*besseli(0,L1(S)./xi)+(-L1(S).^2+L2(S).^2).*besseli(1,L1(S)).*besselk(0,L1(S)./xi)+xi.^2.*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)+xi.^2.*((((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L2(S))-besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S)).*besseli(0,L1(S)./xi)-besselk(0,L1(S)./xi).*(((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)))))./((-(L1(S)-L2(S)).*(-4.*besseli(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)-4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+(2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besseli(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L1(S)).*L1(S)).*(L1(S)+L2(S)).*xi.*besseli(1,L1(S)./xi)+(L1(S)-L2(S)).*(4.*besseli(1,L2(S)).*xi.^2.*L1(S).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(0,L1(S)./xi)-4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+(2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*L1(S)).*(L1(S)+L2(S)).*xi.*besselk(1,L1(S)./xi)+((-4.*besselk(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+(-2.*((-sl.^2-S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besseli(0,L2(S)).*L2(S)).*besselk(1,L1(S))-(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L1(S)).*L1(S)).*besseli(0,L1(S)./xi)+(-4.*besseli(1,L1(S)).*xi.^2.*L2(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+(-2.*((-sl.^2-S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*L1(S)).*besselk(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*((sl.^2-L1(S).^2+S).*L2(S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*besseli(0,L2(S)./xi)+4.*besseli(1,L2(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1)).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S)))).*L1(S)).*besselk(1,L2(S)./xi)+(-L2(S).*(-4.*besseli(1,L2(S)).*xi.^2.*L1(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besselk(0,L1(S)./xi)-4.*besselk(1,L1(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(1,L2(S)./xi)+(2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)).*L1(S)).*besselk(0,L2(S)./xi)+(sl.^2-L2(S).^2+S).*(-4.*besselk(1,L2(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besseli(1,L2(S)./xi)+L2(S).*(4.*besselk(1,L2(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besseli(0,L2(S)./xi)+(sl.^2-L1(S).^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))))).*L1(S).*besselk(0,L1(S)./xi)+((2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L1(S)).*L1(S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.*besseli(1,L2(S)./xi)-L2(S).*(((2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)).*L1(S)).*besseli(0,L2(S)./xi)-4.*besselk(1,L1(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S))).*besseli(1,L1(S)./xi)+(L2(S).*(4.*besseli(1,L2(S)).*xi.^2.*L1(S).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L2(S).^2+S).*besseli(0,L1(S)./xi)-4.*besseli(1,L1(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(sl.^2-L1(S).^2+S).*besseli(1,L2(S)./xi)+(2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besseli(1,L2(S))+(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)).*L2(S)).*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*L1(S)).*besselk(0,L2(S)./xi)+(sl.^2-L2(S).^2+S).*(-4.*besselk(1,L2(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besseli(1,L2(S)./xi)+L2(S).*(4.*besselk(1,L2(S)).*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*besseli(0,L2(S)./xi)+(sl.^2-L1(S).^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))))).*L1(S).*besseli(0,L1(S)./xi)+((-2.*xi.^2.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*((sl.^2+S+4.*c-4).*xi.^2-sl.^2-S).*L1(S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.*besseli(1,L2(S)./xi)-(((-2.*(-2.*xi.^2.*(c-1).*L1(S).^2+(-L2(S).^2+(2.*c-2).*xi.^2+sl.^2+S).*(sl.^2+S)).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))+(sl.^2-L1(S).^2+S).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S))+besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(-2.*xi.^2.*(xi-1).*(xi+1).*(c-1).*L1(S).^2+2.*(c-1).*(sl.^2+S).*xi.^4+((-sl.^2-S-4.*c+4).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).*(sl.^2-L2(S).^2+S)).*L1(S)).*besseli(0,L2(S)./xi)+4.*besseli(1,L1(S)).*xi.^4.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*L2(S)).*besselk(1,L1(S)./xi)+(L2(S).*(-2.*(2.*besselk(1,L1(S)).*xi.*(sl.^2-L1(S).^2+S).*besseli(1,L2(S)./xi)+((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besselk(0,L1(S)).*besseli(1,L2(S)).*L1(S).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*(L1(S)-L2(S)).*(c-1).*(L1(S)+L2(S)).*xi.*besseli(0,L1(S)./xi)-4.*(besseli(1,L1(S)).*xi.*(sl.^2-L1(S).^2+S).*besseli(1,L2(S)./xi)+((xi.^2.*L1(S).^2-xi.^2.*L2(S).^2).*besseli(1,L2(S))+1./2.*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-1./2.*besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*(L1(S)-L2(S)).*(c-1).*(L1(S)+L2(S)).*xi.*besselk(0,L1(S)./xi)+(sl.^2-L2(S).^2+S).*((sl.^2-L1(S).^2+S).*((2.*c-2).*xi.^4+(sl.^2+S+2.*c-2).*xi.^2-sl.^2-S).*besseli(1,L2(S)./xi)+4.*besseli(1,L2(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1)).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S)))).*besselk(0,L2(S)./xi)+(((2.*((-sl.^2-S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L2(S)).*L2(S)).*besselk(1,L1(S))+(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L1(S)).*L1(S)).*besseli(1,L2(S)./xi)-2.*(L1(S)-L2(S)).*L2(S).*((((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))-besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besselk(0,L1(S)).*besselk(1,L2(S)).*L1(S).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+2.*besselk(1,L1(S)).*xi.^2.*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*(c-1).*(L1(S)+L2(S)).*xi).*besseli(0,L1(S)./xi)+(((2.*((-sl.^2-S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^2+(sl.^2+S).^2).*(L1(S)-L2(S)).*(L1(S)+L2(S)).*xi.^2.*besselk(1,L2(S))-(sl.^2-L1(S).^2+S).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*besselk(0,L2(S)).*L2(S)).*besseli(1,L1(S))-besseli(0,L1(S)).*(sl.^2-L2(S).^2+S).*besselk(1,L2(S)).*(((-sl.^2-S-4.*c+4).*xi.^2+sl.^2+S).*L1(S).^2+2.*(c-1).*(sl.^2-L2(S).^2+S).*xi.^4+((2.*c-2).*L2(S).^2+(sl.^2+S).*(sl.^2+S+2.*c-2)).*xi.^2-(sl.^2+S).^2).*L1(S)).*besseli(1,L2(S)./xi)-4.*(L1(S)-L2(S)).*L2(S).*((((xi.^2.*L1(S).^2-xi.^2.*L2(S).^2).*besselk(1,L2(S))-1./2.*besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-1./2.*besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+besseli(1,L1(S)).*xi.^2.*(besselk(0,L2(S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besselk(1,L2(S))).*(sl.^2-L1(S).^2+S)).*(c-1).*(L1(S)+L2(S)).*xi).*besselk(0,L1(S)./xi)+4.*besselk(1,L2(S)).*xi.^3.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(c-1).*(besseli(0,L2(S)./xi).*L2(S)-xi.*besseli(1,L2(S)./xi)).*(besselk(0,L1(S)).*besseli(1,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S)));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
AK =@(S)( -U_x.*(((L1(S)+L2(S)).*(sl.^2+S).*(((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(sl.^2-L2(S).^2+S)).*(L1(S)-L2(S)).*xi.*besseli(1,L2(S)./xi)+(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L1(S))+L1(S).*besselk(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*(L1(S)+L2(S)).*(sl.^2+S).*(L1(S)-L2(S)).*xi.*besselk(1,L2(S)./xi)-((((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+L1(S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)-(((2.*L1(S).^2-2.*L2(S).^2).*besselk(1,L1(S))+L1(S).*besselk(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besseli(1,L1(S)./xi)+(-(((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*(L1(S)+L2(S)).*(sl.^2+S).*(L1(S)-L2(S)).*xi.*besselk(1,L1(S)./xi)+(sl.^2-L1(S).^2+S).*((((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(0,L1(S)./xi)+(((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)./xi)-besselk(0,L2(S)./xi).*L2(S).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S).*L2(S).^2).*besseli(1,L2(S)./xi)+(-(L1(S)+L2(S)).*(sl.^2+S).*(((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*(L1(S)-L2(S)).*xi.*besselk(1,L2(S)./xi)+(-L1(S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+(((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)-(((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besselk(1,L1(S)./xi)-((((-2.*L1(S).^2+2.*L2(S).^2).*besselk(1,L1(S))-L1(S).*besselk(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))-besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)+L2(S).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S).*besseli(0,L2(S)./xi)-besselk(0,L1(S)./xi).*(((2.*L1(S).^2-2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S))).*(sl.^2-L1(S).^2+S).*besselk(1,L2(S)./xi).*L1(S).*L2(S).^2)./(((-2.*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+2.*besselk(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(sl.^2+S).*(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besseli(1,L2(S)./xi)+2.*(besselk(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+besseli(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+1./2.*(sl.^2+S).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))+L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L2(S)./xi)-((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+L1(S).*(xi-1).*(xi+1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)-(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))+L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besseli(1,L1(S)./xi)+(-(2.*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+2.*besseli(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*(sl.^2+S)).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L1(S)./xi)+(sl.^2-L1(S).^2+S).*((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L1(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)./xi)-besselk(0,L2(S)./xi).*L2(S).*(xi-1).*(xi+1).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S).*L2(S).^2).*besseli(1,L2(S)./xi)+(-2.*(-besseli(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+besseli(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+1./2.*(sl.^2+S).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L2(S)./xi)+(-L1(S).*(xi-1).*(xi+1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)-(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besselk(1,L1(S)./xi)-(sl.^2-L1(S).^2+S).*besselk(1,L2(S)./xi).*L1(S).*((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))-L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))-besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)+L2(S).*(xi-1).*(xi+1).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S).*besseli(0,L2(S)./xi)-besselk(0,L1(S)./xi).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*L2(S).^2));
BK =@(S)( (((L1(S)+L2(S)).*(sl.^2+S).*(besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(sl.^2-L2(S).^2+S)).*(L1(S)-L2(S)).*xi.*besseli(1,L2(S)./xi)+(L1(S).*besselk(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*(L1(S)+L2(S)).*(sl.^2+S).*(L1(S)-L2(S)).*xi.*besselk(1,L2(S)./xi)-((besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+L1(S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)-(L1(S).*besselk(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besseli(1,L1(S)./xi)+(-(L1(S)+L2(S)).*(sl.^2+S).*(L1(S)-L2(S)).*xi.*(besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(1,L1(S)./xi)+(sl.^2-L1(S).^2+S).*((besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(0,L1(S)./xi)+(besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)./xi)-besselk(0,L2(S)./xi).*L2(S).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S).*L2(S).^2).*besseli(1,L2(S)./xi)+(-(L1(S)+L2(S)).*(L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*(sl.^2+S).*(L1(S)-L2(S)).*xi.*besselk(1,L2(S)./xi)+(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S).*(-L1(S).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+(besselk(0,L2(S)).*L2(S).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)-(L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi))).*besselk(1,L1(S)./xi)-(sl.^2-L1(S).^2+S).*besselk(1,L2(S)./xi).*L1(S).*L2(S).^2.*((-L1(S).*besselk(0,L1(S)).*(sl.^2-L2(S).^2+S).*besseli(1,L2(S))-besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)+L2(S).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S).*besseli(0,L2(S)./xi)-(L2(S).*besseli(0,L2(S)).*(sl.^2-L1(S).^2+S).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)./xi))).*U_x./(((-2.*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+2.*besselk(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(sl.^2+S).*(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besseli(1,L2(S)./xi)+2.*(besselk(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+besseli(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+1./2.*(sl.^2+S).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))+L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L2(S)./xi)-((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+L1(S).*(xi-1).*(xi+1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)-(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))+L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besseli(1,L1(S)./xi)+(-(2.*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+2.*besseli(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*(sl.^2+S)).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L1(S)./xi)+(sl.^2-L1(S).^2+S).*((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L1(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)./xi)-besselk(0,L2(S)./xi).*L2(S).*(xi-1).*(xi+1).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S).*L2(S).^2).*besseli(1,L2(S)./xi)+(-2.*(-besseli(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+besseli(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+1./2.*(sl.^2+S).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L2(S)./xi)+(-L1(S).*(xi-1).*(xi+1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)-(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besselk(1,L1(S)./xi)-(sl.^2-L1(S).^2+S).*besselk(1,L2(S)./xi).*L1(S).*((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))-L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))-besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)+L2(S).*(xi-1).*(xi+1).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S).*besseli(0,L2(S)./xi)-besselk(0,L1(S)./xi).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*L2(S).^2));
CK =@(S)( -2.*(sl.^2-L2(S).^2+S).*U_x.*((besselk(1,L2(S)).*xi.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L2(S)./xi)-besseli(1,L2(S)).*xi.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L2(S)./xi)-L1(S).^2.*L2(S).*(besseli(1,L2(S)).*besselk(0,L2(S)./xi)+besselk(1,L2(S)).*besseli(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S)).*besseli(1,L1(S)./xi)+(sl.^2-L1(S).^2+S).*((besselk(0,L2(S)./xi).*besseli(1,L1(S)).*L2(S)+besselk(1,L2(S)).*besseli(0,L1(S)./xi).*L1(S)).*besseli(1,L2(S)./xi)+besselk(1,L2(S)./xi).*(besseli(0,L2(S)./xi).*besseli(1,L1(S)).*L2(S)-besseli(0,L1(S)./xi).*besseli(1,L2(S)).*L1(S))).*L2(S).^2)./(((2.*besselk(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)-2.*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+(sl.^2+S).*(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))-L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besselk(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besseli(1,L2(S)./xi)+2.*(besselk(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+besseli(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+1./2.*(sl.^2+S).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))+L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L2(S)./xi)-((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))-L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besselk(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(0,L2(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))-L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))-besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)./xi)+L1(S).*(xi-1).*(xi+1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besseli(1,L1(S)./xi)+(-(2.*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+2.*besseli(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))+L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*(sl.^2+S)).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L1(S)./xi)+(sl.^2-L1(S).^2+S).*((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))-L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besselk(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)-besselk(0,L2(S)./xi).*L2(S).*(xi-1).*(xi+1).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))+L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)./xi)).*L1(S).*L2(S).^2).*besseli(1,L2(S)./xi)+(-2.*(-besseli(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+besseli(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+1./2.*(sl.^2+S).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))-L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L1(S)./xi)-(sl.^2-L1(S).^2+S).*L1(S).*((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))-L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))-besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)+L2(S).*(xi-1).*(xi+1).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S).*besseli(0,L2(S)./xi)-besselk(0,L1(S)./xi).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))-L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S))).*L2(S).^2).*besselk(1,L2(S)./xi)+besselk(1,L1(S)./xi).*(-L1(S).*(xi-1).*(xi+1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))+L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(0,L2(S)./xi)-(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))-L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)));
DK =@(S)( 2.*(sl.^2-L1(S).^2+S).*((besselk(1,L1(S)).*xi.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L1(S)./xi)-besseli(1,L1(S)).*xi.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L1(S)./xi)+L1(S).*L2(S).^2.*(besseli(1,L1(S)).*besselk(0,L1(S)./xi)+besselk(1,L1(S)).*besseli(0,L1(S)./xi)).*(sl.^2-L1(S).^2+S)).*besseli(1,L2(S)./xi)+(sl.^2-L2(S).^2+S).*L1(S).^2.*((-besselk(0,L1(S)./xi).*besseli(1,L2(S)).*L1(S)-besseli(0,L2(S)./xi).*besselk(1,L1(S)).*L2(S)).*besseli(1,L1(S)./xi)+besselk(1,L1(S)./xi).*(besseli(0,L2(S)./xi).*besseli(1,L1(S)).*L2(S)-besseli(0,L1(S)./xi).*besseli(1,L2(S)).*L1(S)))).*U_x./(((-2.*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+2.*besselk(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(sl.^2+S).*(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besseli(1,L1(S)./xi)-(2.*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+2.*besseli(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*(sl.^2+S)).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L1(S)./xi)+(sl.^2-L1(S).^2+S).*((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L1(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)./xi)-besselk(0,L2(S)./xi).*L2(S).*(xi-1).*(xi+1).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)).*L1(S).*L2(S).^2).*besseli(1,L2(S)./xi)+(2.*(besselk(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+besseli(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+1./2.*(sl.^2+S).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besseli(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L2(S)./xi)-((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)+L1(S).*(xi-1).*(xi+1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)-(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besseli(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besseli(1,L1(S)./xi)+(-2.*(-besseli(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+besseli(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+1./2.*(sl.^2+S).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L2(S)./xi)+(-L1(S).*(xi-1).*(xi+1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)-(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besselk(1,L1(S)./xi)-(sl.^2-L1(S).^2+S).*besselk(1,L2(S)./xi).*L1(S).*((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))-L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besseli(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L1(S)./xi)+L2(S).*(xi-1).*(xi+1).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S).*besseli(0,L2(S)./xi)-besselk(0,L1(S)./xi).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*L2(S).^2));
EK =@(S)( 2.*((besselk(1,L2(S)).*xi.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L2(S)./xi)-besseli(1,L2(S)).*xi.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L2(S)./xi)-L1(S).^2.*L2(S).*(besseli(1,L2(S)).*besselk(0,L2(S)./xi)+besselk(1,L2(S)).*besseli(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S)).*besselk(1,L1(S)./xi)+(sl.^2-L1(S).^2+S).*((besselk(0,L2(S)./xi).*besselk(1,L1(S)).*L2(S)-besselk(1,L2(S)).*besselk(0,L1(S)./xi).*L1(S)).*besseli(1,L2(S)./xi)+besselk(1,L2(S)./xi).*(besselk(0,L1(S)./xi).*besseli(1,L2(S)).*L1(S)+besseli(0,L2(S)./xi).*besselk(1,L1(S)).*L2(S))).*L2(S).^2).*(sl.^2-L2(S).^2+S).*U_x./((-(2.*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+2.*besseli(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))+L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*(sl.^2+S)).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besseli(1,L2(S)./xi)-2.*(-besseli(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+besseli(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+1./2.*(sl.^2+S).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))-L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besseli(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L2(S)./xi)+((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))+L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(0,L2(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))+L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))-besseli(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)./xi)-L1(S).*(xi-1).*(xi+1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besselk(1,L1(S)./xi)+((2.*besselk(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)-2.*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+(sl.^2+S).*(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))-L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besselk(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besseli(1,L1(S)./xi)+(sl.^2-L1(S).^2+S).*((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))-L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besselk(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)-besselk(0,L2(S)./xi).*L2(S).*(xi-1).*(xi+1).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))+L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)./xi)).*L1(S).*L2(S).^2).*besseli(1,L2(S)./xi)+(2.*(besselk(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+besseli(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+1./2.*(sl.^2+S).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))+L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besseli(1,L1(S)./xi)-(sl.^2-L1(S).^2+S).*L1(S).*((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L1(S))+L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))-besseli(0,L2(S)).*besseli(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L1(S)./xi)+L2(S).*(xi-1).*(xi+1).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S).*besseli(0,L2(S)./xi)-(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))+L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(0,L1(S)./xi)).*L2(S).^2).*besselk(1,L2(S)./xi)-(L1(S).*(xi-1).*(xi+1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))-L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(1,L2(S))+besselk(0,L2(S)).*besselk(1,L1(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(0,L2(S)./xi)-(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besselk(1,L1(S))+L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(1,L2(S))+besselk(1,L1(S)).*L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S).*besseli(1,L1(S)./xi)));
FK =@(S)( -2.*(sl.^2-L1(S).^2+S).*U_x.*((besselk(1,L1(S)).*xi.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besseli(1,L1(S)./xi)-besseli(1,L1(S)).*xi.*(L1(S)-L2(S)).*(L1(S)+L2(S)).*(sl.^2+S).*besselk(1,L1(S)./xi)+L1(S).*L2(S).^2.*(besseli(1,L1(S)).*besselk(0,L1(S)./xi)+besselk(1,L1(S)).*besseli(0,L1(S)./xi)).*(sl.^2-L1(S).^2+S)).*besselk(1,L2(S)./xi)-((-besselk(0,L2(S)./xi).*besselk(1,L1(S)).*L2(S)+besselk(1,L2(S)).*besselk(0,L1(S)./xi).*L1(S)).*besseli(1,L1(S)./xi)+besselk(1,L1(S)./xi).*(besselk(0,L2(S)./xi).*besseli(1,L1(S)).*L2(S)+besselk(1,L2(S)).*besseli(0,L1(S)./xi).*L1(S))).*(sl.^2-L2(S).^2+S).*L1(S).^2)./((2.*(besselk(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+besseli(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+1./2.*(sl.^2+S).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))+besseli(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besseli(1,L1(S)./xi)-2.*(-besseli(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+besseli(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L2(S)./xi)+1./2.*(sl.^2+S).*(((2.*xi.^2.*L1(S).^2-2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))+L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))-besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besselk(1,L1(S)./xi)-(sl.^2-L1(S).^2+S).*L1(S).*((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))-L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besseli(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L1(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))-L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L1(S)./xi)+L2(S).*(xi-1).*(xi+1).*(besseli(1,L1(S)).*besselk(0,L1(S))+besselk(1,L1(S)).*besseli(0,L1(S))).*(sl.^2-L2(S).^2+S).*besseli(0,L2(S)./xi)).*L2(S).^2).*besselk(1,L2(S)./xi)+((-2.*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+2.*besselk(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(sl.^2+S).*(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S))).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besseli(1,L2(S)./xi)-(L1(S).*(xi-1).*(xi+1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besselk(0,L1(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))-L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besseli(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besseli(1,L1(S)./xi)+(-(2.*besselk(1,L2(S)).*L1(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+2.*besseli(1,L1(S)).*L2(S).*(sl.^2-L2(S).^2+S).*(sl.^2-L1(S).^2+S).*besselk(0,L2(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*(sl.^2+S)).*(L1(S)+L2(S)).*(L1(S)-L2(S)).*xi.*besseli(1,L2(S)./xi)+(-L1(S).*(xi-1).*(xi+1).*(besselk(0,L2(S)).*besseli(1,L2(S))+besselk(1,L2(S)).*besseli(0,L2(S))).*(sl.^2-L1(S).^2+S).*besseli(0,L1(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besseli(1,L2(S))-L2(S).*besseli(0,L2(S)).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besseli(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besselk(0,L2(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,L1(S))+besselk(1,L2(S)).*L1(S).*besseli(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L2(S)./xi)).*(sl.^2-L2(S).^2+S).*L1(S).^2.*L2(S)).*besselk(1,L1(S)./xi)+(sl.^2-L1(S).^2+S).*besseli(1,L2(S)./xi).*((((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besselk(1,L1(S))-besselk(1,L2(S)).*L1(S).*besselk(0,L1(S)).*(xi-1).*(xi+1).*(sl.^2-L2(S).^2+S)).*besseli(0,L1(S)./xi)+(((-2.*xi.^2.*L1(S).^2+2.*xi.^2.*L2(S).^2).*besselk(1,L2(S))+besselk(0,L2(S)).*L2(S).*(xi-1).*(xi+1).*(sl.^2-L1(S).^2+S)).*besseli(1,...

Wow! And nary a comment in sight. Debugging that will be a nightmare, for sure.