function[I31Upper, I31ttUpper, I31xxUpper, I31xtUpper, I31txUpper, I31x1Upper, I31x2Upper, I31t1Upper, I31t2Upper]=IDataTable31Upper(x1, x2, t1, t2, mR, mL, kappa)
I31Upper=@(omega,t)(1/4).*mR.^2.*omega.^(-1).*(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.*(mL.^2+omega.^2) ...
  .^(1/2)+2.*omega.*(omega.*(mL.^2+omega.^2).^(1/2)+((mL.^2+ ...
  omega.^2).*((-1).*mR.^2+omega.^2)).^(1/2))).*cos(omega.*(t1+(-1).* ...
  t2+(-1).*xstar1)).*cos((mL.^2+omega.^2).^(1/2).*xstar2)+mR.^2.*( ...
  mL.^2+omega.^2).^(1/2).*cos(omega.*(t1+(-1).*t2+xstar1)).*cos(( ...
  mL.^2+omega.^2).^(1/2).*xstar2)+omega.*(((-1).*mR.^2+2.*omega.*( ...
  omega+((-1).*mR.^2+omega.^2).^(1/2))).*sin(omega.*(t1+(-1).*t2+( ...
  -1).*xstar1))+mR.^2.*sin(omega.*(t1+(-1).*t2+xstar1))).*sin(( ...
  mL.^2+omega.^2).^(1/2).*xstar2)).*sinh(2.*kappa.^(-1).*omega.*pi); ...
  

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

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

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

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

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

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

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

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