m=10;n=5;
for L=1:30
myfun=@(u,L) quadgk(@(v) chi2pdf(v,(m*(n-1)))./((1-chi2cdf(v.*u./(m*(n-1)),n-1))*(1-chi2cdf(v.*u./(m*(n-1)),n-1).^L)),0,inf)-370.34
fun=@(u) myfun(u,L)
a0=[0+eps,10000+eps]
result=fzero(fun,a0)
result(L)=u;
end
、
大家帮我看看程序在哪儿出错了，谢谢！