Student Support Forum: 'solving poisson equation with mathematica' topicStudent Support Forum > General > "solving poisson equation with mathematica"

 Next Comment > Help | Reply To Topic
 Author Comment/Response Benjamin 03/21/08 06:30am Hi, I'm currently trying to solve the following pde with rectangle boundary (I better already use mathematica code here): The equation (Poisson equation): D[u[y, z], y, y] + D[u[y, z], z, z] + 1 == 0 The boundary conditions on the rectangle(y in [0,0.1] and z in [-0.4,0.4]): u[y, 0.4] == 0, u[y, -0.4] == 0, u[0, z] == 0, Derivative[1, 0][u][0.1, z] == 0 I tried the following in mathematica: eqn = D[u[y, z], y, y] + D[u[y, z], z, z] + 1 == 0; //defining equation NDSolve[{eqn, u[y, 0.4] == 0, u[y, -0.4] == 0, u[0, z] == 0, Derivative[1, 0][u][0.1, z] == 0 }, u, {y, 0, 0.1}, {z, -0.4, 0.4}] //using NDSolve to solve the boundary problem But, as you might guess, I get an error using NDSolve the way above: "NDSolve::ivone: Boundary values may only be specified for one independent variable. Initial values may only be specified at one value of the other independent variable." When I lookup the description of that error I realize that in my example the following error condition matches: "This input generates a message because the equations specify values for the solution on all sides of the solution region." Does this mean, that mathematica can't approximate the solution to my problem with NDSolve or am I doing something wrong? So overall question is: How can I solve my problem using mathematica without writing my own numerical pde solver for my problem? Thanks in advance! URL: ,

 Subject (listing for 'solving poisson equation with mathematica') Author Date Posted solving poisson equation with mathematica Benjamin 03/21/08 06:30am Re: solving poisson equation with mathematica Henry Lamb 03/23/08 12:04pm
 Next Comment > Help | Reply To Topic