我自己编了程序,能否请你帮我看看哪儿出错了
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
result