How to solve coupled systems of PDEs and ODEs using method of lines

I want to develop an 1D model for PEM fuel cell which includes several PDEs and ODEs. I have read that, it is possible to solve it using method of lines. I would like to get an example code so that i understand the implementation in Matlab

Answer 1

[Torsten]

http://www.scholarpedia.org/article/Method_of_lines/example_implementation

https://www.sciencedirect.com/science/article/abs/pii/S1364815203002111

[Sachin Hegde]

Thank you very much Torsten. I have one other question. Is it possible to solve multipoint boundary problems using MOL (i have regions with different propoerties)? I have seen it for bvp4c (Solve BVP with Multiple Boundary Conditions - MATLAB & Simulink - MathWorks Deutschland) but i would like to know if i can do similarly also using MOL as i have a PDE ODE system.

[Torsten]

Is it possible to solve multipoint boundary problems using MOL (i have regions with different propoerties)?

Including your problem in its mathematical formulation would make it possible for me to answer your question.

[Sachin Hegde]

Hi Torsten. I have attached a pdf which includes the system i am trying to solve. To summerise, it is a system of PDEs and ODEs. The system consists of several regions that are joined together (see the figure in the pdf) and each region has different physical properties however the governing equation remains the same. The boundary conditions are provided at the intersection of each region. I am struggling to implement this in matlab. I kindly ask your support and i thank you in advance.