First time here? Check out the FAQ!
answered 2016-07-28 14:09:48 +0100
Have you tried scipy's nquad function? I think that in your problem special care has to be taken in the line $v_1=v_2$, but it can be handled with the opts parameter.