Skip to content

Instantly share code, notes, and snippets.

@m0nhawk
Last active October 25, 2015 08:05
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save m0nhawk/4dc01861f609f71201d3 to your computer and use it in GitHub Desktop.
Save m0nhawk/4dc01861f609f71201d3 to your computer and use it in GitHub Desktop.
xmax = 3.;
zmin = 100.;
zmax = 121.;
tmax = 1.;
eqlist = {-170.2 (0.005 D[Ne[t, x, z], z] - D[Ne[t, x, z], {z, 2}]) -
3436. (0.005 D[Ne[t, x, z], x] - D[Ne[t, x, z], {x, 2}]) +
D[Ne[t, x, z], t] == 0.,
DirichletCondition[Ne[t, x, z] == 1, z == zmin || z == zmax],
DirichletCondition[Ne[t, x, z] == 1, Abs[x] == 3],
Ne[0, x, z] == 1. + 10. E^(-99.9 x^2 - 1. (-110. + z)^2)};
sol2 = NDSolve[
eqlist, {Ne}, {t, 0, tmax}, {x, -xmax, xmax}, {z, zmin, zmax}];
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment