I am trying to solve the coupled differential equation, but i am having trouble finding the constants. I have successfully found the expressions for u[x], p[x] and [Mu][x]. and use the Solve function and the boundary conditions given by eqn 3-7 to find the values of the constants. the result is empty, does that mean the boundary conditions are insufficient, or is there some other reason? Also, could anyone suggest how to simplify the code that I wrote, it looks a bit messy :(