function[I41Upper, I41ttUpper, I41xxUpper, I41xtUpper, I41txUpper, I41x1Upper, I41x2Upper, I41t1Upper, I41t2Upper]=JDataTable41Upper(xstar1, xstar2, t1, t2, mR, mL, kappa)
I41Upper=@(omega, t) (-2).*(mR.^2+(-2).*omega.*(omega+((-1).*mR.^2+omega.^2).^(1/2))) ...
  .^(-1).*pi.^(-1).*(omega.*sin(omega.*(t1+(-1).*t2)).*((-1).*((-1) ...
  .*mR.^2+omega.^2).^(-1/2).*cos((mL.^2+omega.^2).^(1/2).*xstar2).* ...
  sin(((-1).*mR.^2+omega.^2).^(1/2).*xstar1)+(mL.^2+omega.^2).^( ...
  -1/2).*cos(((-1).*mR.^2+omega.^2).^(1/2).*xstar1).*sin((mL.^2+ ...
  omega.^2).^(1/2).*xstar2))+cos(omega.*(t1+(-1).*t2)).*(cos(((-1).* ...
  mR.^2+omega.^2).^(1/2).*xstar1).*cos((mL.^2+omega.^2).^(1/2).* ...
  xstar2)+omega.^2.*((mL.^2+omega.^2).*((-1).*mR.^2+omega.^2)).^( ...
  -1/2).*sin(((-1).*mR.^2+omega.^2).^(1/2).*xstar1).*sin((mL.^2+ ...
  omega.^2).^(1/2).*xstar2)));

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

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

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

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

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

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

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

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