Solving PDE's with Orthogonal Collocation and Matlab

viNcee

Hello people,

I am trying to solve a set of PDE's using Matlab by first applying the orthogonal collocation technique to reduce the PDE's to ODE's and subsequently solve the equations using ODE15s.

The code that I am using is based on what I've found on another website: http://faculty.washington.edu/finlay.../OC/OCcode.htm

I am having problems coding in my boundary conditions for the equations as it is given in the form of derivatives.
At x = 0, dT/dx = dc/dx = 0

Any help on the matter would be very much appreciated. Thank you.