Solving Coupled PDE with zero in denominator

Question

I want to solve a system of coupled 1d pde. but my problem is that the initial values and the boundary conditions cause a zero in the denominator on RHS of the equation.