I am looking at a system of nonlinear simultaneous equations. The two variables are u>0 and b>0. How can I solve this problem in Matlab, in Python, or in Fortran? Thanks.
You can easily eliminate one of those equations by solving #1 for b. Then use that to solve #2 for u.
You're going to have to use an iterative method to do it: guess a solution, calculate an estimate, compare to your guess, adjust and repeat until you converge.
I'd use numerical integration (5th order Runge-Kutta or something else) to calculate the integral.