function[I41Lower, I41ttLower, I41xxLower, I41xtLower, I41txLower, I41x1Lower, I41x2Lower, I41t1Lower, I41t2Lower]=IDataTable41ULower (xstar1, xstar2, t1, t2, mR, mL, kappa)
I41Lower=@(omega, t) (1/2).*exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).* ...
  mR.^(-2).*(mL.^2+omega.^2).^(-1/2).*pi.^(-1).*csch(kappa.^(-1).* ...
  omega.*pi).^2.*(((mR+(-1).*omega).*(mR+omega).*(mL.^2+omega.^2)) ...
  .^(1/2).*cos(omega.*t2).*cos((mL.^2+omega.^2).^(1/2).*xstar2).* ...
  sin(omega.*t1)+omega.*(mL.^2+omega.^2).^(1/2).*cos((mL.^2+ ...
  omega.^2).^(1/2).*xstar2).*sin(omega.*t1).*sin(omega.*t2)+omega.*( ...
  omega.*sin(omega.*(t1+(-1).*t2))+(-1).*((mR+(-1).*omega).*(mR+ ...
  omega)).^(1/2).*sin(omega.*t1).*sin(omega.*t2)).*sin((mL.^2+ ...
  omega.^2).^(1/2).*xstar2)+cos(omega.*t1).*((-1).*((mR+(-1).*omega) ...
  .*(mR+omega).*(mL.^2+omega.^2)).^(1/2).*cos((mL.^2+omega.^2).^( ...
  1/2).*xstar2).*sin(omega.*t2)+omega.*cos(omega.*t2).*((mL.^2+ ...
  omega.^2).^(1/2).*cos((mL.^2+omega.^2).^(1/2).*xstar2)+(-1).*((mR+ ...
  (-1).*omega).*(mR+omega)).^(1/2).*sin((mL.^2+omega.^2).^(1/2).* ...
  xstar2)))).*sinh(2.*kappa.^(-1).*omega.*pi);

I41ttLower=@(omega, t) (1/2).*exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).* ...
  mR.^(-2).*omega.^2.*(mL.^2+omega.^2).^(-1/2).*pi.^(-1).*csch( ...
  kappa.^(-1).*omega.*pi).^2.*(((mR+(-1).*omega).*(mR+omega).*( ...
  mL.^2+omega.^2)).^(1/2).*cos(omega.*t2).*cos((mL.^2+omega.^2).^( ...
  1/2).*xstar2).*sin(omega.*t1)+omega.*(mL.^2+omega.^2).^(1/2).*cos( ...
  (mL.^2+omega.^2).^(1/2).*xstar2).*sin(omega.*t1).*sin(omega.*t2)+ ...
  omega.*(omega.*sin(omega.*(t1+(-1).*t2))+(-1).*((mR+(-1).*omega).* ...
  (mR+omega)).^(1/2).*sin(omega.*t1).*sin(omega.*t2)).*sin((mL.^2+ ...
  omega.^2).^(1/2).*xstar2)+cos(omega.*t1).*((-1).*((mR+(-1).*omega) ...
  .*(mR+omega).*(mL.^2+omega.^2)).^(1/2).*cos((mL.^2+omega.^2).^( ...
  1/2).*xstar2).*sin(omega.*t2)+omega.*cos(omega.*t2).*((mL.^2+ ...
  omega.^2).^(1/2).*cos((mL.^2+omega.^2).^(1/2).*xstar2)+(-1).*((mR+ ...
  (-1).*omega).*(mR+omega)).^(1/2).*sin((mL.^2+omega.^2).^(1/2).* ...
  xstar2)))).*sinh(2.*kappa.^(-1).*omega.*pi);

I41xxLower=@(omega, t) (1/2).*exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).* ...
  mR.^(-2).*pi.^(-1).*csch(kappa.^(-1).*omega.*pi).^2.*((mR+(-1).* ...
  omega).*omega.*(mR+omega).*cos(omega.*(t1+(-1).*t2)).*cos((mL.^2+ ...
  omega.^2).^(1/2).*xstar2)+cos(omega.*t2).*((-1).*omega.^2.*((mR+( ...
  -1).*omega).*(mR+omega)).^(1/2).*cos((mL.^2+omega.^2).^(1/2).* ...
  xstar2).*sin(omega.*t1)+(omega.*((mR+(-1).*omega).*(mR+omega).*( ...
  mL.^2+omega.^2)).^(1/2).*cos(omega.*t1)+(mR+(-1).*omega).*(mR+ ...
  omega).*(mL.^2+omega.^2).^(1/2).*sin(omega.*t1)).*sin((mL.^2+ ...
  omega.^2).^(1/2).*xstar2))+sin(omega.*t2).*(omega.*((mR+(-1).* ...
  omega).*(mR+omega).*(mL.^2+omega.^2)).^(1/2).*sin(omega.*t1).*sin( ...
  (mL.^2+omega.^2).^(1/2).*xstar2)+cos(omega.*t1).*(omega.^2.*((mR+( ...
  -1).*omega).*(mR+omega)).^(1/2).*cos((mL.^2+omega.^2).^(1/2).* ...
  xstar2)+(mL.^2+omega.^2).^(1/2).*((-1).*mR.^2+omega.^2).*sin(( ...
  mL.^2+omega.^2).^(1/2).*xstar2)))).*sinh(2.*kappa.^(-1).*omega.* ...
  pi);

