关于积分的一个问题

在仰望星空

clc;clear
m=40;LSL=38.5;USL=41.5;
$ W/ h8 ?8 O  O: w) d+ C7 s: Q) Sk1=0.6;k2=0.3;k3=1.0;sig0=0.5;
sig1=4.0;
x=[m,sig0];
% c! ], j4 i- P! Q- x
" A& x8 A6 C* ~. I6 B( F3 G, G

  F- O4 z. u. Y# H$ W5 g/ ^# R: u/ BEL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失
/ P+ ~& }, [: s: w6 R
9 W  P/ `/ v  f- Y: f  S
function g=myfun(y,x,m)
mu=x(1);sig=x(2);
" n, K- L1 A5 d" nf=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
g=(y-m).^2.*f;
4 d" O: {0 Z3 t# T5 N4 p, a) d/ ~& s
/ h& _; @( c" K6 C; X% G
' a( M8 e+ q, b7 y6 Y??? Attempted to access y(13); index out of bounds because numel(y)=2.
$ s6 e1 V( {& L. ]
Error in ==> quadl at 78
$ r  P( M" ~1 P% J. x( d) q1 Eif ~isfinite(y(13))
  d$ N7 [$ [9 k" S* U' m7 F7 S

- Q2 O6 F8 u/ R& x- L这个问题怎么解决啊  求助

在仰望星空

帖子怎么删除  这个发错了

hehainan86

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

hehainan86

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