Advertisement
Guest User

Untitled

a guest
Jun 20th, 2019
56
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.62 KB | None | 0 0
  1. timewj2 = AbsoluteTiming[
  2. Wjlambda5 = (1/8.952*10^8)*
  3. NIntegrate[
  4. Dot[Re[{0, -8.952*10^8*(15.675459922348419449)*(((-8.
  5. 9652950534407867999303353030387276032035253437793*10^768 -
  6. 2.1141452007443239580927869423477693294460547800052
  7. 7*10^769*I)*
  8. BesselJ[
  9. 1, (75062.4870217452581789128 +
  10. 75062.4983440564559965704*I )*
  11. x] + (2.1141452007443239580927869423477693294*10^
  12. 769 - 8.965295053440786799930335303038727603*10^768*I)*
  13. BesselY[
  14. 1, (75062.4870217452581789128 +
  15. 75062.4983440564559965704*I )*x])*
  16. Exp[I*(1)*(41.2281675906100480252)*s]), 0}],
  17. Re[{0, -8.952*10^8*(15.675459922348419449)*(((-8.
  18. 9652950534407867999303353030387276032035253437793*10^768 -
  19. 2.1141452007443239580927869423477693294460547800052
  20. 7*10^769*I)*
  21. BesselJ[
  22. 1, (75062.4870217452581789128 +
  23. 75062.4983440564559965704*I )*
  24. x] + (2.1141452007443239580927869423477693294*10^
  25. 769 - 8.965295053440786799930335303038727603*10^768*I)*
  26. BesselY[
  27. 1, (75062.4870217452581789128 +
  28. 75062.4983440564559965704*I )*x])*
  29. Exp[I*(1)*(41.2281675906100480252)*s]), 0}]]*
  30. x, {x, (0.02357124714249428777766010023597686995),
  31. (0.04022197675031113473009371147375231401)}, {s,
  32. 0, (0.15240030480060962059241091992589645088)},
  33. WorkingPrecision -> (2000 - 1), MaxPoints -> 80];
  34. ];
  35. Print["Time Wj5 (s) = ", N[timewj2[[1]],10]];
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement