How to solve a system of multiple integral equations in Python
I need help implementing the Mathcad code shown in the picture into python. This mathcad code is for Conjuctive Partial Absorbtion operators, and the equations that I am trying to find solutions for have multiple integrals contained within. I am at a loss as to how I could find values that satisfy the given requirements for this. The input parameters I am solving for are as follows:
r1 = 1 r2 = -3.5100640888436105 penalty = -25 reward = 15
Mathcad Code