MathCAD和Maple解方程的问题！

stargazer

用MathCAD解方程，调用root函数，发现设不同初始值，得到结果差异很大，到底哪一个才对阿，抓狂！！！

令X0=5000,解得 x=5054；
令X0=6000,解得 x=6420；
令X0=8000,解得 x=3833；
令X0=10000,解得 x=6420；
令X0=12000,解得 x=5054；

用Maple解的结果，等于5054
但是把方程中的0.0072换成0.0064，发现结果出不来了

                               
登录/注册后可看大图

[ 本帖最后由 stargazer 于 2007-1-22 23:49 编辑 ]

stargazer

用Maple解的结果，等于5054
但是把方程中的0.0072换成0.0064，发现结果出不来了

                               
登录/注册后可看大图

stargazer

算的结果是疲劳寿命，但是结果这么不一致，很头疼，方程到底应该怎么解？？
大家帮帮忙

dingd

应该是精度误差判断引起的，如下：
x               f(x)
2527        5.5374201E-6
3285.3        1.8366879E-6
4043.6        4.3810763E-7
4801.9        2.0968216E-8
5560.2        6.7445897E-8
6318.5        3.4482538E-7
7076.8        7.3853063E-7
7835.1        1.1884995E-6
8593.4        1.6618859E-6
9351.7        2.1402621E-6
10110        2.6131922E-6

f(x)的变化很小，如终止判断误差取为1E-6左右，就可的任意多的解。实际上只有一个最优解。下面是用1stOpt算的：
1：Function 0.007634063*x^(-0.10448)+0.685*x^(-0.6)-0.007238722;
结果：x: 5053.95687520354

2：Function 0.007634063*x^(-0.10448)+0.685*x^(-0.6)-0.006415;
结果：x: 6969.7310450444

stargazer

确实如此，本来函数值就小，默认精度0.001根本不够，改成1e-9后差别就很小了
谢谢楼上

stargazer

但是Maple为什么只列出方程，而不能计算结果呢？

suffer

fsolve(0.76341e-2/n^.10448+.575/n^.6 = 0.6415e-2, n, complex)

5419.531242