Advanced Search

Question: Derivative from numeric PDE soln

Question:Derivative from numeric PDE soln

raggapuffin 89 Maple
April 09 2010
1
> PDE := diff(u(x, t), t) = k*(diff(u(x, t), x))+diff(u(x, t), x, x); > a := 1/70; > IBC := {u(0, t) = 0, u(5, t) = 0, u(x, 0) = exp(-(x-2.5)^2/(2*a^2))/((2*Pi)^(1/2)*a)}; > pds := pdsolve(PDE, IBC, numeric, time = t, range = 0 .. 5, timestep = 1/2000, spacestep = 1/2000); >R := pds:-value(output = listprocedure); >uu := subs(R, u(x, t)) I now have a numeric solution for the PDE, uu(x,t) which I can integrate but does anyone know how to obtain diff(uu,t)? It seems to do it automatically for an ODE soln. The goal of this being int(diff(uu,t),x=0..5), if anyone has any ideas I'd like to hear them, Ryan
Please Wait...