关于积分的一个问题

在仰望星空

clc;clear
% [7 i2 g3 @3 O+ _* Lm=40;LSL=38.5;USL=41.5;
/ _  D2 j0 p' i/ M" ^% nk1=0.6;k2=0.3;k3=1.0;sig0=0.5;
sig1=4.0;
' o2 r$ a5 y/ t5 L9 v% X& w9 Xx=[m,sig0];



8 a9 d' U& |3 B1 q) k! @& N! jEL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失
) m, v( a: ^1 o* s6 W' g  {% T" o
1 T! I% N# p$ B2 Y. B9 b$ D
' |+ {, j2 R' D: a2 U8 n) ofunction g=myfun(y,x,m)
mu=x(1);sig=x(2);
f=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
# N, y: N  j/ q$ ]3 }g=(y-m).^2.*f;


7 \0 g+ _  d* M. L7 c' ^??? Attempted to access y(13); index out of bounds because numel(y)=2.

Error in ==> quadl at 78
# {: E8 K# s7 R/ F& x: hif ~isfinite(y(13))


这个问题怎么解决啊  求助

在仰望星空

帖子怎么删除  这个发错了

hehainan86

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

hehainan86

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