I41xtLower=@(omega, t) (1/2).*exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).* ...
  mR.^(-2).*omega.*(mL.^2+omega.^2).^(-1/2).*pi.^(-1).*csch(kappa.^( ...
  -1).*omega.*pi).^2.*(cos(omega.*t1).*(omega.*sin(omega.*t2).*((( ...
  mR+(-1).*omega).*(mR+omega).*(mL.^2+omega.^2)).^(1/2).*cos((mL.^2+ ...
  omega.^2).^(1/2).*xstar2)+((-1).*mR+omega).*(mR+omega).*sin(( ...
  mL.^2+omega.^2).^(1/2).*xstar2))+cos(omega.*t2).*((mR+(-1).*omega) ...
  .*(mR+omega).*(mL.^2+omega.^2).^(1/2).*cos((mL.^2+omega.^2).^(1/2) ...
  .*xstar2)+omega.^2.*((mR+(-1).*omega).*(mR+omega)).^(1/2).*sin(( ...
  mL.^2+omega.^2).^(1/2).*xstar2)))+sin(omega.*t1).*(omega.*cos( ...
  omega.*t2).*((-1).*((mR+(-1).*omega).*(mR+omega).*(mL.^2+omega.^2) ...
  ).^(1/2).*cos((mL.^2+omega.^2).^(1/2).*xstar2)+(mR+(-1).*omega).*( ...
  mR+omega).*sin((mL.^2+omega.^2).^(1/2).*xstar2))+sin(omega.*t2).*( ...
  (mR+(-1).*omega).*(mR+omega).*(mL.^2+omega.^2).^(1/2).*cos((mL.^2+ ...
  omega.^2).^(1/2).*xstar2)+omega.^2.*((mR+(-1).*omega).*(mR+omega)) ...
  .^(1/2).*sin((mL.^2+omega.^2).^(1/2).*xstar2)))).*sinh(2.*kappa.^( ...
  -1).*omega.*pi);

I41txLower=@(omega, t) (1/2).*exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).* ...
  mR.^(-2).*omega.*pi.^(-1).*csch(kappa.^(-1).*omega.*pi).^2.*( ...
  omega.^2.*cos(omega.*(t1+(-1).*t2)).*cos((mL.^2+omega.^2).^(1/2).* ...
  xstar2)+cos(omega.*t2).*(omega.*((mR+(-1).*omega).*(mR+omega)).^( ...
  1/2).*cos((mL.^2+omega.^2).^(1/2).*xstar2).*sin(omega.*t1)+((-1).* ...
  ((mR+(-1).*omega).*(mR+omega).*(mL.^2+omega.^2)).^(1/2).*cos( ...
  omega.*t1)+omega.*(mL.^2+omega.^2).^(1/2).*sin(omega.*t1)).*sin(( ...
  mL.^2+omega.^2).^(1/2).*xstar2))+(-1).*sin(omega.*t2).*(((mR+(-1) ...
  .*omega).*(mR+omega).*(mL.^2+omega.^2)).^(1/2).*sin(omega.*t1).* ...
  sin((mL.^2+omega.^2).^(1/2).*xstar2)+omega.*cos(omega.*t1).*(((mR+ ...
  (-1).*omega).*(mR+omega)).^(1/2).*cos((mL.^2+omega.^2).^(1/2).* ...
  xstar2)+(mL.^2+omega.^2).^(1/2).*sin((mL.^2+omega.^2).^(1/2).* ...
  xstar2)))).*sinh(2.*kappa.^(-1).*omega.*pi);

I41x1Lower=@(omega, t) (-1/4).*exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).* ...
  mR.^(-2).*(mL.^2+omega.^2).^(-1/2).*pi.^(-1).*csch(kappa.^(-1).* ...
  omega.*pi).^2.*(omega.*(omega.*((mR+(-1).*omega).*(mR+omega)).^( ...
  1/2)+((mR+(-1).*omega).*(mR+omega).*(mL.^2+omega.^2)).^(1/2)).* ...
  cos(omega.*t1+(-1).*omega.*t2+(-1).*(mL.^2+omega.^2).^(1/2).* ...
  xstar2)+omega.*((-1).*omega.*((mR+(-1).*omega).*(mR+omega)).^(1/2) ...
  +((mR+(-1).*omega).*(mR+omega).*(mL.^2+omega.^2)).^(1/2)).*cos( ...
  omega.*t1+(-1).*omega.*t2+(mL.^2+omega.^2).^(1/2).*xstar2)+(-2).*( ...
  mR+(-1).*omega).*(mR+omega).*((-1).*(mL.^2+omega.^2).^(1/2).*cos(( ...
  mL.^2+omega.^2).^(1/2).*xstar2).*sin(omega.*(t1+(-1).*t2))+omega.* ...
  cos(omega.*(t1+(-1).*t2)).*sin((mL.^2+omega.^2).^(1/2).*xstar2))) ...
  .*sinh(2.*kappa.^(-1).*omega.*pi);

