关于积分的一个问题

在仰望星空

clc;clear
m=40;LSL=38.5;USL=41.5;
% j  B% m: h- {- B2 s$ hk1=0.6;k2=0.3;k3=1.0;sig0=0.5;
1 P) \: L# ?. |sig1=4.0;
. A  c( G# I- w0 Xx=[m,sig0];
. A- k+ Q9 J7 k! l1 G3 [' y/ ~
- b8 [: B8 ^' I5 p9 w! C
1 t  _2 `8 d, f% o6 ~# `
EL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失
* l! j/ t9 p( t5 P1 _

6 o# V' M- x, ~function g=myfun(y,x,m)
mu=x(1);sig=x(2);
f=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
, o' ^) e5 N5 bg=(y-m).^2.*f;
, r9 U1 B4 F5 U# I* f9 c
  g- O8 T. E+ G; w. T
; D% @# I5 H" |2 p- S??? Attempted to access y(13); index out of bounds because numel(y)=2.

Error in ==> quadl at 78
if ~isfinite(y(13))

& U: h" e& ~. K/ `7 r
这个问题怎么解决啊  求助

在仰望星空

帖子怎么删除  这个发错了

hehainan86

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

hehainan86

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