function[I42Upper, I42ttUpper, I42xxUpper, I42xtUpper, I42txUpper, I42x1Upper, I42x2Upper, I42t1Upper, I42t2Upper]=IDataTable42Upper(x1, x2, t1, t2, mR, mL, kappa)
I42Upper=@(omega,t)((-1).*omega+((-1).*mR.^2+omega.^2).^(1/2)).*(mR.^2+(-2).*omega.*( ...
  omega+((-1).*mR.^2+omega.^2).^(1/2))).^(-1).*pi.^(-1).*cos(((-1).* ...
  mR.^2+omega.^2).^(1/2).*xstar1+omega.*((-1).*t1+t2+xstar2)).*coth( ...
  kappa.^(-1).*omega.*pi);

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

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

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

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

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

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

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

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