如何用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)&]