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

I tried to input these values and expected values for w1 and w2 that satisfied the equations. However, I was unable to get them and need a way to replicate the pictured mathcad code efficiently in python.