function[I21Upper, I21ttUpper, I21xxUpper, I21xtUpper, I21txUpper, I21x1Upper, I21x2Upper, I21t1Upper, I21t2Upper]=IDataTable21Upper(x1, x2, t1, t2, mR, mL, kappa)
I21Upper=@(omega,t)(1/2).*omega.^(-1).*(mL.^2+omega.^2).^(-1/2).*(mR.^2+(-2).*omega.* ...
  (omega+((-1).*mR.^2+omega.^2).^(1/2))).^(-2).*pi.^(-1).*coth( ...
  kappa.^(-1).*omega.*pi).*((mR.^4.*(mL.^2+omega.^2).^(1/2)+8.* ...
  omega.^3.*(omega.*(mL.^2+omega.^2).^(1/2)+((mL.^2+omega.^2).*((-1) ...
  .*mR.^2+omega.^2)).^(1/2))+(-4).*mR.^2.*omega.*(2.*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.^4.*(mL.^2+omega.^2).^(1/2).*cos(omega.*(t1+( ...
  -1).*t2+xstar1)).*cos((mL.^2+omega.^2).^(1/2).*xstar2)+omega.*(( ...
  -1).*(mR.^4+8.*omega.^3.*(omega+((-1).*mR.^2+omega.^2).^(1/2))+( ...
  -4).*mR.^2.*omega.*(2.*omega+((-1).*mR.^2+omega.^2).^(1/2))).*sin( ...
  omega.*(t1+(-1).*t2+(-1).*xstar1))+mR.^4.*sin(omega.*(t1+(-1).*t2+ ...
  xstar1))).*sin((mL.^2+omega.^2).^(1/2).*xstar2));

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

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

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

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

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

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

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

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