Advertisement
Guest User

Untitled

a guest
Aug 17th, 2017
45
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.45 KB | None | 0 0
  1. ClearAll["Global`*"]
  2. Ri = 10.5/10^3;
  3. Ro = 11.5/10^3;
  4. h = 1/10^3;
  5. c = 30/10^3;
  6. l = 3/10^3;
  7. [Mu] = (4*Pi)/10^7;
  8. i0 = 1;
  9. f = 0.159*10^9;
  10. [Omega] = 2*Pi*f;
  11. [Sigma] = 1;
  12. Subscript[[CapitalLambda], 2] = Sqrt[[Gamma]^2 + I*[Omega]*[Mu]*[Sigma]];
  13. [ScriptCapitalL] = (([Mu]*Pi)/(h*(Ri - Ro))^2)*NIntegrate[(r*BesselJ[1, [Gamma]*r])*
  14. (R0*BesselJ[1, [Gamma]*R0])*((1 - E^((-h)*[Gamma]))^2/(E^(2*l*[Gamma])*[Gamma]^2))*
  15. ((E^(2*c*[Gamma] + 2*(c + 2*h + l)*Subscript[[CapitalLambda], 2])*
  16. (-1 + E^(2*c*Subscript[[CapitalLambda], 2]))*([Gamma] - Subscript[[CapitalLambda], 2])*
  17. ([Gamma] + Subscript[[CapitalLambda], 2])*
  18. ((-(E^(2*(2*h + l)*[Gamma]) + E^(2*c*Subscript[[CapitalLambda], 2])))*
  19. ([Gamma] - Subscript[[CapitalLambda], 2])^2 +
  20. (1 + E^(4*h*[Gamma] + 2*l*[Gamma] + 2*c*Subscript[[CapitalLambda], 2]))*([Gamma] + Subscript[[CapitalLambda], 2])^
  21. 2))/(E^(2*(c + 2*h + l)*([Gamma] + Subscript[[CapitalLambda], 2]))*
  22. ([Gamma] - Subscript[[CapitalLambda], 2])^4 +
  23. E^(2*(c + 2*h + l)*[Gamma] + 2*(3*c + 2*h + l)*Subscript[[CapitalLambda], 2])*
  24. ([Gamma] + Subscript[[CapitalLambda], 2])^4 - E^(2*c*[Gamma] + 2*(c + 2*h + l)*Subscript[[CapitalLambda], 2])*
  25. (2*E^(4*h*[Gamma] + 2*l*[Gamma] + 2*c*Subscript[[CapitalLambda], 2]) +
  26. (-1 + E^(2*c*Subscript[[CapitalLambda], 2]))^2)*([Gamma]^2 - Subscript[[CapitalLambda], 2]^2)^2)),
  27. {[Gamma], 0, Infinity}, {r, Ri, Ro}, {R0, Ri, Ro}]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement