Parameters
Double[] x
Sample points (N+1), sorted ascending
Double[] c0
Sample values (N or N+1) at the corresponding points; intercept, zero order coefficients
Double[] c1
Slopes (N) at the sample points (first order coefficients): N