anchor.diff:              let gcvfinal = gcvmin.[eminIndex]
               let lambdafinal = lambdafin.[eminIndex]
               let cfinal = cfin.[eminIndex]
  -            createHermiteResultSimple afinal cfinal gcvfinal efinal lambdafinal
  -
  -
  +            createHermiteResultSimple afinal cfinal gcvfinal efinal lambdafinal ctemp nan
   
       let splineIncreasing (x:Vector<float>) (y:Vector<float>) (W:Matrix<float>) lambd con= 
           spline (~-) x y W lambd con
  
positive.code:                 let efin = e.[gcvminIndex]
                  let afin = a.[gcvminIndex]
                  let cfin = c.[gcvminIndex]
                  let lambdafin = 0.01*(1.2**(float gcvminIndex))
                  ((efin,gcvmin),afin,lambdafin,cfin)
                      )
          |> unzip4List
  
negative.code:                     if maxs.[0] < mins.[0] then 3
                      else -3
  
          createHermiteResult fit.TraceA fit.TraceC fit.Error fit.GCV fit.Lambda mins maxs clNew fit.Ctemp fit.Variance
  
          
  
  