关于积分的一个问题

在仰望星空

clc;clear
6 Y2 g# P. s5 ?! |, @m=40;LSL=38.5;USL=41.5;
& x8 P1 X* l- N6 }5 p! bk1=0.6;k2=0.3;k3=1.0;sig0=0.5;
sig1=4.0;
2 C  k: g& _! D% nx=[m,sig0];

$ B. c' U- p/ o( g

EL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失
# k! _" L2 `# ]4 `
8 t+ [0 Y- ^3 ^/ j3 Q( a* s  N
function g=myfun(y,x,m)
mu=x(1);sig=x(2);
  X; m/ _% s5 C" `' `" [7 K% gf=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
g=(y-m).^2.*f;
, h7 c% ^* G* h" ~; s3 L% x( \

" ]! t2 }+ ?+ A  k5 e& n) y??? Attempted to access y(13); index out of bounds because numel(y)=2.

0 N1 h% i9 g- ~; z( [Error in ==> quadl at 78
1 M, H( }: R/ C& r" k- nif ~isfinite(y(13))
9 p$ Q9 K% N8 {+ y

这个问题怎么解决啊  求助

在仰望星空

帖子怎么删除  这个发错了

hehainan86

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

hehainan86

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