Solve two first order partial differential equations numerically

Hello everyone,

I have to numerically solve a system of coupled first order partial differential equations. I am not posting the actual question here as it has large number of equations. But, a problem of similar nature is posted here. The system is a two coupled first order PDEs.

Where, s and t are the independant variables. The initial and boundary conditions are:

How can I solve the system of equations using Finite Difference Method or Method of Lines (MOL)? Which among these methods is more preferred for this problem?

Thank you.

[J. Alex Lee]

If you discretize s with finite differeneces, and advance the solution in t using an integrator, I believe that whole package is referred to as a method of lines (it's been a while I could be wrong) so it seems your question is not really an either-or.

I think matlab's pdepe will solve such problems.

[Sreejath S]

Thank you very much for the suggestion. I figured out that my equations are elliptic in nature. So pdepe should work. Am I correct?

[J. Alex Lee]

sounds likely, have you given it a shot?

[Sreejath S]

Yes. It worked thanks for the help!