function[I41Upper, I41ttUpper, I41xxUpper, I41xtUpper, I41txUpper, I41x1Upper, I41x2Upper, I41t1Upper, I41t2Upper]=IDataTable41Upper(xstar1, xstar2, t1, t2, mR, mL, kappa)
I41Upper=@(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)+(-1).*omega.*(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)+(omega.^2+((mL.^2+ ...
  omega.^2).*((-1).*mR.^2+omega.^2)).^(1/2)+(-1).*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)).*coth(kappa.^(-1).*omega.*pi);

I41ttUpper=@(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).*(((-1).*(( ...
  mL.^2+omega.^2).*((-1).*mR.^2+omega.^2)).^(1/2)+omega.*(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)+(-1).*(omega.^2+(( ...
  mL.^2+omega.^2).*((-1).*mR.^2+omega.^2)).^(1/2)+(-1).*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);

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

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

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

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

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

I41t1Upper=@(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).*coth( ...
  kappa.^(-1).*omega.*pi).*((((mL.^2+omega.^2).*((-1).*mR.^2+ ...
  omega.^2)).^(1/2)+(-1).*omega.*(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)+(omega.^2+((mL.^2+omega.^2).*((-1).*mR.^2+ ...
  omega.^2)).^(1/2)+(-1).*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));

I41t2Upper=@(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).*coth( ...
  kappa.^(-1).*omega.*pi).*(((-1).*((mL.^2+omega.^2).*((-1).*mR.^2+ ...
  omega.^2)).^(1/2)+omega.*(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)+(-1).*(omega.^2+((mL.^2+omega.^2).*((-1).*mR.^2+ ...
  omega.^2)).^(1/2)+(-1).*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));
