一使用两次一重积分%%二重积分f=(x,ｙ)exp（ｓin(ｘ）)*ln（ｙ),y从5*ｘ积分到ｘ＾2,ｘ从1０积分到2０1（7、X后版本才有此函数quａd2d)y1＝ｑuaｄ2d(（x,y）exｐ(siｎ(x）)、*log(y)，1０,20,(x)5*ｘ,(x)x、＾２)２y2＝ｑuａdl((x)aｒrａｙfuｎ((ｘ)ｑuadl（(y)exｐ(sin（ｘ))、*ｌｏg（ｙ)，5*x,x、^２),ｘ),１0,2０)3y3＝dblquad（(x,y)ｅｘｐ（ｓｉn(x)）、*log(ｙ)、*（y>=5*x＆y<=x、^2)，10，20,５0,４00)详细请瞧吴鹏老师得文章HＹPERLINＫ""＼t"＿blank"二使用ｄblｑuad函数q=dblｑｕad(fuｎ,xmiｎ，ｘmax,ymｉｎ,ymaｘ，tｏｌ，mｅｔｈod)该函数求f(x,y)在[ａ，b]×［c,d]区域上得二重定积分。参数ｔol,tｒａｃe得用法与函数quaｄ完全相同。例8－5计算二重定积分(1)建立一个函数文件fxy、m:funcｔioｎf＝ｆxy(ｘ,y)globaｌkｉ;ｋｉ＝ｋi+1；%ki用于统计被积函数得调用次数ｆ＝ｅxp(－ｘ、^2/2)、*sｉn(x、＾2+y）;(２)调用ｄblquaｄ函数求解。ｇloｂaｌｋi;ki=０；Ｉ=dblｑuad('ｆxｙ',-2，2,-1,1)kｉＩ＝1、574４931８97４４94ｋi=１0３８来源精通MATLAB科学计算一书王正林,龚纯,何倩编写,电子工业出版社三复合辛普森公式(矩形积分区域)funｃtionq=DblSimpson(f，a，A，ｂ,B,m，ｎ)ｉｆ（m==1&&n==1）％辛普森公式q＝((B-b)*(Ａ－a）/9）＊（suｂs(sym(f),ｆinｄsｙm(sym(ｆ)),{a，b})+、、、ｓuｂs（sｙm（f)，finｄsym（sym(f）),{a,B})+、、、subs(sｙｍ（f）,ｆiｎｄsym(ｓym(f)），｛A,b})+、、、subs(sym(ｆ),fiｎｄsｙm(sｙｍ(ｆ)),｛Ａ,Ｂ})+、、、4＊ｓuｂｓ(syｍ(f),findｓｙm(ｓym(f)）,｛(A-a)／2,b})+、、、4*subs(sym(f）,finｄsｙm（ｓym(f）),{(A-ａ)／2,B})+、、、4＊sｕbｓ(sｙm（ｆ),finｄsym（sym(f）),{a,（B-b)/2})+、、、4*ｓｕbs(ｓym(ｆ),ｆinｄｓym(syｍ(ｆ)),｛A,(Ｂ-b)/２})＋、、、1６*sｕbs（sym(f),findｓym(sym（ｆ)),{(A-a)/2，(B－b)／２}));elsｅ%复合辛普森公式ｑ=0;fori=０:ｎ－1forj=0:ｍ-１x＝a+2*i＊(A-ａ)/２/n;y=b＋2*j*(Ｂ-b)/2／m;ｘ1=ａ+(２*i＋1)*(A-a)/2/n;y1=ｂ+(2*ｊ+1)＊(B-b)/2/m;x2＝a＋2*(i+1)*(Ａ-a）／2/ｎ;ｙ２=ｂ+2＊（j+１)*(B-b）/2/ｍ;q＝q+suｂs(sym(f）,fｉnｄsym（sym(f）),{x,y})+、、、subs（sym(f),fｉndsym(sym(f)),｛x,y2})＋、、、suｂｓ(sｙm(f）,findｓyｍ(syｍ(f)),{x2,y})+、、、subｓ(sｙm(f),findsym(ｓym(f)）,{x２，y2｝)+、、、4＊ｓｕbｓ(ｓｙm(f),fｉndｓym(sym(ｆ)),{x，y1})+、、、4*subｓ(sym(ｆ），ｆｉndsym(ｓym(ｆ）)，{ｘ2,y1})+、、、４*subｓ(sｙm(f),fｉｎｄsym(ｓｙm（f）)，｛x1,y})+、、、4＊subs(ｓym(f）,findsyｍ（sｙm(ｆ)），{x1,y2})+、、、16*sｕbs（syｍ（ｆ),fｉndｓyｍ(ｓｙm(f））,{x１,ｙ1｝）；eｎdendｅｎdq=((B-b）＊(A-a)/３6/m／n)*q;四复合梯形公式(矩形积分区域)fuｎctionq=DｂlTrａpｒｌ(f，ａ,A,ｂ,B,m，n)ｉf(ｍ=＝１&&n=＝1)％梯形公式q=((Ｂ-b）*(A-a)/4)*(ｓuｂs(sｙｍ(f)，findsyｍ（syｍ(f)),{a,b})+、、、subｓ(sｙｍ(f),findsｙm(syｍ(f）),{a,B})+、、、subs(ｓyｍ(f),fiｎdsym(sym(f）),{A，b｝)+、、、subs(sｙm(f),fｉndｓｙm(sym(ｆ)),｛Ａ,B}));els