Hi !

In my last thread I think I couldn't make myself clear. So, when I replace the 1D 10-point gauss in the 3D quadrature with an adaptive algorithm the program is running, but it provides false results. I tried it with trapeziod, Simpson and Romberg algorithm but they all calculate the same false result. What can be the problem? I think it should work.
Can anybody help?