目前状况是
宣告两变量后
求联立方程式
得到的两解ansm及ansv
求位在0.95积率的beta积分
但每次执行都会卡住
http://imgur.com/m7TPfIx
也找不到是什么原因
syms a b
m=1.002139;
v=1.004282;
[eqm,eqv]=solve('m=a/(a+b)','v=a*b/((a+b)^2*(a+b+1))');
ansm=subs(eqm,{a,b},{m,v});
ansv=subs(eqv,{a,b},{m,v});
prob_6=betainv(0.95,ansm,ansv)
p.s其中的m,v是我在前面计算完后所求得的值
所以才会设这两个值