Items tagged with pdsolve

hello,I want to solve a quesstion about heat equation,that the quesstion like this:

I use the code like this

but the results is wrong obviously and what's wrong with this code?

anxious for your help,thanks.

 

 

I want to solve the following non linear PDE

SS := [diff(u(x, y, t), t)-0.625e-1*(diff(u(x, y, t), x, x)+diff(u(x, y, t), y, y))-6*(diff(u(x, y, t)*(diff(v(x, y, t), x)), x)+diff(u(x, y, t)*(diff(v(x, y, t), y)), y)) - 2*(u(x, y, t))(1-u(x, y, t))=0, diff(v(x, y, t), t)-(diff(v(x, y, t), x, x))-(diff(v(x, y, t), y, y))+16*v(x, y, t) -u(x, y, t)=0]

when i use the command

sol := pdsolve(SS, [u, v], singsol = false)

maple give the error message

Error, (in pdsolve) found the independent variables {t, x, y} also present in the names of the functions of the system []

 

how i can save results on text file from pdsolve?

 

thanks

maple11.mw
 

PDE := diff(u(y, t), t) = diff(u(y, t), y, y)-0.9e-1*(diff(u(y, t), y, y, y, y))

diff(u(y, t), t) = diff(diff(u(y, t), y), y)-0.9e-1*(diff(diff(diff(diff(u(y, t), y), y), y), y))

(1)

IBC := {(D[1, 1](u))(0, t)-0.9e-1*(D[1, 1, 1, 1](u))(0, t) = 0, u(0, t) = 1, u(100, t) = 0, u(y, 0) = 0, (D[1, 1](u))(100, t) = 0}

{(D[1, 1](u))(0, t)-0.9e-1*(D[1, 1, 1, 1](u))(0, t) = 0, u(0, t) = 1, u(100, t) = 0, u(y, 0) = 0, (D[1, 1](u))(100, t) = 0}

(2)

 

 

 

 

``

 

pds := pdsolve(PDE, IBC, numeric, u(y, t), spacestep = 1/5, timestep = 1/5)

_m2865547563008

(3)

 

     

   

NULL

_m2865547563008

(4)

 

 

fname1 := "C:/Users/test/Desktop/txtop1.txt"; seq(fprintf(fname1, "%12.8f, %12.8f\n", y, u(y, .1), y = 0 .. 10, 1)); fclose(fname1)

Error, (in fprintf) file or directory does not exist

 

p1 := pds:-plot(t = 0):p2 := pds:-plot(t = 1/10):

plots[display]({p2}, title = `\` profile at t=0,0.1`)

 

NULL


 

Download maple11.mw

 

Hi all,

 

I have a partial differential equation similar to the following:

Equation: f_x(x,y) + f_y(x,y) = f(x,y) + f(x,0),
Boundary value conditions: f(x,10) = f(10,y) = 0.

The solution is that f is identically equal to 0.

 

However, I am having trouble solving this equation in Maple. I type the following:

pde := diff(f(x, y), x)+diff(f(x, y), y) = f(x, y)+f(x, 0);

bv1 := f(x, 10) = 0;

bv2 := f(10, y) = 0;

solution := pdsolve(pde, {bv1, bv2}, numeric, time = x, range = 0 .. 10);

 

When Maple tries to evaluate the last expression, I get the error

Error, (in pdsolve/numeric/process_PDEs) PDEs can only contain dependent variables with direct dependence on the independent variables of the problem, got {f(x, 0)}

 

It seems to have difficulties with the expression "f(x,0)". Is there some trick to typing this in a way that makes Maple interpret it correctly?

 

Edit: I encounter the same problem, when I try to solve the ODE f'(x) = f(x) + f(0), where f(10) = 0.

 

Best regards.

Hello all

 

I am new to Maple, and I am solving a system of two coupled partial differential equations using pdsolve, but I am having a hard time retrieving the solution evaluated at some point from the output. The output of pdsolve is a module, which appears to have different "methods" on it, including "plot3d" and "value". I can easily get a plot of my solution by using plot3d, but I don't know how to get a meaningful value out. For instance, if my solution is (f(x,y), g(x,y)), I would like to define H(x,y) = (f(x,y), g(x,y)), and be able to type H(10,10) into Maple to have my solution evaluated at that point. The result should be (1,1).

Here is a toy example:

