home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #26 / NN_1992_26.iso / spool / sci / math / symbolic / 2918 < prev    next >
Encoding:
Internet Message Format  |  1992-11-07  |  3.9 KB

  1. Path: sparky!uunet!ogicse!emory!swrinde!cs.utexas.edu!hermes.chpc.utexas.edu!news.utdallas.edu!aicklen
  2. From: aicklen@utdallas.edu (Greg Aicklen)
  3. Newsgroups: sci.math.symbolic
  4. Subject: Inconsistent Mma results across platforms?
  5. Message-ID: <1992Nov6.154109.6056@utdallas.edu>
  6. Date: 6 Nov 92 15:41:09 GMT
  7. Article-I.D.: utdallas.1992Nov6.154109.6056
  8. Sender: aicklen@utdallas.edu (Gregory H. Aicklen)
  9. Organization: Univ. of Texas at Dallas
  10. Lines: 110
  11. Nntp-Posting-Host: drips.utdallas.edu
  12.  
  13. I am running the same version of Mma (2.1) on two different platforms and I
  14. get different results in a symbolic operation.  Note that I am running the
  15. Windows version of Mma 2.1, Student Edition on a 386SX with 10MB of RAM.
  16.  
  17. I seem to get different answers from ComplexExpand[] between the SPARC and
  18. PC (windows) implementations of Mma 2.1.  Here is an example:
  19.  
  20. Mathematica 2.1 for SPARC
  21. Copyright 1988-92 Wolfram Research, Inc.
  22.  -- X11 windows graphics initialized --
  23.  
  24. In[1]:= eta = .
  25.  
  26. In[2]:= xbar = .
  27.  
  28. In[3]:= q[x_,t_] := 2 eta Sech[2 eta (x -xbar)] E^(4 I eta^2 t)
  29.  
  30. In[4]:= eta = 1;
  31.  
  32. In[5]:= xbar = 0;
  33.  
  34. In[6]:= qt[x_,t_] = D[q[x,t], t]
  35.  
  36.              4 I t
  37. Out[6]= 8 I E      Sech[2 x]
  38.  
  39. In[7]:= qxx[x_,t_] = D[q[x,t], {x,2}]
  40.  
  41.             4 I t          3      4 I t                    2
  42. Out[7]= -8 E      Sech[2 x]  + 8 E      Sech[2 x] Tanh[2 x]
  43.  
  44. In[8]:= fnls[x_,t_] = qt[x,t] - I qxx[x,t] - 2 I q[x,t]^2 Conjugate[q[x,t]]
  45.  
  46.              4 I t                  8 I t              4 I t
  47. Out[8]= 8 I E      Sech[2 x] - 8 I E      Conjugate[2 E      Sech[2 x]]
  48.  
  49.                2          4 I t          3      4 I t                    2
  50. >     Sech[2 x]  - I (-8 E      Sech[2 x]  + 8 E      Sech[2 x] Tanh[2 x] )
  51.  
  52. In[9]:= nls = fnls[x,t]
  53.  
  54.              4 I t                  8 I t              4 I t
  55. Out[9]= 8 I E      Sech[2 x] - 8 I E      Conjugate[2 E      Sech[2 x]]
  56.  
  57.                2          4 I t          3      4 I t                    2
  58. >     Sech[2 x]  - I (-8 E      Sech[2 x]  + 8 E      Sech[2 x] Tanh[2 x] )
  59.  
  60. In[10]:= nls = ComplexExpand[nls, TargetFunctions->{Re, Im}]
  61.  
  62.                       3
  63.          -64 Cosh[2 x]  Sin[4 t]   16 Cosh[2 x] Sin[4 t]
  64. Out[10]= ----------------------- - --------------------- +
  65.                            3           1 + Cosh[4 x]
  66.             (1 + Cosh[4 x])
  67.  
  68.                  2  -32 Cos[8 t] Cosh[2 x] Sin[4 t]
  69. >    (4 Cosh[2 x]  (------------------------------- +
  70.                              1 + Cosh[4 x]
  71.  
  72.           32 Cos[4 t] Cosh[2 x] Sin[8 t]                    2
  73. >         ------------------------------)) / (1 + Cosh[4 x])  +
  74.                   1 + Cosh[4 x]
  75.  
  76.                                     2
  77.      16 Cosh[2 x] Sin[4 t] Sinh[4 x]
  78. >    -------------------------------- +
  79.                             3
  80.              (1 + Cosh[4 x])
  81.  
  82.                              3
  83.         64 Cos[4 t] Cosh[2 x]    16 Cos[4 t] Cosh[2 x]
  84. >    I (---------------------- + --------------------- +
  85.                           3          1 + Cosh[4 x]
  86.            (1 + Cosh[4 x])
  87.  
  88.                     2  -32 Cos[4 t] Cos[8 t] Cosh[2 x]
  89. >       (4 Cosh[2 x]  (------------------------------- -
  90.                                 1 + Cosh[4 x]
  91.  
  92.              32 Cosh[2 x] Sin[4 t] Sin[8 t]                    2
  93. >            ------------------------------)) / (1 + Cosh[4 x])  -
  94.                      1 + Cosh[4 x]
  95.  
  96.                                        2
  97.         16 Cos[4 t] Cosh[2 x] Sinh[4 x]
  98. >       --------------------------------)
  99.                                3
  100.                 (1 + Cosh[4 x])
  101.  
  102. In[11]:= Simplify[nls]
  103.  
  104. Out[11]= 0
  105.  
  106. In[12]:=Quit
  107.  
  108. For the Windows version 2.1 (Student Version), the result of the ComplexExpand[] is
  109.  
  110.          3
  111. Sech[2 x]  (8 (-2 + 2 Cos[8 t]) Sin[4 t] + 16 Cos[4 t] Sin[8 t]) +
  112.  
  113.             3
  114.  I Sech[2 x]  (-8 Cos[4 t] (-2 + 2 Cos[8 t]) + 16 Sin[4 t] Sin[8 t])
  115.  
  116. and I (obviously) can't reduce this to zero.
  117.  
  118.  
  119. Any ideas?  Anybody else able to duplicate this or seen something similar?
  120.  
  121. Greg Aicklen
  122. aicklen@utdallas.edu
  123.