Guest User

Untitled

a guest
Feb 19th, 2018
82
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.87 KB | None | 0 0
  1. gg[x_] := 1 - x^2/2
  2. region = ImplicitRegion[x > -1 && x < 1 && y < gg[x] && y > -gg[x], {x, y}];
  3. NDEigensystem[{-Laplacian[{ux[x, y], uy[x, y]}, {x, y}],
  4. DirichletCondition[uy[x, y] == 0, x == -1],
  5. DirichletCondition[uy[x, y] == 0, x == 1],
  6. DirichletCondition[{ux[x, y], uy[x, y]}.{1, D[gg[x], x]} == 0,
  7. y > 0 && x > -1 && x < 1],
  8. DirichletCondition[{ux[x, y], uy[x, y]}.{1, -D[gg[x], x]} == 0,
  9. y > 0 && x > -1 && x < 1]}, {ux[x, y],
  10. uy[x, y]}, {x, y} [Element] region, 6]
  11.  
  12. region = ImplicitRegion[x^2 + y^2 + z^2 <= 1, {x, y, z}];
  13. g1 = (x* Ex[x, y, z] + y *Ey[x, y, z] + z* Ez[x, y, z])/(x^2 + y^2 + z^2);
  14. NDEigensystem[{-Laplacian[{Ex[x, y, z], Ey[x, y, z], Ez[x, y, z]}, {x,
  15. y, z}], DirichletCondition[{Ex[x, y, z] - x*g1 == 0,
  16. Ey[x, y, z] - y*g1 == 0, Ez[x, y, z] - z*g1 == 0}, True]}, {Ex[x, y, z], Ey[x, y, z], Ez[x, y, z]}, {x, y, z} [Element] region, 6]
Add Comment
Please, Sign In to add comment