应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Mathematica 简单 四元一次方程组 ，求解问题~
101/1返回列表
查看: 2317  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

chineserice

禁虫 (正式写手)
本帖内容被屏蔽

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-06-14 00:16:11
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mshwangg

至尊木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
chineserice: 金币+10, ★★★★★最佳答案, 感谢，我在研究一下下~ 2012-06-14 01:47:46
chineserice: 金币+20, ★★★★★最佳答案, 我也搞定了~ 非常感谢了，就是不知道怎么和老板汇报，在国外也不能说我用了盗版软件吧。 2012-06-14 13:21:00
很快就出结果了呀，你在未知量{A,B,C,D}前面的逗号之前少了一个}
就是方程需要用{}括起来，你只写了{，没写}
结果如下：
{{A -> 0. + (1. (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
               1. (0.5 xi + xim1)^3) -
            1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
               1. (0.5 xi + xim1)^3)) (-1. (1.09091 +
               0. (0.5 xi + xim1)) (-0.5 xi^2 +
               1. (0.5 xi + xim1)^2) + (0. xi -
               1. xim1) (0. (0.5 xi + xim1)^2 +
               2. (0.0454545 (-0.5 xi - 0.5 xim1) +
                  0.0454545 (0.5 xi + 0.5 xip1)))) + (-0.25 xi^3 +
            1. xi^2 xim1 + 1. xi xim1^2 +
            0. xim1^3) (-1. (1.09091 +
               0. (0.5 xi + xim1)) (0.5 xi^3 -
               1. (0.5 xi + xim1)^3) + (0. xi -
               1. xim1) (0. (0.5 xi + xim1)^3 +
               3. (0.0454545 (0.5 xi + 0.5 xim1)^2 +
                  0.0454545 (0.5 xi + 0.5 xip1)^2)))))/(-1. ((0. xi -
             1. xim1) (-0.5 xi^3 - 1. (0.5 xi + xim1)^3) -
          1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
             1. (0.5 xi + xim1)^3)) (-1. (-0.5 xi^2 +
             1. (0.5 xi + xim1)^2) (-1. (0.5 xi + xim1) -
             1. (0.5 xi + xip1)) + (0. xi -
             1. xim1) (1. (0.5 xi + xim1)^2 -
             1. (0.5 xi + xip1)^2)) + (-0.25 xi^3 + 1. xi^2 xim1 +
          1. xi xim1^2 +
         
          0. xim1^3) (-1. (0.5 xi^3 -
             1. (0.5 xi + xim1)^3) (-1. (0.5 xi + xim1) -
             1. (0.5 xi + xip1)) + (0. xi -
             1. xim1) (-1. (0.5 xi + xim1)^3 -
             1. (0.5 xi + xip1)^3))) +
    1. ((-1. (1.09091 + 0. (0.5 xi + xim1)) (-0.5 xi^2 +
             1. (0.5 xi + xim1)^2) + (0. xi -
             1. xim1) (0. (0.5 xi + xim1)^2 +
             2. (0.0454545 (-0.5 xi - 0.5 xim1) +
                0.0454545 (0.5 xi + 0.5 xip1))))/(-0.25 xi^3 +
          1. xi^2 xim1 + 1. xi xim1^2 +
          0. xim1^3) - (1. (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
                  1. (0.5 xi + xim1)^3) -
               1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
                  1. (0.5 xi + xim1)^3)) (-1. (1.09091 +
                  0. (0.5 xi + xim1)) (-0.5 xi^2 +
                  1. (0.5 xi + xim1)^2) + (0. xi -
                  1. xim1) (0. (0.5 xi + xim1)^2 +
                  2. (0.0454545 (-0.5 xi - 0.5 xim1) +
                    0.0454545 (0.5 xi + 0.5 xip1)))) + (-0.25 xi^3 +
               1. xi^2 xim1 + 1. xi xim1^2 +
               0. xim1^3) (-1. (1.09091 +
                  0. (0.5 xi + xim1)) (0.5 xi^3 -
                  1. (0.5 xi + xim1)^3) + (0. xi -
                  1. xim1) (0. (0.5 xi + xim1)^3 +
                  3. (0.0454545 (0.5 xi + 0.5 xim1)^2 +
                    
                    0.0454545 (0.5 xi +
                    0.5 xip1)^2)))) (-1. (-0.5 xi^2 +
               1. (0.5 xi + xim1)^2) (-1. (0.5 xi + xim1) -
               1. (0.5 xi + xip1)) + (0. xi -
               1. xim1) (1. (0.5 xi + xim1)^2 -
               1. (0.5 xi + xip1)^2)))/((-0.25 xi^3 + 1. xi^2 xim1 +
            1. xi xim1^2 +
            0. xim1^3) (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
                  1. (0.5 xi + xim1)^3) -
               1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
                  1. (0.5 xi + xim1)^3)) (-1. (-0.5 xi^2 +
                  1. (0.5 xi + xim1)^2) (-1. (0.5 xi + xim1) -
                  1. (0.5 xi + xip1)) + (0. xi -
                  1. xim1) (1. (0.5 xi + xim1)^2 -
                  1. (0.5 xi + xip1)^2)) + (-0.25 xi^3 +
               1. xi^2 xim1 + 1. xi xim1^2 +
               0. xim1^3) (-1. (0.5 xi^3 -
                  1. (0.5 xi + xim1)^3) (-1. (0.5 xi + xim1) -
                  1. (0.5 xi + xip1)) + (0. xi -
                  1. xim1) (-1. (0.5 xi + xim1)^3 -
                  1. (0.5 xi + xip1)^3))))) +
    1. ((1.09091 + 0. (0.5 xi + xim1))/(
       0. xi - 1. xim1) - (1. (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
                  1. (0.5 xi + xim1)^3) -
               1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
                  1. (0.5 xi + xim1)^3)) (-1. (1.09091 +
                  0. (0.5 xi + xim1)) (-0.5 xi^2 +
                  1. (0.5 xi + xim1)^2) + (0. xi -
                  1. xim1) (0. (0.5 xi + xim1)^2 +
                  2. (0.0454545 (-0.5 xi - 0.5 xim1) +
                    0.0454545 (0.5 xi + 0.5 xip1)))) + (-0.25 xi^3 +
               1. xi^2 xim1 + 1. xi xim1^2 +
               0. xim1^3) (-1. (1.09091 +
                  0. (0.5 xi + xim1)) (0.5 xi^3 -
                  1. (0.5 xi + xim1)^3) + (0. xi -
                  1. xim1) (0. (0.5 xi + xim1)^3 +
                  3. (0.0454545 (0.5 xi + 0.5 xim1)^2 +
                    0.0454545 (0.5 xi + 0.5 xip1)^2)))) (-1. (0.5 xi +
                xim1) - 1. (0.5 xi + xip1)))/((0. xi -
            1. xim1) (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
                  1. (0.5 xi + xim1)^3) -
               1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
                  1. (0.5 xi + xim1)^3)) (-1. (-0.5 xi^2 +
                  1. (0.5 xi + xim1)^2) (-1. (0.5 xi + xim1) -
                  1. (0.5 xi + xip1)) + (0. xi -
                  1. xim1) (1. (0.5 xi + xim1)^2 -
                  1. (0.5 xi + xip1)^2)) + (-0.25 xi^3 +
               1. xi^2 xim1 + 1. xi xim1^2 +
               0. xim1^3) (-1. (0.5 xi^3 -
                  1. (0.5 xi + xim1)^3) (-1. (0.5 xi + xim1) -
                  1. (0.5 xi + xip1)) + (0. xi -
                  1. xim1) (-1. (0.5 xi + xim1)^3 -
                  1. (0.5 xi + xip1)^3)))) -
       1/(0. xi - 1. xim1)
         1. (-0.5 xi -
          1. (0.5 xi +
             xim1)) ((-1. (1.09091 + 0. (0.5 xi + xim1)) (-0.5 xi^2 +
                1. (0.5 xi + xim1)^2) + (0. xi -
                1. xim1) (0. (0.5 xi + xim1)^2 +
                2. (0.0454545 (-0.5 xi - 0.5 xim1) +
                   0.0454545 (0.5 xi + 0.5 xip1))))/(-0.25 xi^3 +
             1. xi^2 xim1 + 1. xi xim1^2 +
             0. xim1^3) - (1. (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
                    1. (0.5 xi + xim1)^3) -
                  1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
                    1. (0.5 xi + xim1)^3)) (-1. (1.09091 +
                    0. (0.5 xi + xim1)) (-0.5 xi^2 +
                    1. (0.5 xi + xim1)^2) + (0. xi -
                    1. xim1) (0. (0.5 xi + xim1)^2 +
                    2. (0.0454545 (-0.5 xi - 0.5 xim1) +
                    0.0454545 (0.5 xi + 0.5 xip1)))) + (-0.25 xi^3 +
                  1. xi^2 xim1 + 1. xi xim1^2 +
                  0. xim1^3) (-1. (1.09091 +
                    0. (0.5 xi + xim1)) (0.5 xi^3 -
                    1. (0.5 xi + xim1)^3) + (0. xi -
                    1. xim1) (0. (0.5 xi + xim1)^3 +
                    3. (0.0454545 (0.5 xi + 0.5 xim1)^2 +
                    
                    0.0454545 (0.5 xi +
                    0.5 xip1)^2)))) (-1. (-0.5 xi^2 +
                  1. (0.5 xi + xim1)^2) (-1. (0.5 xi + xim1) -
                  1. (0.5 xi + xip1)) + (0. xi -
                  1. xim1) (1. (0.5 xi + xim1)^2 -
                  1. (0.5 xi + xip1)^2)))/((-0.25 xi^3 +
               1. xi^2 xim1 + 1. xi xim1^2 +
               0. xim1^3) (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
                    1. (0.5 xi + xim1)^3) -
                  1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
                    1. (0.5 xi + xim1)^3)) (-1. (-0.5 xi^2 +
                    1. (0.5 xi + xim1)^2) (-1. (0.5 xi + xim1) -
                    1. (0.5 xi + xip1)) + (0. xi -
                    1. xim1) (1. (0.5 xi + xim1)^2 -
                    1. (0.5 xi + xip1)^2)) + (-0.25 xi^3 +
                  1. xi^2 xim1 + 1. xi xim1^2 +
                  0. xim1^3) (-1. (0.5 xi^3 -
                    1. (0.5 xi + xim1)^3) (-1. (0.5 xi + xim1) -
                    1. (0.5 xi + xip1)) + (0. xi -
                    1. xim1) (-1. (0.5 xi + xim1)^3 -
                    1. (0.5 xi + xip1)^3)))))),
  B -> -((1. (1.09091 + 0. (0.5 xi + xim1)))/(
     0. xi - 1. xim1)) + (1. (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
               1. (0.5 xi + xim1)^3) -
            1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
               1. (0.5 xi + xim1)^3)) (-1. (1.09091 +
               0. (0.5 xi + xim1)) (-0.5 xi^2 +
               1. (0.5 xi + xim1)^2) + (0. xi -
               1. xim1) (0. (0.5 xi + xim1)^2 +
               2. (0.0454545 (-0.5 xi - 0.5 xim1) +
                  0.0454545 (0.5 xi + 0.5 xip1)))) + (-0.25 xi^3 +
            1. xi^2 xim1 + 1. xi xim1^2 +
            0. xim1^3) (-1. (1.09091 +
               0. (0.5 xi + xim1)) (0.5 xi^3 -
               1. (0.5 xi + xim1)^3) + (0. xi -
               1. xim1) (0. (0.5 xi + xim1)^3 +
               3. (0.0454545 (0.5 xi + 0.5 xim1)^2 +
                  0.0454545 (0.5 xi + 0.5 xip1)^2)))) (-1. (0.5 xi +
            xim1) - 1. (0.5 xi + xip1)))/((0. xi -
         1. xim1) (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
               1. (0.5 xi + xim1)^3) -
            1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
               1. (0.5 xi + xim1)^3)) (-1. (-0.5 xi^2 +
               1. (0.5 xi + xim1)^2) (-1. (0.5 xi + xim1) -
               1. (0.5 xi + xip1)) + (0. xi -
               1. xim1) (1. (0.5 xi + xim1)^2 -
               1. (0.5 xi + xip1)^2)) + (-0.25 xi^3 + 1. xi^2 xim1 +
            1. xi xim1^2 +
            0. xim1^3) (-1. (0.5 xi^3 -
               1. (0.5 xi + xim1)^3) (-1. (0.5 xi + xim1) -
               1. (0.5 xi + xip1)) + (0. xi -
               1. xim1) (-1. (0.5 xi + xim1)^3 -
               1. (0.5 xi + xip1)^3)))) +
    1/(0. xi - 1. xim1)
      1. (-0.5 xi -
       1. (0.5 xi +
          xim1)) ((-1. (1.09091 + 0. (0.5 xi + xim1)) (-0.5 xi^2 +
             1. (0.5 xi + xim1)^2) + (0. xi -
             1. xim1) (0. (0.5 xi + xim1)^2 +
             2. (0.0454545 (-0.5 xi - 0.5 xim1) +
                0.0454545 (0.5 xi + 0.5 xip1))))/(-0.25 xi^3 +
          1. xi^2 xim1 + 1. xi xim1^2 +
          0. xim1^3) - (1. (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
                  1. (0.5 xi + xim1)^3) -
               1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
                  1. (0.5 xi + xim1)^3)) (-1. (1.09091 +
                  0. (0.5 xi + xim1)) (-0.5 xi^2 +
                  1. (0.5 xi + xim1)^2) + (0. xi -
                  1. xim1) (0. (0.5 xi + xim1)^2 +
                  2. (0.0454545 (-0.5 xi - 0.5 xim1) +
                    0.0454545 (0.5 xi + 0.5 xip1)))) + (-0.25 xi^3 +
               1. xi^2 xim1 + 1. xi xim1^2 +
               0. xim1^3) (-1. (1.09091 +
                  0. (0.5 xi + xim1)) (0.5 xi^3 -
                  1. (0.5 xi + xim1)^3) + (0. xi -
                  1. xim1) (0. (0.5 xi + xim1)^3 +
                  3. (0.0454545 (0.5 xi + 0.5 xim1)^2 +
                    
                    0.0454545 (0.5 xi +
                    0.5 xip1)^2)))) (-1. (-0.5 xi^2 +
               1. (0.5 xi + xim1)^2) (-1. (0.5 xi + xim1) -
               1. (0.5 xi + xip1)) + (0. xi -
               1. xim1) (1. (0.5 xi + xim1)^2 -
               1. (0.5 xi + xip1)^2)))/((-0.25 xi^3 + 1. xi^2 xim1 +
            1. xi xim1^2 +
            0. xim1^3) (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
                  1. (0.5 xi + xim1)^3) -
               1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
                  1. (0.5 xi + xim1)^3)) (-1. (-0.5 xi^2 +
                  1. (0.5 xi + xim1)^2) (-1. (0.5 xi + xim1) -
                  1. (0.5 xi + xip1)) + (0. xi -
                  1. xim1) (1. (0.5 xi + xim1)^2 -
                  1. (0.5 xi + xip1)^2)) + (-0.25 xi^3 +
               1. xi^2 xim1 + 1. xi xim1^2 +
               0. xim1^3) (-1. (0.5 xi^3 -
                  1. (0.5 xi + xim1)^3) (-1. (0.5 xi + xim1) -
                  1. (0.5 xi + xip1)) + (0. xi -
                  1. xim1) (-1. (0.5 xi + xim1)^3 -
                  1. (0.5 xi + xip1)^3))))),
  C -> -(1. (-1. (1.09091 + 0. (0.5 xi + xim1)) (-0.5 xi^2 +
             1. (0.5 xi + xim1)^2) + (0. xi -
             1. xim1) (0. (0.5 xi + xim1)^2 +
             2. (0.0454545 (-0.5 xi - 0.5 xim1) +
                0.0454545 (0.5 xi + 0.5 xip1)))))/(-0.25 xi^3 +
       1. xi^2 xim1 + 1. xi xim1^2 +
       0. xim1^3) + (1. (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
               1. (0.5 xi + xim1)^3) -
            1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
               1. (0.5 xi + xim1)^3)) (-1. (1.09091 +
               0. (0.5 xi + xim1)) (-0.5 xi^2 +
               1. (0.5 xi + xim1)^2) + (0. xi -
               1. xim1) (0. (0.5 xi + xim1)^2 +
               2. (0.0454545 (-0.5 xi - 0.5 xim1) +
                  0.0454545 (0.5 xi + 0.5 xip1)))) + (-0.25 xi^3 +
            1. xi^2 xim1 + 1. xi xim1^2 +
            0. xim1^3) (-1. (1.09091 +
               0. (0.5 xi + xim1)) (0.5 xi^3 -
               1. (0.5 xi + xim1)^3) + (0. xi -
               1. xim1) (0. (0.5 xi + xim1)^3 +
               3. (0.0454545 (0.5 xi + 0.5 xim1)^2 +
                  0.0454545 (0.5 xi +
                    0.5 xip1)^2)))) (-1. (-0.5 xi^2 +
            1. (0.5 xi + xim1)^2) (-1. (0.5 xi + xim1) -
            1. (0.5 xi + xip1)) + (0. xi -
            1. xim1) (1. (0.5 xi + xim1)^2 -
            1. (0.5 xi + xip1)^2)))/((-0.25 xi^3 + 1. xi^2 xim1 +
         1. xi xim1^2 +
         0. xim1^3) (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
               1. (0.5 xi + xim1)^3) -
            1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
               1. (0.5 xi + xim1)^3)) (-1. (-0.5 xi^2 +
               1. (0.5 xi + xim1)^2) (-1. (0.5 xi + xim1) -
               1. (0.5 xi + xip1)) + (0. xi -
               1. xim1) (1. (0.5 xi + xim1)^2 -
               1. (0.5 xi + xip1)^2)) + (-0.25 xi^3 + 1. xi^2 xim1 +
            1. xi xim1^2 +
            0. xim1^3) (-1. (0.5 xi^3 -
               1. (0.5 xi + xim1)^3) (-1. (0.5 xi + xim1) -
               1. (0.5 xi + xip1)) + (0. xi -
               1. xim1) (-1. (0.5 xi + xim1)^3 -
               1. (0.5 xi + xip1)^3)))),
  D -> -(1. (-1. ((0. xi - 1. xim1) (-0.5 xi^3 -
               1. (0.5 xi + xim1)^3) -
            1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
               1. (0.5 xi + xim1)^3)) (-1. (1.09091 +
               0. (0.5 xi + xim1)) (-0.5 xi^2 +
               1. (0.5 xi + xim1)^2) + (0. xi -
               1. xim1) (0. (0.5 xi + xim1)^2 +
               2. (0.0454545 (-0.5 xi - 0.5 xim1) +
                  0.0454545 (0.5 xi + 0.5 xip1)))) + (-0.25 xi^3 +
            1. xi^2 xim1 + 1. xi xim1^2 +
            0. xim1^3) (-1. (1.09091 +
               0. (0.5 xi + xim1)) (0.5 xi^3 -
               1. (0.5 xi + xim1)^3) + (0. xi -
               1. xim1) (0. (0.5 xi + xim1)^3 +
               3. (0.0454545 (0.5 xi + 0.5 xim1)^2 +
                  0.0454545 (0.5 xi + 0.5 xip1)^2)))))/(-1. ((0. xi -
            1. xim1) (-0.5 xi^3 - 1. (0.5 xi + xim1)^3) -
         1. (-0.5 xi - 1. (0.5 xi + xim1)) (0.5 xi^3 -
            1. (0.5 xi + xim1)^3)) (-1. (-0.5 xi^2 +
            1. (0.5 xi + xim1)^2) (-1. (0.5 xi + xim1) -
            1. (0.5 xi + xip1)) + (0. xi -
            1. xim1) (1. (0.5 xi + xim1)^2 -
            1. (0.5 xi + xip1)^2)) + (-0.25 xi^3 + 1. xi^2 xim1 +
         1. xi xim1^2 +
         0. xim1^3) (-1. (0.5 xi^3 -
            1. (0.5 xi + xim1)^3) (-1. (0.5 xi + xim1) -
            1. (0.5 xi + xip1)) + (0. xi -
            1. xim1) (-1. (0.5 xi + xim1)^3 -
            1. (0.5 xi + xip1)^3)))}}
