Guest User

Untitled

a guest
Feb 17th, 2018
74
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.63 KB | None | 0 0
  1. R = ((1 - x/L)*W^2*(16*10^-20)^2*Pi*Sqrt[me]*LnLumbda)/( T^(3/2)*e0);
  2.  
  3. Simplify[1/Sqrt[[Pi]]
  4. E^-S^2/(S - (W + 2 I R - Om)/(
  5. vt (W/CC (1 - (1 - x/L)/(1 + (Wc)/W))^(1/2))))]
  6. -((0.5641895835477563` E^-S^2 Sqrt[
  7. x])/((0.016006834477809786` + 0.00008178287140291831` I) -
  8. 1.` S Sqrt[x] - (0.` + 8.178287140291829` I) x))
  9. A1[x_] :=
  10. NIntegrate[-((0.5641895835477563` E^-S^2 Sqrt[
  11. x])/((0.016006834477809786` + 0.00008178287140291831` I) -
  12. 1.` S Sqrt[
  13. x] - (0.` + 8.178287140291829` I) x)), {S, -Infinity,
  14. Infinity}]
  15. yWKB[x_] :=
  16. b/Sqrt[Abs[[Kappa][x]]] Exp[NIntegrate[[Kappa][xp], {xp}]]
  17. [Kappa][x_] := Sqrt[W^2/CC^2 + (W wp^2)/(CC^2 vt k) A1[x]]
  18. Simplify[[Kappa][x]]
  19. Sqrt[(3.947841760435743`*^13 Sqrt[
  20. x] + (6.319244960388013`*^11 -
  21. 6.319244960388012`*^16 x) NIntegrate[-((
  22. 0.5641895835477563` E^-S^2 Sqrt[
  23. x])/((0.016006834477809786` + 0.00008178287140291831` I) -
  24. 1.` S Sqrt[
  25. x] - (0.` +
  26. 8.178287140291829` I) x)), {S, -[Infinity],
  27. [Infinity]}])/Sqrt[x]]
  28. [Kappa][x_] := [Sqrt](3.947841760435743`*^13 +
  29. 1/Sqrt[x] (6.319244960388013`*^11 -
  30. 6.319244960388012`*^16 x) NIntegrate[-((0.5641895835477563`
  31. E^-S^2 Sqrt[
  32. x])/((0.016006834477809786` + 0.00008178287140291831` I) -
  33. 1.` S Sqrt[
  34. x] - (0.` +
  35. 8.178287140291829` I) x)), {S, -[Infinity],
  36. [Infinity]}])
  37. yWKB[-L]
  38. NIntegrate::ilim
  39. 0.0004027134076113933` b E^NIntegrate[[Kappa][xp], {xp}]
  40. Solve[0.0004027134076113933` b E^NIntegrate[[Kappa][xp], {xp}] ==
  41. 10, b]
  42. NIntegrate::ilim
  43. {{b -> 24831.554676346186` E^(-1.` NIntegrate[[Kappa][xp], {xp}])}}
  44.  
  45. Plot[Abs[yWKB[x]], {x, -L, L}]
Add Comment
Please, Sign In to add comment