Ask Your Question

Revision history [back]

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.