MathGroup Archive 1990

[Date Index] [Thread Index] [Author Index]

Search the Archive

Does anyone have an apparent fix for this mathematica problem?

  • To: mathgroup at yoda.ncsa.uiuc.edu
  • Subject: Does anyone have an apparent fix for this mathematica problem?
  • From: uunet!oregon.uoregon.edu!JOE
  • Date: Thu, 7 Jun 90 09:10 PDT

Has anyone else heard of the problem my users describe below? Is a patch
available to make Mathematica on the NeXT yield the same results as 
Mathematica on the Mac?

If you've any queries Fred needs to see or respond to, please CC: him directly
since he's not on the list right now.

Thanks,

Joe St Sauver (JOE at OREGON.UOREGON.EDU or JOE at OREGON)
Statistical Programmer and Consultant
University of Oregon Computing Center

>From:	OREGON::OLNESS       "Fredrick Olness (503) 346-5206 OLNESS at OREGON"  
>To:	OREGON::JOE
>
>Joe,
>
> .... One of Bob Zimmerman's students
>found what appers to be a bug in Mathematica on the Next.
>(Yes, I did  not believe it until I tried it myself.)
>
>I attach two examples of the problem below. 
>(For reference, the Mac gave the correct answer,
>so it is not a subtle error pertaining to branch cuts or 
>anything.)  I think he was going to try it on a different Next 
>machine for comparison.
>
>Thanks,
>
>Fred.
>
>--- example one --
>1% math
>Mathematica 1.2 (August 30, 1989) [With pre-loaded data]
>by S. Wolfram, D. Grayson, R. Maeder, H. Cejtin,
>   S. Omohundro, D. Ballman and J. Keiper
>with I. Rivin and D. Withoff
>Copyright 1988,1989 Wolfram Research Inc.
>
>In[1]:= InputForm[(b*c*(a + d - (a^2 + 4*b*c - 2*a*d + d^2)^(1/2))^3)/
>           (4*(a - d + (a^2 + 4*b*c - 2*a*d + d^2)^(1/2))*
>               (a^2 + 4*b*c - 2*a*d + d^2)^(1/2)) -
>        (b*c*(a +                 d + (a^2 + 4*b*c - 2*a*d + d^2)^(1/2))^3)/
>           (4*(a - d                 - (a^2 + 4*b*c - 2*a*d + d^2)^(1/2))*
>               (a^2 + 4*b        *c - 2*a*d + d^2)^(1/2))]
>
>Out[1]//InputForm=
>  (b*c*(a + d - (a^2 + 4*b*c - 2*a*d + d^2)^(1/2))^3)/
>    (4*(a^2 + 4*b*c - 2*a*d + d^2)^(1/2)*
>      (a - d + (a^2 + 4*b*c - 2*a*d + d^2)^(1/2))) -
>   (b*c*(a + d + (a^2 + 4*b*c - 2*a*d + d^2)^(1/2))^3)/
>    (4*(a^2 + 4*b*c - 2*a*d + d^2)^(1/2)*
>      (a - d - (a^2 + 4*b*c - 2*a*d + d^2)^(1/2)))
>
>In[2]:= rule={a->1,b->2,c->3,d->4}
>
>Out[2]= {a -> 1, b -> 2, c -> 3, d -> 4}
>
>In[3]:= N[res//.rule]
>
>Out[3]= res
>
>In[4]:= res=%1;
>
>In[5]:= N[res//.rule]
>
>Out[5]= 37.
>
>In[6]:= N[Together[Expand[res]]//.rule]
>
>Out[6]= 47.2208
>
>In[7]:= N[Simplify[res]//.rule]
>
>Out[7]= 47.2208
>
>-- example two follows --
>
>1 % math
>Mathematica 1.2 (August 30, 1989) [With pre-loaded data]
>by S. Wolfram, D. Grayson, R. Maeder, H. Cejtin,
>   S. Omohundro, D. Ballman and J. Keiper
>with I. Rivin and D. Withoff
>Copyright 1988,1989 Wolfram Research Inc.
>
>In[1]:= InputForm[(b*c*(a + d - (a^2 - 2*a*c + 4*b*c + d^2)^(1/2))^3)/
>           (4*(a - d + (a^2 - 2*a*c + 4*b*c + d^2)^(1/2))*
>               (a^2 - 2*a*c + 4*b*c + d^2)^(1/2)) -
>        (b*c*(a + d + (a^2 - 2*a*c + 4*b*c + d^2)^(1/2))^3)/
>           (4*(a - d - (a^2 - 2*a*c + 4*b*c + d^2)^(1/2))*
>               (a^2 - 2*a        *c + 4*b*c + d^2)^(1/2))]
>
>
>Out[1]//InputForm=
>  (b*c*(a + d - (a^2 - 2*a*c + 4*b*c + d^2)^(1/2))^3)/
>    (4*(a^2 - 2*a*c + 4*b*c + d^2)^(1/2)*
>      (a - d + (a^2 - 2*a*c + 4*b*c + d^2)^(1/2))) -
>   (b*c*(a + d + (a^2 - 2*a*c + 4*b*c + d^2)^(1/2))^3)/
>    (4*(a^2 - 2*a*c + 4*b*c + d^2)^(1/2)*
>      (a - d - (a^2 - 2*a*c + 4*b*c + d^2)^(1/2)))
>
>In[2]:= res=%;
>
>In[3]:= rule={a->1,b->2,c->3,d->4}
>
>Out[3]= {a -> 1, b -> 2, c -> 3, d -> 4}
>
>In[4]:= N[res//.rule]
>
>Out[4]= 36.9231
>
>In[5]:= N[Together[Expand[res]]//.rule]
>
>Out[5]= 46.9906
>
>In[6]:= N[Simplify[res]//.rule]
>
>Out[6]= 46.9906


  • Prev by Date: Re: Gaussians
  • Next by Date: Re:Does anyone have an apparent fix for this mathematica problem?
  • Previous by thread: mathematical relations
  • Next by thread: Re:Does anyone have an apparent fix for this mathematica problem?