On 3/23/11, michael.kogan at gmx.net <michael.kogan at gmx.net> wrote: > foo (x, y) := quag_qags (1/sqrt(1+(1/z-1)*x+(z^2-1)*y),z,0,1) [1]; I think you want quad_qags here, not quag_qags. best Robert Dodier