Error solving non-linear BVP problem using NDSolve?

Hello everyone,

I am trying to solve the equations for a particulate filter by considering variable density with the pressure. I have attached the script.

Can someone help me?

The authors who solved this system used Matlab’s bvp5c. How can I convert this to Mathematica?

prova_densità_variabile.nb (25.1 KB)