目前狀況是
宣告兩變數後
求聯立方程式
得到的兩解ansm及ansv
求位在0.95積率的beta積分
但每次執行都會卡住
也找不到是什麼原因
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是我在前面計算完後所求得的值
所以才會設這兩個值