Evaluate numerical solution of PDE using output of pdepe

Syntax

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

Arguments

m

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

x

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.

ui

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

xout

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

Description

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

Notepdeval evaluates the partial derivative
∂u_{i}/∂x rather
than the flux f. Although the flux is continuous,
the partial derivative may have a jump at a material interface.