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
>
>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
>
>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?