Advertisement
Guest User

Untitled

a guest
Aug 4th, 2015
147
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.76 KB | None | 0 0
  1. $Version
  2. (* 10.1.0 for Microsoft Windows (64-bit) (March 24, 2015) *)
  3.  
  4. Integrate[1/(Sqrt[x] Sqrt[1 - x + x^2]), {x, 1, 2}]
  5. (* Integrate[1/(Sqrt[x] Sqrt[1 - x + x^2]), {x, 1, 2}] *)
  6.  
  7. Integrate[1/(Sqrt[x] Sqrt[1 - x + x^2]), x]
  8. (* (2 (-1)^(1/6) Sqrt[1 - (-1)^(1/3)/x] Sqrt[1 + (-1)^(2/3)/x] x
  9. EllipticF[I ArcSinh[(-1)^(1/3)/Sqrt[x]], (-1)^(2/3)])/Sqrt[1 - x + x^2] *)
  10.  
  11. FullSimplify[(% /. x -> 2) - (% /. x -> 1)]
  12. (* 2 (-1)^(1/6) (-EllipticF[I ArcSinh[(-1)^(1/3)], (-1)^(2/3)] +
  13. EllipticF[I ArcSinh[(-1)^(1/3)/Sqrt[2]], (-1)^(2/3)]) *)
  14.  
  15. % - NIntegrate[1/(Sqrt[x] Sqrt[1 - x + x^2]), {x, 1, 2}, WorkingPrecision -> 30]
  16. (* 0. 10^-31 + 0. 10^-47 I *)
  17.  
  18. $Version
  19.  
  20. (* Out[1]= "8.0 for Microsoft Windows (64-bit) (October 7, 2011)" *)
  21.  
  22. f[x_] = 1/(Sqrt[x] Sqrt[1 - x + x^2]);
  23.  
  24. Timing[Integrate[1/(Sqrt[x] Sqrt[1 - x + x^2]), {x, 1, 2}] ]
  25.  
  26. (*
  27. Out[3]= {3.011, 2 (-1)^(
  28. 1/6) (-EllipticF[I ArcSinh[(-1)^(1/3)], (-1)^(2/3)] +
  29. EllipticF[I ArcSinh[(-1)^(1/3)/Sqrt[2]], (-1)^(2/3)])}
  30. *)
  31.  
  32. % // N
  33.  
  34. (* Out[4]= {3.011, 0.646172 - 5.55112*10^-17 I} *)
  35.  
  36. NIntegrate[f[x], {x, 1, 2}]
  37.  
  38. (*
  39. Out[5]= 0.646172
  40. *)
  41.  
  42. Integrate[1/(Sqrt[x] Sqrt[1 - x + x^2]), x]
  43.  
  44. (*
  45. Out[6]= (2 (-1)^(1/6) Sqrt[1 - (-1)^(1/3)/x] Sqrt[
  46. 1 + (-1)^(2/3)/x] x EllipticF[I ArcSinh[(-1)^(1/3)/Sqrt[x]], (-1)^(
  47. 2/3)])/Sqrt[1 - x + x^2]
  48. *)
  49.  
  50. $Version
  51.  
  52. (* Out[2]=5.2 for Microsoft Windows x86 (64 bit) (June 20, 2005) *)
  53.  
  54. Timing[Integrate[1/(Sqrt[x] Sqrt[1 - x + x^2]), {x, 1, 2}] ]
  55.  
  56. (*
  57. Out[4]=
  58. {0.2030*Second, 2*(-1)^(1/6)*(-EllipticF[I*ArcSinh[(-1)^(1/3)], (-1)^(2/3)] +
  59. EllipticF[I*ArcSinh[(-1)^(1/3)/Sqrt[2]], (-1)^(2/3)])}
  60. *)
  61.  
  62. N[InverseJacobiCN[-1/3, 3/4] - EllipticK[3/4], 20]
  63. 0.64617199330515618196
  64.  
  65. NIntegrate[1/Sqrt[x (1 - x + x^2)], {x, 1, 2}, WorkingPrecision -> 20]
  66. 0.64617199330515618293
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement