关于积分的一个问题

在仰望星空

clc;clear
5 s4 t- F( f% V) g; y* r. |( {m=40;LSL=38.5;USL=41.5;
% `5 B4 B, t/ W  J5 {* g( ]/ n9 D3 bk1=0.6;k2=0.3;k3=1.0;sig0=0.5;
3 u3 \& D4 D  a0 I+ }9 wsig1=4.0;
x=[m,sig0];
; n* }& Z, A3 x# r& i; l8 w


EL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失
9 Q/ q$ q: {. @' X

, q3 X  ^6 P& w/ W4 Q) |function g=myfun(y,x,m)
) k6 \8 C8 ~4 M6 rmu=x(1);sig=x(2);
f=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
g=(y-m).^2.*f;
; G' s) c, S3 c

: Q3 g4 P6 O- Y1 A8 i; R4 e! }??? Attempted to access y(13); index out of bounds because numel(y)=2.
# }# Z9 X" F  o5 \1 L2 o
Error in ==> quadl at 78
if ~isfinite(y(13))
3 X2 I, b* a9 g

这个问题怎么解决啊  求助

hehainan86

错误提示的意思是：你访问的y已经超出了边界，y的维数为2。

hehainan86

你的函数myfun（x,y,m）在前面的一个式子出现两次里跟后面的定义的不一样，x与y调换了，是故意为之？？

在仰望星空

帖子怎么删除  这个发错了