function[I41LRLower, I41LRttLower, I41LRxxLower, I41LRxtLower, I41LRtxLower, I41LRx1Lower, I41LRx2Lower, I41LRt1Lower, I41LRt2Lower]=IDataTable41LRULower (xstar1, xstar2, t1, t2, mR, mL, kappa)
I41LRLower=@(omega, t) 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).*( ...
  (-1).*((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)+(-1).*omega.*(omega.*sin( ...
  omega.*(t1+(-1).*t2))+((mR+(-1).*omega).*(mR+omega)).^(1/2).*sin( ...
  omega.*t1).*sin(omega.*t2)).*sin((mL.^2+omega.^2).^(1/2).*xstar2)+ ...
  cos(omega.*t1).*(((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))));

I41LRttLower=@(omega, t) 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).*((-1).*((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)+(-1).* ...
  omega.*(omega.*sin(omega.*(t1+(-1).*t2))+((mR+(-1).*omega).*(mR+ ...
  omega)).^(1/2).*sin(omega.*t1).*sin(omega.*t2)).*sin((mL.^2+ ...
  omega.^2).^(1/2).*xstar2)+cos(omega.*t1).*(((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))) ...
  );

I41LRxxLower=@(omega, t) exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).*mR.^(-2).*( ...
  (mR+(-1).*omega).*(mR+omega)).^(1/2).*pi.^(-1).*csch(kappa.^(-1).* ...
  omega.*pi).*(omega.*cos((mL.^2+omega.^2).^(1/2).*xstar2).*(omega.* ...
  sin(omega.*(t1+(-1).*t2))+((mR+(-1).*omega).*(mR+omega)).^(1/2).* ...
  sin(omega.*t1).*sin(omega.*t2))+sin(omega.*t1).*((-1).*((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).*(((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))));

I41LRxtLower=@(omega, t) (-1).*exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).*mR.^( ...
  -2).*omega.*((mR+(-1).*omega).*(mR+omega).*(mL.^2+omega.^2).^(-1)) ...
  .^(1/2).*pi.^(-1).*csch(kappa.^(-1).*omega.*pi).*(((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).*((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))+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).*((-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))));

I41LRtxLower=@(omega, t) exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).*mR.^(-2).* ...
  omega.*pi.^(-1).*csch(kappa.^(-1).*omega.*pi).*((-1).*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)+(((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).*((-1).*((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))));

I41LRx1Lower=@(omega, t) exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).*mR.^(-2).*( ...
  (mR+(-1).*omega).*(mR+omega).*(mL.^2+omega.^2).^(-1)).^(1/2).* ...
  pi.^(-1).*csch(kappa.^(-1).*omega.*pi).*(((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)+(-1).*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) ...
  )+((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).*((-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))));

I41LRx2Lower=@(omega, t) (-1).*exp(1).^((-1).*(mR.^2+(-1).*omega.^2).^(1/2).*xstar1).*mR.^( ...
  -2).*pi.^(-1).*csch(kappa.^(-1).*omega.*pi).*(omega.*cos((mL.^2+ ...
  omega.^2).^(1/2).*xstar2).*(omega.*sin(omega.*(t1+(-1).*t2))+((mR+ ...
  (-1).*omega).*(mR+omega)).^(1/2).*sin(omega.*t1).*sin(omega.*t2))+ ...
  sin(omega.*t1).*((-1).*((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).*(((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))));

I41LRt1Lower=@(omega, t) (-1).*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).*(((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).*( ...
  (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))+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).*((-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))));

I41LRt2Lower=@(omega, t) 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).*(((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).*( ...
  (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))+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).*((-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))));