赞一下 回复此楼
2楼2012-06-14 00:45:27
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

chineserice

禁虫 (正式写手)
本帖内容被屏蔽
3楼2012-06-14 01:56:53
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mshwangg

至尊木虫 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by chineserice at 2012-06-14 01:56:53
奇了怪了~ 当时我检查没有符号颜色提示。
不好意思，还有不少问题，我要变成分数形式怎么调整，因为这个是超级计算机计算用，对精度要求很高，一般都是双精度~
还有就是0. + (1. (-1. ((0. xi - 1. xim1) (-0.5  ...

用Simplify或者FullSimplify试试看能不能获得你想要的结果
如果你在其他计算机上也用mathematica的话，不必考虑数据精度问题，mathematica中也没有数据类型啊，除非你用fortran。
用fortran的话还真得琢磨下分子分母的数量级，否则舍入误差是个问题。要是没有太大问题的话，把表达式直接用FortranForm[expr]，再输出文本，copy进fortran就差不多可以了，省得再type一遍。
0.就当作是0就可以了，这是mathematica的特殊癖好之一。
输出用Write呀，比如说要输出y1,y2随t的变化，那就Write[strm,t,y1,y2]，再对这个语句作Do循环不就行了。不知道是不是准确理解了你的问题
赞一下 回复此楼
4楼2012-06-14 11:13:34
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zyt2001
5楼2012-06-14 11:56:13
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

chineserice

禁虫 (正式写手)
本帖内容被屏蔽
6楼2012-06-14 12:00:10
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mshwangg

至尊木虫 (正式写手)
引用回帖:
6楼: Originally posted by chineserice at 2012-06-14 12:00:10
TXS 再试验一下，因为我还是要转到fortran上面计算，最好精确到小数点后15位，双精度要求~...

这是FortranForm[A]的结果，数字部分还是满足双精度要求的
0. + (0.*xi**7 + xi**6*(0. + 0.*xim1) + 0.*xim1**6 + 0.*xim1**7 +
     xi**5*(-0.2556818181818182*xim1 + 0.*xim1**2 + 0.*xip1) +
     0.*xim1**5*xip1 + 0.*xim1**4*xip1**2 +
     xi**4*(1.0227272727272727*xim1**2 + 0.*xim1**3 +
        0.03409090909090909*xim1*xip1 + 0.*xip1**2) +
     xi**3*(0.9517045454545452*xim1**3 + 0.*xim1**4 -
        0.05681818181818181*xim1**2*xip1 +
        0.008522727272727272*xim1*xip1**2) +
     xi**2*(-0.03409090909090917*xim1**4 + 0.*xim1**5 -
        0.13636363636363635*xim1**3*xip1 -
        0.03409090909090909*xim1**2*xip1**2) +
     xi*(0.011363636363636243*xim1**5 + 0.*xim1**6 -
        0.045454545454545456*xim1**4*xip1 -
        0.03409090909090909*xim1**3*xip1**2))/
   (0.*xi**7 + 0.*xim1**7 + xi**6*(0.1875*xim1 + 0.*xip1) +
     0.*xim1**6*xip1 + 0.*xim1**5*xip1**2 + 0.*xim1**4*xip1**3 +
     xi**5*(-0.3125*xim1**2 - 0.3125*xim1*xip1 + 0.*xip1**2) +
     xi**3*xim1*(-0.25*xim1**3 + 1.25*xim1**2*xip1 + 1.25*xim1*xip1**2 -
        0.25*xip1**3) + xi**4*(-0.75*xim1**3 - 0.5*xim1**2*xip1 -
        0.75*xim1*xip1**2 + 0.*xip1**3) +
     xi**2*(0.*xim1**5 + 1.*xim1**4*xip1 + 3.*xim1**3*xip1**2 +
        1.*xim1**2*xip1**3) + xi*
      (0.*xim1**6 + 0.*xim1**5*xip1 + 1.*xim1**4*xip1**2 +
        1.*xim1**3*xip1**3)) + 1.*
   ((0.*xi**3 + xi**2*(0.2727272727272727 + 0.*xim1) -
        1.0454545454545454*xim1**2 + 0.*xim1**3 +
        xi*(-1.0909090909090908*xim1 + 0.*xim1**2 + 0.*xip1) -
        0.045454545454545456*xim1*xip1)/
      (-0.25*xi**3 + 1.*xi**2*xim1 + 1.*xi*xim1**2 + 0.*xim1**3) -
     (1.*(-0.25*xi**3 + 0.*xim1**3 + xi**2*(0.75*xim1 - 0.25*xip1) +
          1.*xim1**2*xip1 + 1.*xim1*xip1**2 +
          xi*(1.*xim1**2 + 2.*xim1*xip1 + 0.*xip1**2))*
        (0.*xi**7 + xi**6*(0. + 0.*xim1) + 0.*xim1**6 + 0.*xim1**7 +
          xi**5*(-0.2556818181818182*xim1 + 0.*xim1**2 + 0.*xip1) +
          0.*xim1**5*xip1 + 0.*xim1**4*xip1**2 +
          xi**4*(1.0227272727272727*xim1**2 + 0.*xim1**3 +
             0.03409090909090909*xim1*xip1 + 0.*xip1**2) +
          xi**3*(0.9517045454545452*xim1**3 + 0.*xim1**4 -
             0.05681818181818181*xim1**2*xip1 +
             0.008522727272727272*xim1*xip1**2) +
          xi**2*(-0.03409090909090917*xim1**4 + 0.*xim1**5 -
             0.13636363636363635*xim1**3*xip1 -
             0.03409090909090909*xim1**2*xip1**2) +
          xi*(0.011363636363636243*xim1**5 + 0.*xim1**6 -
             0.045454545454545456*xim1**4*xip1 -
             0.03409090909090909*xim1**3*xip1**2)))/
      ((-0.25*xi**3 + 1.*xi**2*xim1 + 1.*xi*xim1**2 + 0.*xim1**3)*
        (0.*xi**7 + 0.*xim1**7 + xi**6*(0.1875*xim1 + 0.*xip1) +
          0.*xim1**6*xip1 + 0.*xim1**5*xip1**2 + 0.*xim1**4*xip1**3 +
          xi**5*(-0.3125*xim1**2 - 0.3125*xim1*xip1 + 0.*xip1**2) +
          xi**3*xim1*(-0.25*xim1**3 + 1.25*xim1**2*xip1 +
             1.25*xim1*xip1**2 - 0.25*xip1**3) +
          xi**4*(-0.75*xim1**3 - 0.5*xim1**2*xip1 - 0.75*xim1*xip1**2 +
             0.*xip1**3) + xi**2*
           (0.*xim1**5 + 1.*xim1**4*xip1 + 3.*xim1**3*xip1**2 +
             1.*xim1**2*xip1**3) +
          xi*(0.*xim1**6 + 0.*xim1**5*xip1 + 1.*xim1**4*xip1**2 +
             1.*xim1**3*xip1**3)))) +
  1.*((1.0909090909090908 + 0.*xi + 0.*xim1)/(0.*xi - 1.*xim1) -
     (1.*(-1.*xi - 1.*xim1 - 1.*xip1)*
        (0.*xi**7 + xi**6*(0. + 0.*xim1) + 0.*xim1**6 + 0.*xim1**7 +
          xi**5*(-0.2556818181818182*xim1 + 0.*xim1**2 + 0.*xip1) +
          0.*xim1**5*xip1 + 0.*xim1**4*xip1**2 +
          xi**4*(1.0227272727272727*xim1**2 + 0.*xim1**3 +
             0.03409090909090909*xim1*xip1 + 0.*xip1**2) +
          xi**3*(0.9517045454545452*xim1**3 + 0.*xim1**4 -
             0.05681818181818181*xim1**2*xip1 +
             0.008522727272727272*xim1*xip1**2) +
          xi**2*(-0.03409090909090917*xim1**4 + 0.*xim1**5 -
             0.13636363636363635*xim1**3*xip1 -
             0.03409090909090909*xim1**2*xip1**2) +
          xi*(0.011363636363636243*xim1**5 + 0.*xim1**6 -
             0.045454545454545456*xim1**4*xip1 -
             0.03409090909090909*xim1**3*xip1**2)))/
      ((0.*xi - 1.*xim1)*(0.*xi**7 + 0.*xim1**7 +
          xi**6*(0.1875*xim1 + 0.*xip1) + 0.*xim1**6*xip1 +
          0.*xim1**5*xip1**2 + 0.*xim1**4*xip1**3 +
          xi**5*(-0.3125*xim1**2 - 0.3125*xim1*xip1 + 0.*xip1**2) +
          xi**3*xim1*(-0.25*xim1**3 + 1.25*xim1**2*xip1 +
             1.25*xim1*xip1**2 - 0.25*xip1**3) +
          xi**4*(-0.75*xim1**3 - 0.5*xim1**2*xip1 - 0.75*xim1*xip1**2 +
             0.*xip1**3) + xi**2*
           (0.*xim1**5 + 1.*xim1**4*xip1 + 3.*xim1**3*xip1**2 +
             1.*xim1**2*xip1**3) +
          xi*(0.*xim1**6 + 0.*xim1**5*xip1 + 1.*xim1**4*xip1**2 +
             1.*xim1**3*xip1**3))) -
     (1.*(-1.*xi - 1.*xim1)*((0.*xi**3 +
             xi**2*(0.2727272727272727 + 0.*xim1) -
             1.0454545454545454*xim1**2 + 0.*xim1**3 +
             xi*(-1.0909090909090908*xim1 + 0.*xim1**2 + 0.*xip1) -
             0.045454545454545456*xim1*xip1)/
           (-0.25*xi**3 + 1.*xi**2*xim1 + 1.*xi*xim1**2 + 0.*xim1**3) -
          (1.*(-0.25*xi**3 + 0.*xim1**3 + xi**2*(0.75*xim1 - 0.25*xip1) +
               1.*xim1**2*xip1 + 1.*xim1*xip1**2 +
               xi*(1.*xim1**2 + 2.*xim1*xip1 + 0.*xip1**2))*
             (0.*xi**7 + xi**6*(0. + 0.*xim1) + 0.*xim1**6 + 0.*xim1**7 +
               xi**5*(-0.2556818181818182*xim1 + 0.*xim1**2 + 0.*xip1) +
               0.*xim1**5*xip1 + 0.*xim1**4*xip1**2 +
               xi**4*(1.0227272727272727*xim1**2 + 0.*xim1**3 +
                  0.03409090909090909*xim1*xip1 + 0.*xip1**2) +
               xi**3*(0.9517045454545452*xim1**3 + 0.*xim1**4 -
                  0.05681818181818181*xim1**2*xip1 +
                  0.008522727272727272*xim1*xip1**2) +
               xi**2*(-0.03409090909090917*xim1**4 + 0.*xim1**5 -
                  0.13636363636363635*xim1**3*xip1 -
                  0.03409090909090909*xim1**2*xip1**2) +
               xi*(0.011363636363636243*xim1**5 + 0.*xim1**6 -
                  0.045454545454545456*xim1**4*xip1 -
                  0.03409090909090909*xim1**3*xip1**2)))/
           ((-0.25*xi**3 + 1.*xi**2*xim1 + 1.*xi*xim1**2 + 0.*xim1**3)*
             (0.*xi**7 + 0.*xim1**7 + xi**6*(0.1875*xim1 + 0.*xip1) +
               0.*xim1**6*xip1 + 0.*xim1**5*xip1**2 + 0.*xim1**4*xip1**3 +
               xi**5*(-0.3125*xim1**2 - 0.3125*xim1*xip1 + 0.*xip1**2) +
               xi**3*xim1*(-0.25*xim1**3 + 1.25*xim1**2*xip1 +
                  1.25*xim1*xip1**2 - 0.25*xip1**3) +
               xi**4*(-0.75*xim1**3 - 0.5*xim1**2*xip1 -
                  0.75*xim1*xip1**2 + 0.*xip1**3) +
               xi**2*(0.*xim1**5 + 1.*xim1**4*xip1 + 3.*xim1**3*xip1**2 +
                  1.*xim1**2*xip1**3) +
               xi*(0.*xim1**6 + 0.*xim1**5*xip1 + 1.*xim1**4*xip1**2 +
                  1.*xim1**3*xip1**3)))))/(0.*xi - 1.*xim1))
赞一下(1人) 回复此楼
7楼2012-06-14 12:23:52
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

chineserice

禁虫 (正式写手)
本帖内容被屏蔽
8楼2012-06-14 13:27:22
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

mshwangg

至尊木虫 (正式写手)
引用回帖:
8楼: Originally posted by chineserice at 2012-06-14 13:27:22
小木虫真是牛人辈出，昨天提出问题需找软件，发了两个帖子，两天解决问题~ o yeah~...

人外有人，每个人都有他的长处和短板。何必与别人计较短短长长？够用就好
赞一下 回复此楼
9楼2012-06-15 01:28:18
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lfc3583

铁杆木虫 (著名写手)
小木虫真是牛人辈出
赞一下 回复此楼
10楼2014-12-07 18:00:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chineserice 的主题更新
101/1返回列表
信息提示
关闭
请填处理意见
关闭