24小时热门版块排行榜    

查看: 2321  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

chineserice

禁虫 (正式写手)

本帖内容被屏蔽

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注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的回帖
查看全部 10 个回答

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的回帖

chineserice

禁虫 (正式写手)

本帖内容被屏蔽

6楼2012-06-14 12:00:10
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见