dblquad or similar in Octave?

[Jens Benecke]

Hi everybody,

I need to perform a double quadrature in Octave. In Matlab, this is the
code:

function retval=induk(z1,z2,r1,th)
    tol  = [];
    retval = r1* dblquad(@(p1,p2) ...
        (sin(p1).*sin(p2).*cos(th)+cos(p1).*cos(p2))./ ...
        sqrt( (z2*sin(th)+cos(p2)*cos(th)-r1*cos(p1)).^2.+...
              (sin(p2)-r1*sin(p1)).^2. +...
              (z2*cos(th)-cos(p2)*sin(th)-z1).^2. ), ...
   0, pi, 0, 2*pi, tol, @quadl)/(2.*pi)
end


I have searched the net but all I could find was several references
to "missing functions" in Octave (compared to Matlab).

I use Octave 2.1.71 under Linux (SuSE and Ubuntu). How would I perform this
kind of quadrature with Octave?

Thank you!


-- 
Jens Benecke