foozhencheng 发表于 2016-5-1 10:29 
这应该就只能数值了。。。
我自己编了程序，能否请你帮我看看哪儿出错了
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
扫码加好友，拉您进群
全部版块
我的主页
收藏
![]GB37)3Z{LPY3N_A2RR0V{2.png](https://bbs-pic.datacourse.cn/forum/201604/29/1859533ed2xsujuuzutu3d.png "]GB37)3Z{LPY3N_A2RR0V{2.png")
