| 1 | initial version |
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.
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.