Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- timewj2 = AbsoluteTiming[
- Wjlambda5 = (1/8.952*10^8)*
- NIntegrate[
- Dot[Re[{0, -8.952*10^8*(15.675459922348419449)*(((-8.
- 9652950534407867999303353030387276032035253437793*10^768 -
- 2.1141452007443239580927869423477693294460547800052
- 7*10^769*I)*
- BesselJ[
- 1, (75062.4870217452581789128 +
- 75062.4983440564559965704*I )*
- x] + (2.1141452007443239580927869423477693294*10^
- 769 - 8.965295053440786799930335303038727603*10^768*I)*
- BesselY[
- 1, (75062.4870217452581789128 +
- 75062.4983440564559965704*I )*x])*
- Exp[I*(1)*(41.2281675906100480252)*s]), 0}],
- Re[{0, -8.952*10^8*(15.675459922348419449)*(((-8.
- 9652950534407867999303353030387276032035253437793*10^768 -
- 2.1141452007443239580927869423477693294460547800052
- 7*10^769*I)*
- BesselJ[
- 1, (75062.4870217452581789128 +
- 75062.4983440564559965704*I )*
- x] + (2.1141452007443239580927869423477693294*10^
- 769 - 8.965295053440786799930335303038727603*10^768*I)*
- BesselY[
- 1, (75062.4870217452581789128 +
- 75062.4983440564559965704*I )*x])*
- Exp[I*(1)*(41.2281675906100480252)*s]), 0}]]*
- x, {x, (0.02357124714249428777766010023597686995),
- (0.04022197675031113473009371147375231401)}, {s,
- 0, (0.15240030480060962059241091992589645088)},
- WorkingPrecision -> (2000 - 1), MaxPoints -> 80];
- ];
- Print["Time Wj5 (s) = ", N[timewj2[[1]],10]];
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement