Example: The Kochanek-Bartels cubic spline interpolant

import java.text.*;
import com.imsl.math.*;

public class CsTCBEx1 {
    public static void main(String args[]) {

        double[]    xdata = {0., 1., 2., 3., 4., 5.};
        double[]    ydata = {5., 2., 3., 5., 1., 2.};
        double[]    continuity = {0., 0., 0., -1., 0., 0.};
        double[]    tension = {0., 0., 0., -1., 0., 0.};
        
        CsTCB cs = new CsTCB(xdata, ydata);
        cs.setContinuity(continuity);
        cs.setTension(tension);
        cs.compute();
        NumberFormat nf = NumberFormat.getInstance();
        nf.setMaximumFractionDigits(4);
        nf.setMinimumFractionDigits(4);
        System.out.println("      x    " +"   value   ");
        
        for (int k = 0;  k < 11;  k++) {
            double x = (double)k/2.0;
            double y = cs.value(x);
            System.out.println("   "+nf.format(x) + "     " +nf.format(y));
        }  
    }
}

Output

      x       value   
   0.0000     5.0000
   0.5000     3.4375
   1.0000     2.0000
   1.5000     2.1875
   2.0000     3.0000
   2.5000     3.6875
   3.0000     5.0000
   3.5000     2.1875
   4.0000     1.0000
   4.5000     1.2500
   5.0000     2.0000