如何用mathematica得到方程组的实数根?

mathematica

真不知道这个符号解能干什么！

mathematica

(*清除所有变量*)
Clear["Global`*"];
(*求解方程,并得到10位有效数字*)
eq1=NSolve[{x^2+c^2==4^2,y^2+c^2==5^2,
            2/x==b/c,2/y==a/c,
            a+b==c},{a,b,c,x,y},10]
(*只保留数据,去掉变量符号*)
eq2={a,b,c,x,y}/.eq1
(*提取结果都是实数的结果*)
eq3=Cases[eq2,{_Real,_Real,_Real,_Real,_Real}]
(*提取都大于零的结果*)
eq4=Select[eq3,(#[[1]]>0&&#[[2]]>0&&#[[3]]>0&&#[[3]]>0&&#[[4]]>0&&#[[5]]>0)&]

wayne

4# mathematica

1. a=List@ToRules@Reduce[{x^2 + y^3 == 1, x^4 + y^4 == 2}, {x, y}, Reals]

复制代码

1. N[{x, y} //. a, 40]

复制代码

wayne

8# mathematica

1. aa=List@ToRules@Reduce[{x^2 + c^2 == 4^2, y^2 + c^2 == 5^2, 2/x == b/c, 2/y == a/c, a + b == c}, {a, b, c, x, y}, Reals]

复制代码

1. N[{a, b, c, x, y} //. aa, 20]

复制代码

wayne

如果是Mathematica 8的话，还简单些，直接用Solve

1. NSolve[{x^2 + c^2 == 4^2, y^2 + c^2 == 5^2, 2/x == b/c, 2/y == a/c, a + b == c}, {a, b, c, x, y}, Reals, WorkingPrecision -> 60]

复制代码

mathematica

如果是Mathematica 8的话，还简单些，直接用SolveNSolve[{x^2 + c^2 == 4^2, y^2 + c^2 == 5^2, 2/x == b/c, 2/y == a/c, a + b == c}, {a, b, c, x, y}, Reals, WorkingPrecision -> 60]
wayne 发表于 2011-2-26 23:21

我觉得mathematica有不如maple的地方，maple这个功能早就有了，
mathematica到了8才有这个功能，当然maple在群论这方面似乎也
比mathematica好。

mathematica

既然把我屏蔽了，我只好把我以前的代码重新发一下了：
(*清除所有变量*)
Clear["Global`*"];
(*求解方程,并得到10位有效数字*)
eq1=NSolve[{x^2+c^2==4^2,y^2+c^2==5^2,
            2/x==b/c,2/y==a/c,
            a+b==c},{a,b,c,x,y},10]
(*只保留数据,去掉变量符号*)
eq2={a,b,c,x,y}/.eq1
(*根根据虚部都等于零的结果*)
eq3=Select[eq2,Im[#]=={0,0,0,0,0}&]
(*提取结果都大于零的结果*)
eq4=Select[eq3,(#[[1]]>0&&#[[2]]>0&&#[[3]]>0&&#[[3]]>0&&#[[4]]>0&&#[[5]]>0)&]