MathGroup Archive 1992

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

Search the Archive


  • To: mathgroup at
  • Subject: AllRoots
  • From: akyildiz at ( Yilmaz Akyildiz tel 2-73558)
  • Date: Mon, 19 Oct 92 11:27:13 BST


The previous  version of the programme does not work for the initial
interval [0, 2 Pi].

It must be because of the limitations with RealInterval.

One way to get around this problem is to fatten the initial interval 
a bit in case it involves symbolic numbers.

The following version is OK now.

akyildiz at


lefthalf[RealInterval[{x_, y_}]]    := RealInterval[{x,(x+y)/2.}];             
righthalf[RealInterval[{x_, y_}]]   := RealInterval[{(x+y)/2.,y}];             
containszero[RealInterval[{x_,y_}]] :=  (x <= 0 <= y);                         
bisect[int_RealInterval]            := {lefthalf[int],righthalf[int]};         
intbisec[f_, a_, b_, tol_] :=                                         
      {step =   Select[(bisect/@#)//Flatten, containszero[N[f[#]]]&]& },          
      Nest[step, {RealInterval[N[{a-tol/2.,b+tol/2.}]]}, 

  • Prev by Date: Transparency in 3D objects?
  • Next by Date: Vvedensky's "PDE's with Mma" ?
  • Previous by thread: Re:AllRoots
  • Next by thread: quaternions