firstEq := diff(f(x, y), x)+diff(f(x, y), y) = f(x, y)+g(x, y);
secondEq := diff(g(x, y), x)+diff(g(x, y), y) = 2*f(x, y)+g(x, y);
pdsystem := {firstEq, secondEq};

bv11 := f(10, y) = 1;
bv12 := f(x, 10) = 1;
bv21 := g(10, y) = 1;
bv22 := g(x, 10) = 1;
bvs := {bv11, bv12, bv21, bv22};
 
pdsolution := pdsolve(pdsystem, bvs, numeric, time = x, range = 0 .. 10);
 
pdsolution:-plot3d(x = 1 .. 10, y = 0 .. 10);
pdsolution:-value(10, 10);
Error, (in pdsolve/numeric/value) got additional unknown arguments {2}
 

Best regards.

 

 

Hi all. I need to solve a large number of PDEs (partial differential equations) symbolically and simultaneously, to find the linearly independant answers for all of them, I use

ans := pdsolve({seq(PDE[i]=0,i=1..d)});

The PDEs are all linear first order and it is very easy to solve them one by one by hand, but in some cases I have 100 PDEs or more, so Maple is either very slow or doesn't work. For d=120, it was evaluating for hours without a result.

 

For example I have d=120 PDEs, and 200 variables. It works for d=30 of them (takes 13 minutes on my 16GB RAM  windows 7 computer). So if I do this:

ans1 := pdsolve({seq(PDE[i]=0,i=1..30)});

ans2 := pdsolve({seq(PDE[i]=0,i=31..60)});

ans3 := pdsolve({seq(PDE[i]=0,i=61..90)});

ans4 := pdsolve({seq(PDE[i]=0,i=91..120)});

Then how can I have only one vector of the linearly independent answers of all of them?

 

And in general, is Maple supposed to do this kind of calculations at all?

If yes, do you have any ideas on how to improve this procedure? 

If not, do you know in which software or programming package I can solve a large number of PDEs symbolically?

Your help is much appreciated. 

How can I find the numeric value of this PDE:

Where

I want the numeric values i.e theta(2,2) , theta (1,3) and the likes.

Hi! 

Tried to solve the PDE below (q and p are time-dependent variabels, q(t),p(t)):

pde := diff(rho(t, q, p), t) = -(diff(rho(t, q, p), q))*p+(diff(rho(t, q, p), p))*(2*q+2);
          

pdsolve(pde, rho(t, q, p));
           

And got the answer: 

rho(t, q, p) = _F1(p^2+2*q^2+4*q, -(1/2)*sqrt(2)*arctan((q+1)*sqrt(2)*(1/sqrt(p^2)))+t)

But I'm not sure how to interpret the result. I understand that  _F1 is an arbitrary function, but then I get confused with the comma? I thought that I'd get a function of q and p, where they depend on t. 

Best regards
Sannis

 

I have the following code in Maple 13:

SYS := {diff(T(x, t), t) = diff(T(x, t), x, x)};
IBC := {T(1, t) = 0, T(x, 0) = 1, (D[1](T))(0, t) = -exp(t)};
SOL := pdsolve(SYS, IBC, numeric, time = t, timestep = 1/10);
R := SOL:-value(output = listprocedure); temperature := subs(R, T(x, t));

It integrates the heat equation in the interval x=0..1 and it seems to work ok.
However, I have problems trying to obtain the temperature derivative at the boundaries (this is at x=0 and x=1). I'm using different commands and it seems the derivative is evaluated ok inside the domain but not at the boundaries.
If I try
fdiff(temperature(x, t), [x], {x = 1,t=0.5});

or

evalf((D[1](temperature))(1, .5));

I don't get any numerical answer. Any idea how I could obtain the value of the derivative at the boundaries?

Thanks in advance,

Javier

 




 

 

Hello everybody,

 

I want to solve this pde.the desire solution is V(r,z). three boundar conditions are written that two of them are related to rhe radial and one is related to the longitudinal coordinate.

I attached the solution for you. but this solution is derived by Matlab. Now, I just want to resolve it by Maple, but I couldn't reach it. Please let me know the correct way asap.

 

Thanks a lot.PDE.mw
 

NULL

NULL

restart 

Q := diff(V(r, z), r, r)+(diff(V(r, z), r))/r-V(r, z)/r^2+diff(V(r, z), z, z)+C/r = 0

