anchor.diff:  
           let nnSpacing =
               unfoldedEgv
  -            |> nearesNeighbourSpacing unfoldedEgv.Length
  +            |> nearestNeighbourSpacing unfoldedEgv.Length
               |> Seq.toArray
           
           let histo =
  
positive.code:     let private evaluateWigner (xValue:float) =
          System.Math.PI*0.5*xValue*System.Math.E**(-System.Math.PI*0.25*xValue*xValue)    
  
      let private nearesNeighbourSpacing egvLength (uEgv:seq<float>) = 
          uEgv
          |> Seq.map (fun x -> x * float egvLength)
          |> Seq.sort
  
negative.code: 
      open FSharp.Stats
  
      let private spectralUnfolding egvalues =  
          let y = [0. .. 0.01 .. 1.] 
          let x = 
              egvalues
              |> Quantile.computePercentiles (Quantile.OfSorted.compute) y
          
          Interpolation.Approximation.approx x y egvalues (Seq.min)
          
  