Re: Why does this crash the kernel?
- To: mathgroup at yoda.physics.unc.edu
- Subject: Re: Why does this crash the kernel?
- From: David Withoff <withoff>
- Date: Fri, 31 Jul 1992 14:59:33 -0500
> The following file crashes a sun kernel and a mac kernel.
>
> b1 = Sqrt[2] x (500000 - 375 2^(5/2) x + 2 x^2 - 1000 p +
> 1.5 Sqrt[2] x p + 0.5 p^2)/(1000 - Sqrt[2] x - p);
>
> a1 = (500000 - 125 2^(7/2) x + x^2 - 2000 p + 2^(3/2) x p +
> 1.5 p^2)/(1000 - Sqrt[2] x - p);
>
> e = 0.125 b1^2 - 0.125 a1^2;
> kk = Simplify[e/.p->190];
>
> Solve[kk==0,x]
>
> Any workarounds?
>
> Howard Fink
If what you need are the zeros of kk, you can get them using FindRoot:
In[7]:= FindRoot[kk, {x, 0}]
Out[7]= {x -> 0.375406}
In[9]:= FindRoot[kk, {x, -1}]
Out[9]= {x -> -0.375349}
The crash problem in Solve (it's actually in NRoots) will be fixed in
the next release.
Dave Withoff
withoff at wri.com