件浮点数和任意精度的浮点数甚至于矩阵的计算等等。总之, Maple 可以进行任意数值计算。
但是, 任何软件或程序毕竟只是人们进行科学研究的一种必要的辅助, 即便它有很多优点,
但也有它的局限性, 为了客观地认识数学软件、认识 Maple, 下面通过两个简单例子予以说明。
第一个简单的数值计算实例想说明 Maple 数值计算的答案的正确性:
>
3!!!;
2601218943565795100204903227081043611191521875016945785727541837850835631156947382240678577958130457082619920575892
2472595366415651620520158737919845877408325291052446903888118841237643411919510455053466586162432719401971139098455
3672727853709934562985558671936977407000370043078375899742067678401696720784628062922903210716166986726054898844551
4257193985499448939594496064045132362140265986193073249369770477606067680670176491669403034819961881455625195592566
9188308255149429475965372748456246288242345265977897377408964665539924359287862125159674832209760295056966999272846
7056374713753301924831358707612541268341586012944756601145542074958995256354306828863463108496565068277155299625679
0845235702552186222358130016700834523443236821935793184701956510729781804354173890560727428048583995919729021726612
2912984205160675790362323376994539641914751755675576953922338030568253085999774416757843528159134613403946049012695
4202883834710136373382448450666009334848444071193129253769465735433737572477223018153403264717753198453734147867432
7048457983786618703257405938924215709695994630557521063203263493209220738320923356309923267504401701760572026010829
2880423356066430898887102973807975780130560495763428386830571906622052911748225105366977566030295740433879834715185
5260280533386635713910104633641976909739743228599421983704697910995630338960467588986579571117656667003915674815311
5943980043625399399731203066490601325311304719028898491856203766669164468791125249193754425845895000311561682974304
6411425380748972817233759553806617198014046779356147936352662656833395097600000000000000000000000000000000000000000
0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000
0000000000000000000000
上述运算结果在 IBM PC机(1G, 128M)上计算只需要 0.01秒, 得到如此复杂的结果(1747位),
一个自然的问题是: 答案正确吗?
为了回答这个问题, 我们借助于数值分析方法, 由 Stiring 公式
)exp(2! nnnn
n
−⋅⋅≈
π
可得:
, 前三位数字与 Maple 输出结果相同, 且两者结果均为 1747 位。
另外, 在 720!的计算中, 5 的因子的个数为:
1746
1060091.2!720 ×≈
178
5
720
5
720
5
720
5
720
432
=
⎥
⎦
⎤
⎢
⎣
⎡
+
⎥
⎦
⎤
⎢
⎣
⎡
+
⎥
⎦
⎤
⎢
⎣
⎡
+
⎥
⎦
⎤
⎢
⎣
⎡
这些 5 与足够多的 2 相乘将得到 178 个 0, 而 Maple 的输出结果中最后 178 位数为零。由此, 可
以相信 Maple 结果的正确性。
另一个例子则想说明 Maple 计算的局限性:
() ()
?8 ?8
6/23/1
=−=−
Maple 在处理问题时, 为了避免失根, 从不求算术式的近似值, 分数则化简为既约分数。因
此, 在 Maple 中很容易得到:
()
)
6/23/1
8 8 −=−
显然这是错误的。这一点可以从代数的角度予以分析。
不妨设
, 则 , 即 , 显然
()
x=−
3/1
8
08
3
=+x
0)42)(2(
2
=+−+ xxx
)
3/1
8−
有 3 个结果,
-2 是其实数结果。
- 4 -
评论0
最新资源