area and volume integrals with quadpack

 whoops, meant to include qvolume code:

qvolume(ex,uL,vL,wL) :=
 block([u,u1,u2,v,v1,v2,w,w1,w2,rL,uvwerr],
      [u,u1,u2] : uL,
      [v,v1,v2] : vL,
      [w,w1,w2] : wL,
   
      rL : qgsz(qgsz(qgsz(ex,w,w1,w2)[1],v,v1,v2)[1],u,u1,u2),
      if not listp(rL) then
               (print(" returns noun form "),return(false)),
      uvwerr : rL[4],
      if uvwerr = 0 then rL[1]
        else (print(" error = ",part(codelist,uvwerr,2)),
                                  false))$

------------------------
Ted Woollett