I am a little confused by the "pdepe" help file, for two reasons:
1. How can I assign continuity conditions at the boundaries where the discontinuities occur? The "bvp4c" help file explains how to do this quite clearly for ODEs (see the "Multipoint Boundary Value Problems" section in http://www.mathworks.com/help/techdoc/math/f1-663004.html), but the "pdepe" help file does not mention how to assign continuity conditions. Also, the example file for discontinuous PDEs ("pdex2", which can be accessed by the command "edit pdex2"), does not assign any continuity conditions at the discontinuous boundaries. Why is this?
2. The "pdepe" help file states that discontinuities "are permitted provided that a mesh point is placed at each interface". I'm assuming this means that the mesh point at the interface does not have to be repeated, as in "bvp4c" (see the "Multipoint Boundary Value Problems" section in http://www.mathworks.com/help/techdoc/math/f1-663004.html).