Advanced Search

Question: pdsolve bug or something else?

Question:pdsolve bug or something else?

J4James 355 Maple
numeric pdsolve
August 07 2013
0

Hello everyone,

restart:with(PDEtools):with(plots):

sys := {diff(T(x, t), t) = diff(T(x, t), x, x)+(diff(u(x, t), x))^2, diff(u(x, t), t) = diff(u(x, t), x, x)};

PDE:={diff(u(x,t),t)=diff(u(x,t),x$2),diff(T(x, t), t) = diff(T2(x,t), x)+(diff(u(x, t), x))^2,diff(T(x, t), x)=T2(x,t)};

BCs := {u(0,t)=sin(t), u(10,t)=0,T(0,t)=1, T(10,t)=0,u(x,0)=0,T(x,0)=0};

pds2 := pdsolve(PDE, BCs, numeric, spacestep=1/50);

I just want to plot T2(x,t) vs t for x=0.

pds2:-plot(T2(x,t), x=0, t= 0..5);

To me and @Carl Love, this seems to be some sort of a bug. 

What you think mates?

Thanks
Please Wait...