各位大侠，帮忙看看啊
r=rand(22,1);<BR>xma=;<BR>xmi=;<BR>for i=1:22<BR> x(i)=xmi(i)+xma(i)*unifrnd(0,1);<BR> syms t real;<BR> a(i)=gamma(r(i)/2);<BR> b(i)=t-xmi(i);<BR> f(i)=(b(i)^(r(i)/2-1)*exp(-b(i)/2))/(a(i)*2^(r(i)/2));<BR> M=int(f(i),t,xmi(i),x(i));<BR> m(i)=double(M(i));<BR> end<BR><BR>为什么m（i）可以循环求取几个后，就出现问题？<BR>如果是程序有错误应该一次循环也不行啊？不懂啊。<BR>请知道的大侠帮帮阿<BR><BR>该怎么改？回复：（ccchen）各位大侠，帮忙看看啊
<DIV class=quote><B>以下是引用<I>ccchen</I>在2006-4-19 15:52:56的发言：</B><BR>r=rand(22,1);<BR>xma=;<BR>xmi=;<BR>for i=1:22<BR> x(i)=xmi(i)+xma(i)*unifrnd(0,1);<BR> syms t real;<BR> a(i)=gamma(r(i)/2);<BR> <FONT color=#ff0000>b(i)=t-xmi(i);</FONT><BR> f(i)=(b(i)^(r(i)/2-1)*exp(-b(i)/2))/(a(i)*2^(r(i)/2));<BR> M=int(f(i),t,xmi(i),x(i));<BR> m(i)=double(M(i));<BR>end<BR><BR>为什么m（i）可以循环求取几个后，就出现问题？<BR>如果是程序有错误应该一次循环也不行啊？不懂啊。<BR>请知道的大侠帮帮阿<BR><BR>该怎么改？</DIV><P>t是符号型变量，xmi是数值型变量，当然不能做运算了</P>
页:
[1]