PDE Solver for solving different equations in different subdomains

Hi MATLAB Community,

I am looking for a solver that solves different PDEs in 2 different domains (1-D) connected by a BC. I have a diffusion equation (species 1) in domain 1 and another diffusion equation (species 2) in domain 2. Concentration of species 1 and concentration of species 2 is connected to each other through a BC equation at their interface. I was looking into pde solvers in MATLAB that can deal with these scenarios but I couldn't find any solver that can accomodate an internal BC in the domain. Any suggestions/advice to solve these system of equations is very much appreciated.