Expression des conditions aux limites d'un système EDP

bonjour à tous,
En voulant résoudre numériquement un système EDP parabolique suivant : http://hpics.li/dc99115
en utilisant Matlab avec la fonction Pdepe, je suis tombé sur le problème suivant: comment remplacer les termes intégrale des conditions aux limites des fonctions ni par une différentielle dni/dx pour que je puisse l'implémenter dans la fonction pdepe de Matlab dont les conditions aux limites sont en fonction de la différentielle et non de l'intégrale?
Merci d'avance de votre aide