I41x2Lower=@(omega, t) (-1/2).*exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).* ...
  mR.^(-2).*pi.^(-1).*csch(kappa.^(-1).*omega.*pi).^2.*(omega.*cos(( ...
  mL.^2+omega.^2).^(1/2).*xstar2).*((-1).*omega.*sin(omega.*(t1+(-1) ...
  .*t2))+((mR+(-1).*omega).*(mR+omega)).^(1/2).*sin(omega.*t1).*sin( ...
  omega.*t2))+sin(omega.*t1).*(((mR+(-1).*omega).*(mR+omega).*( ...
  mL.^2+omega.^2)).^(1/2).*cos(omega.*t2)+omega.*(mL.^2+omega.^2).^( ...
  1/2).*sin(omega.*t2)).*sin((mL.^2+omega.^2).^(1/2).*xstar2)+cos( ...
  omega.*t1).*((-1).*((mR+(-1).*omega).*(mR+omega).*(mL.^2+omega.^2) ...
  ).^(1/2).*sin(omega.*t2).*sin((mL.^2+omega.^2).^(1/2).*xstar2)+ ...
  omega.*cos(omega.*t2).*(((mR+(-1).*omega).*(mR+omega)).^(1/2).* ...
  cos((mL.^2+omega.^2).^(1/2).*xstar2)+(mL.^2+omega.^2).^(1/2).*sin( ...
  (mL.^2+omega.^2).^(1/2).*xstar2)))).*sinh(2.*kappa.^(-1).*omega.* ...
  pi);

I41t1Lower=@(omega, t) (1/2).*exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).* ...
  mR.^(-2).*omega.*(mL.^2+omega.^2).^(-1/2).*pi.^(-1).*csch(kappa.^( ...
  -1).*omega.*pi).^2.*(cos(omega.*t1).*(cos(omega.*t2).*(((mR.^2+( ...
  -1).*omega.^2).*(mL.^2+omega.^2)).^(1/2).*cos((mL.^2+omega.^2).^( ...
  1/2).*xstar2)+omega.^2.*sin((mL.^2+omega.^2).^(1/2).*xstar2))+ ...
  omega.*sin(omega.*t2).*((mL.^2+omega.^2).^(1/2).*cos((mL.^2+ ...
  omega.^2).^(1/2).*xstar2)+(-1).*(mR.^2+(-1).*omega.^2).^(1/2).* ...
  sin((mL.^2+omega.^2).^(1/2).*xstar2)))+sin(omega.*t1).*(sin( ...
  omega.*t2).*(((mR.^2+(-1).*omega.^2).*(mL.^2+omega.^2)).^(1/2).* ...
  cos((mL.^2+omega.^2).^(1/2).*xstar2)+omega.^2.*sin((mL.^2+ ...
  omega.^2).^(1/2).*xstar2))+cos(omega.*t2).*((-1).*omega.*(mL.^2+ ...
  omega.^2).^(1/2).*cos((mL.^2+omega.^2).^(1/2).*xstar2)+omega.*( ...
  mR.^2+(-1).*omega.^2).^(1/2).*sin((mL.^2+omega.^2).^(1/2).*xstar2) ...
  ))).*sinh(2.*kappa.^(-1).*omega.*pi);

I41t2Lower=@(omega, t) (-1/2).*exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).* ...
  mR.^(-2).*omega.*(mL.^2+omega.^2).^(-1/2).*pi.^(-1).*csch(kappa.^( ...
  -1).*omega.*pi).^2.*(((mR+(-1).*omega).*(mR+omega).*(mL.^2+ ...
  omega.^2)).^(1/2).*cos((mL.^2+omega.^2).^(1/2).*xstar2).*sin( ...
  omega.*t1).*sin(omega.*t2)+omega.^2.*cos(omega.*(t1+(-1).*t2)).* ...
  sin((mL.^2+omega.^2).^(1/2).*xstar2)+omega.*cos(omega.*t2).*sin( ...
  omega.*t1).*((-1).*(mL.^2+omega.^2).^(1/2).*cos((mL.^2+omega.^2) ...
  .^(1/2).*xstar2)+((mR+(-1).*omega).*(mR+omega)).^(1/2).*sin(( ...
  mL.^2+omega.^2).^(1/2).*xstar2))+cos(omega.*t1).*(((mR+(-1).* ...
  omega).*(mR+omega).*(mL.^2+omega.^2)).^(1/2).*cos(omega.*t2).*cos( ...
  (mL.^2+omega.^2).^(1/2).*xstar2)+omega.*sin(omega.*t2).*((mL.^2+ ...
  omega.^2).^(1/2).*cos((mL.^2+omega.^2).^(1/2).*xstar2)+(-1).*((mR+ ...
  (-1).*omega).*(mR+omega)).^(1/2).*sin((mL.^2+omega.^2).^(1/2).* ...
  xstar2)))).*sinh(2.*kappa.^(-1).*omega.*pi);
