fun1=0.23(
α
1
+
β
1
)-
2
b
(
α
1
+
β
1
)
1+b
+
0.83
2
b
(
α
1
+
β
1
)
1.21+b
+0.25(
α
2
+
β
2
)-
0.83
2
b
(
α
2
+
β
2
)
1.210+b
+
0.69
2
b
(
α
2
+
β
2
)
1.44+b
fun2=0.12
α
1
+0.13
α
2
+
b(
α
1
-
β
1
)
2(1+b)
-
b(
α
1
-
β
1
)
2(1.21+b)
-Log
1
1.`+b
(
α
1
-
β
1
)+Log
1.1
1.21+b
(
α
1
-
β
1
)+
b(
α
2
-
β
2
)
2(1.21+b)
-
b(
α
2
-
β
2
)
2(1.44`+b)
-Log
1.1
1.21+b
(
α
2
-
β
2
)+Log
1.2
1.44+b
(
α
2
-
β
2
)
b=FindRoot[fun2==0,b]