diff(diff(V(r, z), r), r)+(diff(V(r, z), r))/r-V(r, z)/r^2+diff(diff(V(r, z), z), z)+C/r = 0

(1)

NULL

NULL

pdsolve(Q)

PDESolStruc(V(r, z) = _F1(r)*_F2(z)-(1/2)*(_C1/r+_C2*r+_C3*r*ln(r))*C/_C3, [{diff(diff(_F1(r), r), r) = _F1(r)*_c[1]+(-(diff(_F1(r), r))*r+_F1(r))/r^2, diff(diff(_F2(z), z), z) = -_F2(z)*_c[1]}])

(2)

NULL

 

NULL

NULL

``

``


 

Download PDE.mwPDE.mw

 

 

For my fdiff graph, it seems that the cirtical points appear to be jaggered or not smooth. Anyone nows what seem to be the problem? i tried increasing the numpoints but it did not work:( I am open to all opinions. Thanks:)

 

fyp2.mw

Hi all. I'm using pdsolve to find linearly independent solutions for two partial differentioal equations. I can solve them seperatly by "pdsolve", but when I try to solve them at the same time to get answers that fit in both of the eauations, I get an error message. I need to be able to solve them as a set of equations, as I have bigger examples. Here are the equations and outputs: 

PDE[1] := diff(f(theta[I], theta[A], theta[B], theta[AB]), theta[AB]);

PDE[2] := -(diff(f(theta[I], theta[A], theta[B], theta[AB]), theta[I]))+diff(f(theta[I], theta[A], theta[B], theta[AB]), theta[A])+diff(f(theta[I], theta[A], theta[B], theta[AB]), theta[B]);

pdsolve(PDE[1] = 0)  ----->  f(theta[I], theta[A], theta[B], theta[AB]) = _F1(theta[I], theta[A], theta[B])

pdsolve(PDE[2] = 0) ----->  f(theta[I], theta[A], theta[B], theta[AB]) = _F1(theta[I]+theta[A], theta[I]+theta[B], theta[AB])

pdsolve({seq(PDE[i] = 0, i = 1 .. 2)}) ---->  Error, (in PDEtools:-casesplit) too many levels of recursion

Why doesn't it produce the answer, although I know that it must be:

f(theta[I], theta[A], theta[B], theta[AB])=_F1(theta[I]+theta[A], theta[I]+theta[B])

Can somebody of Maple users execute the following command

restart; pdsolve({diff(w(x, y, z), x)+diff(w(x, y, z), y, y)+2*(diff(v(x, y, z), x)-(diff(u(x, y, z), y))-2*w(x, y, z)) = diff(w(x, y, z), z, z), 3*(diff(u(x, y, z), x, x))+2*(diff(u(x, y, z), y, y))+2*(diff(v(x, y, z), x, y))+2*(diff(w(x, y, z), y)) = diff(u(x, y, z), z, z), 3*(diff(v(x, y, z), y, y))+2*(diff(v(x, y, z), x, x))+2*(diff(u(x, y, z), x, y))-2*(diff(w(x, y, z), x)) = diff(v(x, y, z), z, z)}, {u(x, y, z), v(x, y, z), w(x, y, z)})

restart; pdsolve({diff(w(x, y, z), x)+diff(w(x, y, z), y, y)+2*(diff(v(x, y, z), x)-(diff(u(x, y, z), y))-2*w(x, y, z)) = diff(w(x, y, z), z, z), 3*(diff(u(x, y, z), x, x))+2*(diff(u(x, y, z), y, y))+2*(diff(v(x, y, z), x, y))+2*(diff(w(x, y, z), y)) = diff(u(x, y, z), z, z), 3*(diff(v(x, y, z), y, y))+2*(diff(v(x, y, z), x, x))+2*(diff(u(x, y, z), x, y))-2*(diff(w(x, y, z), x)) = diff(v(x, y, z), z, z)}, {u(x, y, z), v(x, y, z), w(x, y, z)})

Error, (in simplify/normal) Maple was unable to allocate enough memory to complete this computation.  Please see ?alloc

 

 

 

in Maple 2016.1.1 on a powerful comp and report the obtained result as an answer to the question?
 That would be very kind of her/him. Thanks in advance. 

Download pdsolve.mw

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

Hello, 

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! 

test.mw

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