- 注册时间
- 2008-11-26
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 149497
- 在线时间
- 小时
|
楼主 |
发表于 2017-7-3 16:45:01
|
显示全部楼层
- Clear["Global`*"];(*Clear all variables*)
- f=1/Sqrt[1+a]+1/Sqrt[1+b]+1/Sqrt[1+c]+x*(a*b*c-8)
- fa=D[f,a]
- fb=D[f,b]
- fc=D[f,c]
- fx=D[f,x]
- out=NSolve[{fa==0,fb==0,fc==0,fx==0},{a,b,c,x},200]
- cc=(f/.out)
- RootApproximant[cc]
复制代码
根据代码结果,只在实数区域内有a=b=c=2的驻点
{{a -> -1.0000000000000000000 - 1.7320508075688772935 I,
b -> -1.0000000000000000000 - 1.7320508075688772935 I,
c -> -1.0000000000000000000 - 1.7320508075688772935 I,
x -> 0.05296791160077942663 +
0.014192709138191943341 I}, {a -> -1.0000000000000000000 +
1.7320508075688772935 I,
b -> -1.0000000000000000000 + 1.7320508075688772935 I,
c -> -1.0000000000000000000 + 1.7320508075688772935 I,
x -> 0.05296791160077942663 - 0.014192709138191943341 I}, {a -> 2.,
b -> 2., c -> 2.0000, x -> 0.0240563}, {a -> 2., b -> 2.,
c -> 2.0000, x -> 0.0240563}, {a -> 2., b -> 2., c -> 2.0000,
x -> 0.0240563}, {a -> 2., b -> 2., c -> 2.0000, x -> 0.0240563}}
{1.611854897735312879 + 1.611854897735312879 I,
1.611854897735312879 -
1.611854897735312879 I, 1.73205, 1.73205, 1.73205, 1.73205}
{Root[27 + #1^4 &, 4],
Root[27 + #1^4 &, 3], Sqrt[3], Sqrt[3], Sqrt[3], Sqrt[3]}
|
|