quad3d.c 内のUser Definition 内のfunc関数を変更するだけで
任意の関数に対応する．
積分範囲を変えたければz1,z2(最内),yy1,yy2(中間)などを変える．
最外の積分範囲はmain.c 内の quad3d(func,0,2) の0と2を変更する．

//必要ならnr.h とnrutil.h とqgauss.c は Numerical Recipes in C を購入して手に入れてください！