:˹ֵ()

ʽ:
[r,er,m,n]=IntegralxyGauss(f,xa,xb,ya,yb,er0,m0,n0,mN,nN)
[r,er,m,n]=IntegralxyGauss(f,xa,xb,ya,yb,er0,m0,n0,mN)
[r,er,m,n]=IntegralxyGauss(f,xa,xb,ya,yb,er0,m0,n0)
[r,er,m,n]=IntegralxyGauss(f,xa,xb,ya,yb,er0,m0)
[r,er,m,n]=IntegralxyGauss(f,xa,xb,ya,yb,er0)
[r,er,m,n]=IntegralxyGauss(f,xa,xb,ya,yb)

f:ű,洢ֻxyΪƵĻֱʽ
xa:ֵ߾,ʾx
xb:ֵ߾,ʾx
ya:ֵ߾߰xķűʽ,ʾy
yb:ֵ߾߰xķűʽ,ʾy
er0:ʾƵ,ĬΪ1e-12
m0:ʾxлʱʼȷ,ĬΪ8
n0:ʾyлʱʼȷ,Ĭm0һ
mN:ʾxлʱȷ,ĬΪ1024
nN:ʾyлʱȷ,ĬmNһ

r:صĻֵֽ
er:ص
m:ص⵽ֵʱxֱĸ
n:ص⵽ֵʱyֱĸ

˵:ڱxa,xb,ya,ybдʱ,Щάһ,ҷصĽҲǶӦľ

ԭ:

1˫ػ,ɷֱÿһõĻ㷨ֵ
2㷨xyø˹ʽֵ֡

ο://.ֵ̳[M].ߵȽ,,2004,87ҳ3 

:

//֪ f(x,y)=ln(x+2*y)f(x,y)x1.42;yֱ0.7,0.8,0.9,11.5ֵ
//ʾ,f(x,y)xyһλֺɵõ ((x + 2*y)^2*(2*log(x + 2*y) - 3))/8

f="ln(x+2*y)";
ya=[0.7,0.8,0.9,1];
[fy,er,m,n]=IntegralxyGauss(f,1.4,2,ya,1.5)//ִԼسõ½ע,yaÿһԪضӦ½Ӧ±ֵ

fy =
[ 0.64936403492781    0.57967296484631    0.50633319683878    0.42955452754829 ]
er =
[ 1.7763568394E-15    1.4432899320E-15    7.7715611723E-16    7.2164496600E-16 ]
m =
[ 16.0000000000000    16.0000000000000    16.0000000000000    16.0000000000000 ]
n =
[ 8.00000000000000    8.00000000000000    8.00000000000000    8.00000000000000 ]
