answered 8 years ago
Have you tried scipy's nquad function? I think that in your problem special care has to be taken in the line v1=v2, but it can be handled with the opts parameter.