solve coupled nonlinear equations (integral inside)