关于积分的一个问题

在仰望星空

clc;clear
m=40;LSL=38.5;USL=41.5;
# i+ k2 a6 [$ w/ O8 t1 f$ Sk1=0.6;k2=0.3;k3=1.0;sig0=0.5;
8 k; r" j0 m/ e7 M, w/ r% ?sig1=4.0;
% ?, Z- p/ n, j! Bx=[m,sig0];
0 d( @- \5 M& g% y3 A" w


EL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失
8 Z: X/ R! ?/ r7 D, m

$ e; u- ?' o4 vfunction g=myfun(y,x,m)
' t6 C; O+ E8 ~mu=x(1);sig=x(2);
8 E7 e  S& P* m! t7 x& b1 Sf=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
5 N5 J. O* G! ^% l* ]% Hg=(y-m).^2.*f;
; B4 Z3 c" l, }

??? Attempted to access y(13); index out of bounds because numel(y)=2.
6 s% O, j, _, I
Error in ==> quadl at 78
if ~isfinite(y(13))
5 n9 `* t; J6 ^- r9 @4 s. p: B

+ b2 x4 @# i# r! o& ?+ J6 }这个问题怎么解决啊  求助

在仰望星空

帖子怎么删除  这个发错了

hehainan86

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

hehainan86

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