Message Boards Message Boards

GROUPS:

Solve a system of 2nd order PDE with Mathematica

Posted 4 months ago
535 Views
|
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,

Michael&

Wolfram Notebook

Welcome to Wolfram Community! Please make sure you know the rules: https://wolfr.am/READ-1ST
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
Attachments
Remove
or Discard

Group Abstract Group Abstract