谁可以帮我看看程序哪里出了问题

Question

这个计算公式带指数我用matlab计算不出结果，
M2 = int(int((155.44+398.87*((x+y+3840-3828.342)/3828.342-0.002179)^0.154)(abs(-x-y+23.316)/sqrt(2)),y,-3.316,20),x,-3.316,20)+...
int(int((155.44+((x+y+3840-3828.342)/3828.342-0.002179)^0.154398.87)(abs(-x-y+23.316)/sqrt(2)),y,-11.658,-3.316),x,-3.316,5.026)+...
int(int((155.44+((x+3840-3828.342)/3828.342-0.002179)^0.154398.87+(y+3840-3828.342)/3828.34271320)(abs(-x-y+23.316)/sqrt(2)),y,-14.974-x,-11.658),x,-3.316,5.026)+...
int(int((155.44+((x+3840-3828.342)/3828.342-0.002179)^0.154398.87+(y+3840-3828.342)/3828.34271320)*(abs(-x-y+23.316)/sqrt(2)),y,-20,-3.316),x,5.026,20)
谁能帮我看一下啊

Answer 1

这个建议用integral2计算，类似这种

f=@(x,y,a)x.^2+a*y.^2; % xy是积分变量，a是其它参数
a = 1.0;
xmin=0,xmax=1,ymin=-1,ymax=1;%%xmin,xmax,ymin,ymax分别是x最小最大值、y最小最大值
integral2(@(x,y)f(x,y,a),xmin,xmax,ymin,ymax)