function[I31LRUpper, I31ttLRUpper, I31xxLRUpper, I31xtLRUpper, I31txLRUpper, I31x1LRUpper, I31x2LRUpper, I31t1LRUpper, I31t2LRUpper]=IDataTable31LRUpper(x1, x2, t1, t2, mR, mL, kappa)
I31LRUpper=@(omega,t)omega.^(-1).*(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).*(sin(omega.*(t1+(-1).*t2)).*(((-1).* ...
  mR.^2.*(mL.^2+omega.^2).^(1/2)+omega.*(omega.*(mL.^2+omega.^2).^( ...
  1/2)+((mL.^2+omega.^2).*((-1).*mR.^2+omega.^2)).^(1/2))).*cos(( ...
  mL.^2+omega.^2).^(1/2).*xstar2).*sin(omega.*xstar1)+(-1).* ...
  omega.^2.*(omega+((-1).*mR.^2+omega.^2).^(1/2)).*cos(omega.* ...
  xstar1).*sin((mL.^2+omega.^2).^(1/2).*xstar2))+omega.*cos(omega.*( ...
  t1+(-1).*t2)).*((omega.*(mL.^2+omega.^2).^(1/2)+((mL.^2+omega.^2) ...
  .*((-1).*mR.^2+omega.^2)).^(1/2)).*cos(omega.*xstar1).*cos((mL.^2+ ...
  omega.^2).^(1/2).*xstar2)+((-1).*mR.^2+omega.*(omega+((-1).*mR.^2+ ...
  omega.^2).^(1/2))).*sin(omega.*xstar1).*sin((mL.^2+omega.^2).^( ...
  1/2).*xstar2)));

I31LRttUpper=@(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).*(((-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)+(-1).*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));

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

I31LRxtUpper=@(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).*(((-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)+(-1).*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));

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

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

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

I31LRt1Upper=@(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).*csch(kappa.^(-1).* ...
  omega.*pi).*(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));

I31LRt2Upper=@(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).*csch(kappa.^(-1).* ...
  omega.*pi).*(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));
