关于积分的一个问题

在仰望星空

clc;clear
m=40;LSL=38.5;USL=41.5;
# n1 e+ N" A! i+ u# o+ N( ok1=0.6;k2=0.3;k3=1.0;sig0=0.5;
" l, m2 i' }+ ?" W" Ssig1=4.0;
0 n' E0 r* D: [0 C* m6 Z) t  l. fx=[m,sig0];
2 j- |" y) e$ Q. l. C
' j  W: U) J, n* s: d7 n

! F, x5 I8 j9 LEL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失


function g=myfun(y,x,m)
  i: c  u" @0 D+ c( ?$ ~3 [8 Tmu=x(1);sig=x(2);
f=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
; k0 y  X3 r: I& a7 ~- D  z  L4 B7 qg=(y-m).^2.*f;
; x# R: K( E$ `6 a4 h% T/ O! ~

* c" c0 M% v8 w! F. c# E  c??? Attempted to access y(13); index out of bounds because numel(y)=2.

Error in ==> quadl at 78
9 ^1 Z4 w5 Y; |" J2 O& `, a, e6 i! Bif ~isfinite(y(13))
7 [9 P& \; f- D" o0 Z
/ {% w2 V( [; g$ w5 Y: w
这个问题怎么解决啊  求助

hehainan86

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

hehainan86

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

在仰望星空

帖子怎么删除  这个发错了