Guest User

Untitled

a guest
May 24th, 2018
103
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.11 KB | None | 0 0
  1. N@Integrate[Sqrt[16 E^(4 I*t) - 2 E^(I*t)]*2 I*E^(I*t), {t, 0, 2 Pi}]
  2. 6.79746 - 6.40204 I
  3.  
  4. NIntegrate[Sqrt[16 E^(4 I*t) - 2 E^(I*t)]*2 I*E^(I*t), {t, 0, 2 Pi}]
  5. 7.60503*10^-15 + 1.37562 I
  6.  
  7. f[z_] = Sqrt[z^4 - z]
  8.  
  9. integrand[t_] =
  10. f[2 E^(I t)]*2 I*E^(I*t) //
  11. ComplexExpand[#, TargetFunctions -> {Re, Im}] & //
  12. FullSimplify[#, 0 < t < 2 Pi] &
  13.  
  14. (* 2 I Sqrt[2] E^(I t) Sqrt[E^(I t) (-1 + 8 E^(3 I t))]
  15.  
  16. 2 I Sqrt[2] (65 - 16 Cos[3 t])^(1/4) (Cos[
  17. t + 1/2 ArcTan[-Cos[t] + 8 Cos[4 t], -Sin[t] + 8 Sin[4 t]]] +
  18. I Sin[t + 1/2 ArcTan[-Cos[t] + 8 Cos[4 t], -Sin[t] + 8 Sin[4 t]]]) *)
  19.  
  20. {sol = Solve[-Sin[t] + 8 Sin[4 t] == 0 && -Cos[t] + 8 Cos[4 t] < 0 &&
  21. 0 < t < 2 Pi, t], t /. sol // N}
  22.  
  23. (* {{{t -> 2 [Pi] -
  24. 2 ArcTan[Sqrt[
  25. Root[-31 + 227 #1 - 221 #1^2 + 33 #1^3 &, 1]]]}, {t ->
  26. 2 ArcTan[Sqrt[
  27. Root[-31 + 227 #1 - 221 #1^2 + 33 #1^3 &, 1]]]}, {t ->
  28. 2 [Pi] -
  29. 2 ArcTan[Sqrt[
  30. Root[-31 + 227 #1 - 221 #1^2 + 33 #1^3 &, 3]]]}, {t ->
  31. 2 ArcTan[Sqrt[
  32. Root[-31 + 227 #1 - 221 #1^2 + 33 #1^3 &, 3]]]}}, {5.51943,
  33. 0.763757, 3.94961, 2.33357}} *)
  34.  
  35. integrand2[t_] =
  36. Piecewise[{{integrand[t],
  37. 0 <= t <= (t /. sol[[2]])}, {-integrand[t], (t /. sol[[2]]) <
  38. t <= (t /. sol[[4]])}, {integrand[t], (t /. sol[[4]]) <
  39. t <= (t /. sol[[3]])}, {-integrand[t], (t /. sol[[3]]) <
  40. t <= (t /. sol[[1]])}, {integrand[t], (t /. sol[[1]]) < t <=
  41. 2 Pi}}, 0]
  42.  
  43. NIntegrate[integrand2[t], {t, 0, 2 Pi}]
  44.  
  45. (* 1.77636*10^-15 - 3.14159 I *)
  46.  
  47. l1 = Int[integrand2[t][[1, All, 1]], t];
  48.  
  49. l2 = integrand2[t][[1, All, 2]];
  50.  
  51. tp = Transpose[{l1, l2}];
  52.  
  53. rint[t_] = Piecewise[tp, 0]
  54.  
  55. Plot[{rint[t] // Re, rint[t] // Im}, {t, 0, 2 Pi}, PlotPoints -> 200,
  56. PlotStyle -> {Blue, Red}, GridLines -> {{2/3 Pi, 4/3 Pi}, Automatic},
  57. Ticks -> {{2/3 Pi, 4/3 Pi}, Automatic}]
  58.  
  59. lim = Plus @@ {Limit[rint[t], t -> 2 Pi,
  60. Direction -> 1], -Limit[rint[t], t -> 4/3 Pi, Direction -> -1],
  61. Limit[rint[t], t -> 4/3 Pi,
  62. Direction -> 1], -Limit[rint[t], t -> 2/3 Pi, Direction -> -1],
  63. Limit[rint[t], t -> 2/3 Pi,
  64. Direction -> 1], -Limit[rint[t], t -> 0, Direction -> -1]} //
  65. FullSimplify
  66.  
  67. (* - I Pi *)
Add Comment
Please, Sign In to add comment