一个maple的方程组

goalfly

我现在有三个方程
m1:=2*arctan(sqrt(b2/(1-b2)))/(s1*sqrt((r2-1)/(r1-1)*(1-b2)))-v1=0;
m2:=2*arctan(sqrt(b1/(1-b1)))/sqrt(1-b1)-v1=0;
m3:=2*arctan(sqrt((b-b2*(r2-1)/(r1-1))/(b1-b))*(b1*(r1-1)+1)/(b2*(r2-1)+1))/(s2*sqrt(b1-b))-v1=0;
然后我给出r1,r2,s1,s2的值，用solve可以解出b和v1的关系，但是当我输入一个v1的值时，我用solve解，计算机没有任何反应，请问是不是计算机没办法求出值呢。r1:=1.0678;r2:=6.5536;s1:=0.5;s2:=2; solve({m1,m2,m3},{b1,b2,b});我给v1定了5
  

jbzhou

可以得到数值解的。
  
m1:=2*arctan(sqrt(b2/(1-b2)))/(s1*sqrt((r2-1)/(r1-1)*(1-b2)))-v1=0;  
m2:=2*arctan(sqrt(b1/(1-b1)))/sqrt(1-b1)-v1=0;  
m3:=2*arctan(sqrt((b-b2*(r2-1)/(r1-1))/(b1-b))*(b1*(r1-1)+1)/(b2*(r2-1)+1))/(s2*sqrt(b1-b))-v1=0;  
r1:=1.0678;r2:=6.5536;s1:=0.5;s2:=2;
  
sol:=solve({m1,m2,m3},{b1,b2,b});这样得到的解为解析解，但是maple无法表示出，故用Rootof()表示。
  
此时给v1赋值  v1:=5; 而后执行evalf(sol);命令求出sol的近似解。
（例如，用你给的参数得到的解为：
  
  {b2 = 0.9767839216, b1 = -0.156039013 - 0.5350116824 I,
  
        b = -2.048436346 + 0.1945266145 I}

jbzhou

workspace见附件。