:ɭʽֵ()

ʽ:
[r,er,m,n]=IntegralxySimpson(f,xa,xb,ya,yb,er0,m0,n0,mN,nN)
[r,er,m,n]=IntegralxySimpson(f,xa,xb,ya,yb,er0,m0,n0,mN)
[r,er,m,n]=IntegralxySimpson(f,xa,xb,ya,yb,er0,m0,n0)
[r,er,m,n]=IntegralxySimpson(f,xa,xb,ya,yb,er0,m0)
[r,er,m,n]=IntegralxySimpson(f,xa,xb,ya,yb,er0)
[r,er,m,n]=IntegralxySimpson(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,59ҳʽ(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]=IntegralxySimpson(f,1.4,2,ya,1.5)//ִԼسõ½ע,yaÿһԪضӦ½Ӧ±ֵ

fy =
[ 0.64936403489515    0.57967296483187    0.50633319683294    0.42955452754621 ]
er =
[ 5.9074967140E-13    4.4819703504E-13    3.3661962106E-13    2.4824586830E-13 ]
m =
[ 64.0000000000000    64.0000000000000    64.0000000000000    64.0000000000000 ]
n =
[ 32.0000000000000    32.0000000000000    32.0000000000000    32.0000000000000 ]
