:˹ʽֵ()

ʽ:

[r,er,n]=IntergalMultiple(f,x,x1_a,x1_b,x2_a,x2_b,,xn_a,xn_b,er0,n0,n_Max)
[r,er,n]=IntergalMultiple(f,x,x1_a,x1_b,x2_a,x2_b,,xn_a,xn_b,er0,n0)
[r,er,n]=IntergalMultiple(f,x,x1_a,x1_b,x2_a,x2_b,,xn_a,xn_b,er0)
[r,er,n]=IntergalMultiple(f,x,x1_a,x1_b,x2_a,x2_b,,xn_a,xn_b)

f:ű,洢ֱʽ
x:ű,洢ֵı,ÿԶŷָ
x1_a,x2_a,,xn_aֱӦxеĵ12n,ע˱ΪݿΪʽ,ǵһĻx1_aΪ
x1_b,x2_b,,xn_bֱӦxеĵ12n,ע˱ΪݿΪʽ,ǵһĻx1_bΪ
er0:ʾƵ,ĬΪ1e-12
n0:Ϊݻ߾,ʾxÿֵʱ,ĬΪ4
n_Max:Ϊݻ߾,ʾxÿֵʱĸ,ĬΪ4n0

r:صĻֵֽ
er:ص
n:ص⵽ֵʱxиĸ

˵:

1ֺֺ֮ͬһǿɶغֵ,˲һءءء
2ΪػʹñȽ,˱һءػֺֻ֧ΪҲ˵,ÿִһκֻһޡ
3Żά,ԼõĸӼʱԽ׷󾫶Ƚϸߵ,齫n0ƫ

ԭ:

1ػ,ɷֱÿһõĻ㷨ֵ
2㷨xеÿ˹㷨ֵ֡

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

:

//Ժf=x^2*y*z,xyzֱΪ[0,1],[0,2*x+3],[0,x+4*y+1]עⱾ׼ȷΪ183821/504364.724206349206

f="x^2*y*z";x="x,y,z";[a,er,n]=IntergalMultiple(f,x,0,1,0,"2*x+3",0,"x+4*y+1")//سõ½
a =
[ 364.724206289715 ]
er =
[ 0.00000000000000 ]
n =
[ 8.00000000000000    8.00000000000000    4.00000000000000 ]

[a,er,n]=IntergalMultiple(f,x,0,1,0,"2*x+3",0,"x+4*y+1",1e-12,64)//ѳʼΪ64,뱾سõ½,Է־ȷλࡣע1e-12ǿƵ
a =
[ 364.724206349206 ]
er =
[ 1.7053025658E-13 ]
n =
[ 128.000000000000    64.0000000000000    64.0000000000000 ]