function[I41LRUpper, I41LRttUpper, I41LRxxUpper, I41LRxtUpper, I41LRtxUpper, I41LRx1Upper, I41LRx2Upper, I41LRt1Upper, I41LRt2Upper]=IDataTable41LRUpper(xstar1, xstar2, t1, t2, mR, mL, kappa)
I41LRUpper=@(omega, t) (1/2).*(mL.^2+omega.^2).^(-1/2).*((-1).*mR.^2+2.*omega.*(omega+(( ...
  -1).*mR.^2+omega.^2).^(1/2))).^(-1).*pi.^(-1).*((((mL.^2+omega.^2) ...
  .*((-1).*mR.^2+omega.^2)).^(1/2)+omega.*((-1).*omega+(mL.^2+ ...
  omega.^2).^(1/2)+(-1).*((-1).*mR.^2+omega.^2).^(1/2))).*cos( ...
  omega.*t1+(-1).*omega.*t2+(-1).*((-1).*mR.^2+omega.^2).^(1/2).* ...
  xstar1+(-1).*(mL.^2+omega.^2).^(1/2).*xstar2)+(((mL.^2+omega.^2).* ...
  ((-1).*mR.^2+omega.^2)).^(1/2)+omega.*(omega+(mL.^2+omega.^2).^( ...
  1/2)+((-1).*mR.^2+omega.^2).^(1/2))).*cos(omega.*t1+(-1).*omega.* ...
  t2+(-1).*((-1).*mR.^2+omega.^2).^(1/2).*xstar1+(mL.^2+omega.^2).^( ...
  1/2).*xstar2)).*csch(kappa.^(-1).*omega.*pi);

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

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

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

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

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

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

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

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