To solve one dimensional model, we use a finite volume method applied to the PDE written in conservative form. Most of the one dimensional models can be rewritten in this form. We have implemented the method in Matlab taking advantage of the Symbolic toolbox to avoid deriving the jacobians by hand. The method is efficient and versatile. For this reason, we have used it in many applications. If you want to learn more about the method, you can find more details in our publication.