A Mathematica based Version of the CKMfitter Package

68 Appendix D. Source Code
else
k = int( (dble(n-1)*ObsPred +xmax -dble(n)*xmin)/(xmax -xmin) )
SplinePred = ObsPred -tabsave(fn,0,k)
Tableau = tabsave(fn,2,k) +
> tabsave(fn,3,k) *SplinePred +
> tabsave(fn,4,k) *SplinePred**2 +
> tabsave(fn,5,k) *SplinePred**3
end if
End
D.2 dTableauO2
Tableau = Max(Tableau, 0D0)
!-----------------------------------------------------------------------
C
C Derivatives of Tableau function (divided by 2)
C
Double precision Function dTableauO2( ObsPred, fn )
Implicit None
include ’dimarray.f’
integer k, n, fn, lengthTable(maxLUT)
double precision ObsPred, SplinePred
double precision xmin, xmax
double precision tabsave(maxLUT,0:5,maxGranularity)
common /LUT/ lengthTable, tabsave
n = lengthTable(fn)
xmin = tabsave(fn,0,1)
xmax = tabsave(fn,0,n)
if( ObsPred.le.xmin .or. ObsPred.ge.xmax ) then