Example: The cubic spline interpolant with periodic boundary conditions

A cubic spline interpolant to a function is computed. The value of the spline at point 0.23 is printed.
using System;
using Imsl.Math;

public class CsPeriodicEx1
{
	public static void  Main(String[] args)
	{
		int n = 11;
		double[] x = new double[n];
		double[] y = new double[n];
		
		double h = 2.0 * System.Math.PI / 15.0 / 10.0;
		for (int k = 0; k < n; k++)
		{
			x[k] = h * (double) (k);
			y[k] = System.Math.Sin(15.0 * x[k]);
		}
		
		CsPeriodic cs = new CsPeriodic(x, y);
		double csv = cs.Eval(0.23);
		Console.Out.WriteLine("The computed cubic spline value at " +
			                  "point .23 is " + csv);
	}
}

Output

The computed cubic spline value at point .23 is -0.303401472606451

Link to C# source.