Thank you very much for pointing this out to me! q was never meant to take values outside [0,2], so my Vq was ill-defined. I added min(1/E,2) in the upper limit in Vq, and it worked out just fine. Thanks again, Sturla

Thank you very much for pointing this out to me! q was never meant to take values outside [0,2], so my Vq was ill-defined. I added min(1/E,2) in the upper limit in Vq, and it worked out just fine. Thanks again, Sturla