Items tagged with pde

Hi all, I need a heep to solve 6 PDEs from the first order with initial and boundry conditions.

Here is the file

Hello there. I have a simple differential equation (diff(y(x), x)+1)*exp(y(x)) = 4 , in which I want to got a partial solution (find a value of constant) for y(1)=0. How do I do that? I've checked help on pdsolve function and it's too complicated for me to understand something. Examples do not contain anything such simple as that. Thanks in advance.

Hi guys! 

I have a PDE system. The mayority of the equations are equal to zero, but two of them are:


where a, b, c, d are CONSTANT parameters. I know that if a=b=d=c=1 the system is inconsistent. But I also know that if a=-1, b=d=0 and c=1 the system is consistent and exist the solution. I wanna know if there's a way to ask maple to find another selections of my parameter that make my PDE consistent and what it's the solution for that selection of a,b,c,d. 

Here's my PDE system (sys2).

thank you so much for your time! 

I am trying to solve a system of equations of motion of gravitational field and in this way, I deal with a second order differential equation containing Dirac delta function as follows:

Eq9 := -(l[f]^2/l[g]^2+2)*(diff(G(r, t), r, r)-2*(diff(G(r, t), r))/r)+l[f]^2*(diff(F(r, t), r, r)-2*(diff(F(r, t), r))/r)/l[g]^2+2*l[f]^2*(m^2*l[f]^2-3)*(G(r, t)/r^2-F(r, t)/r^2)/l[g]^2 = l[g]*E*exp(Pi*t*l[f]^2*(1-2*kappa)/l[g]^2)*Dirac(t)*Dirac(r)*Dirac(z)


Eq10 := diff(G(r, t), r, r)-2*(diff(G(r, t), r))/r-2*kappa*(diff(F(r, t), r, r)-2*(diff(F(r, t), r))/r)-(2*(m^2*l[f]^2-3))*(G(r, t)/r^2-F(r, t)/r^2) = l[g]*E*exp(Pi*t*l[f]^2*(1-2*kappa)/l[g]^2)*Dirac(t)*Dirac(r)*Dirac(z)

 want to solve these equations with MAPLE software symbolically.
Can anyone guide me in this way, please?

I had problems plotting the solutions of my PDE with the derivative as the axis. Anyone knows what went wrong?

p1 := pds:-plot(diff(f(x,y),y), y = 0 .. 10, x = 0, numpoints = 100)

I seem to be getting a plot of y=0. Any help would be greatly appreciated:))

I know there is some code to compute the symmetry of pde, but i want to compute the symmetry of some differential expression, as f(x,t,ux,ut).

The difference is that when solve the infinitesimal, we don't have the identity f(x,t,ux,ut)=0, we just take each variable in the jet space as independent.

Is there any package to do that?

help and thanks!

I have a system of pde as follow,

PDE := [(x*y+1)*(diff(f(x, y), y, y, y))+(x+(3/4)*f(x, y))*(diff(f(x, y), y, y))-(1/2)*(diff(f(x, y), y))^2+T(x, y) = (1/4)*x*(diff(f(x, y), y))*(diff(diff(f(x, y), y), x))-(1/4)*x*(diff(f(x, y), x))*(diff(f(x, y), y, y)), (x*y+1)*(diff(T(x, y), y, y))/(.733)+(x/(.733)+(3/4)*f(x, y))*(diff(T(x, y), y)) = (1/4)*x*(diff(f(x, y), y))*(diff(T(x, y), x))-(1/4)*x*(diff(f(x, y), x))*(diff(T(x, y), y))]

sys_ode := diff(g(y), y, y, y)+(3/4)*g(y)*(diff(g(y), y, y))-(1/2)*(diff(g(y), y))^2+h(y) = 0, (diff(h(y), y, y))/(.733)+(3/4)*g(y)*(diff(h(y), y)) = 0

ics := g(0) = 0, h(0) = 1, (D(g))(10) = 0, g(10) = 0, h(10) = 0

sol2 := dsolve([sys_ode, ics], numeric)

BC := {T(0, y) = h(y), T(x, 0) = 1, T(x, 10) = 0, f(0, y) = g(y), f(x, 0) = 0, f(x, 10) = 0, (D[2](f))(x, 0) = 0}

pds := pdsolve(PDE, BC, numeric)

module() ... end module

pds:-plot(T, y = 0 .. 10, x = 0);

Error, (in pdsolve/numeric/plot) unable to compute solution for x<HFloat(0.0):
solution becomes undefined, problem may be ill posed or method may be ill suited to solution

When I try to use the solution of the Ode as the boundary condition for PDE, by subbing g(y) and h(y) into BC. The plot returns me the error. Anyone knows the reason behind this and how to solve? Any help would be really greatly appreciated. Thanks



I have a PDE which I would like to group certain terms in common. They are all general functions and do not all share variables in common.

It is easier to explain in the script attached.

Any guidance anyone can offer would be greatly appreciated. 

Thank you


I am trying to substitute an general function pde expression I have derived in maple into another general function pde expression but cannot seem to find a method that works. 

Attached is an example worksheet - any help offered would be much appreciated.



I have a PDE system. When I use pdsolve it gets me the messege " pdsolve->Warning: System is inconsistent". Is there a way I can see which equations breaks the system down? 
For this system, it's difficult to see from ayeball where the problem is. 
Thank you!

Dear all

I have a PDE and its analytical solution. I want to find the numerical solution by Finite Difference Method.

I duscratize the PDE and boundary condition and Could not able to solve them togethe.

Here is the file

Please, I need assistance with this problem.

Here is the problem I am trying to solve:


myPDE1 := D11*diff(w(x,y), x$4) + 2*(D12+2*D66)*diff(w(x,y), y$4) + D22*diff(w(x,y), x$2, y$2) - G*diff(w(x,y), x,y)= 0;


pdsolve(myPDE1, build);

"Boundary conditions";
"(Note:the domain for the problem is a rectangle)";
bc1 := w(0,y) = 0; # @ x=0 edge;
bc2 := w(a,y) = 0;  # @ x=a edge;
bc3 := w(x,0) = 0; # @ y=0 edge;
bc4 := w(x,b) = 0; # @ y=b edge;
bcx1 := -D11*D[2](w)(0,y) - D12*D[2](w)(0,y) = 0; # @ x=0 edge;
bcx2 := -D11*D[2](w)(a,y) - D12*D[2](w)(a,y) = 0; # @ x=a edge;
bcy1 := -D12*D[2](w)(x,0) - D22*D[2](w)(x,0) = 0; # @ y=0 edge;
bcy2 := -D12*D[2](w)(x,b) - D22*D[2](w)(x,b) = 0; # @ y=b edge;

sol := [myPDE1, bc1, bc2, bc3, bc4, bcx1, bcx2, bcy1, bcy2];


and D11, D12, D22, D66 and G are constant.
The intention is to find the critical value for G"

I need help with how I can handle the boundary conditions for the problem. Thanks a million.

Hello everybody!
I have a PDE with initial and boundary conditions. I want to plot its solution by taking "t" as x axis. I have seen the documentation. It only has the space variable on x axis. Please show me a way to achieve what I intend.

here is the file

The value of x can be chosen as 0.16 or 0.21

I am trying to solve a PDE using pdsolve-numeric. I am getting an error related to boundary conditions.
Please see the follwing worksheet and suggest me some solutions

Please help me.I don't know how to achieve the following iteration relation by maple code.iterative relationship



1 2 3 4 5 6 7 Last Page 1 of 13