Message Boards Message Boards


Solve a system of 2nd order PDE with Mathematica

Posted 4 months ago
1 Reply
0 Total Likes

Hi community,

I'm working on a drift-diffusion model which gives equations of the following type:

(a[0, 1] ' ' )[z] + 5.78 a[0, 1] [z] - 0.036 (a[0, 1] ' )[z] + (UnitStep[-1.63265 + z] - UnitStep[1.63265 + z]) == 0

wher a[ i , j ][z] are the expressions of interest (here as an example for i=0, j=1).

The most simple system of equations for a[0,1], a[1,1], a[1,2], a[2,2] can be found in the notebook.

The boundary condition should be n vanishes in both directions (it is a density).

I was trying NSolve, but unfortunately I only get oscillating results.

Is it possible that Mathematica has a problem with this type of DGL? I don't get any meaningful results for uncoupled equations either. What can I do ?

Best regards,


Wolfram Notebook

Welcome to Wolfram Community! Please make sure you know the rules:
Your post is too vague. Please EDIT your post and describe your subject extensively providing your code. You can format your code like that

int = Integrate[1/(x^3 - 1), x];
Map[Framed, int, Infinity]

You can also embed notebook or attach notebook.

enter image description here

Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
or Discard

Group Abstract Group Abstract