Evaluate numerical solution of PDE using output of pdepe


[uout,duoutdx] = pdeval(m,x,ui,xout)



Symmetry of the problem: slab = 0, cylindrical = 1, spherical = 2. This is the first input argument used in the call to pdepe.


A vector [x0, x1, ..., xn] specifying the points at which the elements of ui were computed. This is the same vector with which pdepe was called.


A vector sol(j,:,i) that approximates component i of the solution at time tf and mesh points xmesh, where sol is the solution returned by pdepe.


A vector of points from the interval [x0,xn] at which the interpolated solution is requested.


[uout,duoutdx] = pdeval(m,x,ui,xout) approximates the solution ui and its partial derivative ∂ui/∂x at points from the interval [x0,xn]. The pdeval function returns the computed values in uout and duoutdx, respectively.

    Note   pdeval evaluates the partial derivative ∂ui/∂x rather than the flux f. Although the flux is continuous, the partial derivative may have a jump at a material interface.

See Also

Was this topic helpful?