关于积分的一个问题

在仰望星空

clc;clear
m=40;LSL=38.5;USL=41.5;
& h4 @0 A' a  o: A1 n8 G* p! [k1=0.6;k2=0.3;k3=1.0;sig0=0.5;
# o" i: K1 k# ?0 ?' `sig1=4.0;
x=[m,sig0];
. d" b# P" J( I; N4 L/ Z! J

) h/ }: Q" }! r. ~4 o' H" H! Y
/ I# Z1 B' s( i2 Z6 Y$ i& S6 x, c3 gEL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失
# K5 Q+ o( a" k! o

0 U, w6 u6 ^2 t) p3 `function g=myfun(y,x,m)
% |# [6 _) K) Xmu=x(1);sig=x(2);
9 b" i0 u; i6 n! xf=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
+ y. B4 j0 ~* U7 Sg=(y-m).^2.*f;


??? Attempted to access y(13); index out of bounds because numel(y)=2.
" C2 U! n; J1 v8 G
Error in ==> quadl at 78
if ~isfinite(y(13))
2 N3 v; Y' ?; e- E
" J4 C, k( j1 E6 r& E# a3 U! U
" X. \- C, y6 K这个问题怎么解决啊  求助

在仰望星空

帖子怎么删除  这个发错了

hehainan86

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

hehainan86

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