Simulation of coupled partial differential equation using MATLAB

I am trying to simulate a model for ammonia stripping using a mass transfer system. The equations are as given in the attachment.

I need to find out profile for ammonia concentration in liquid and gas phase. Please suggest a method which I can use. If you have any code which can be used, it would be of great help if you can give me that too.