/*
 * 示例程序Sample6_13: Interpolation類的第一種邊界條件的三次樣條函數插值、微商與積分
 */
package javaalgorithm.sample;

import javaalgorithm.algorithm.Interpolation;

public class Sample6_13 
{
	public static void main(String[] args)
	{
		// 數據
		int n = 12;
		double[] x = {0.52,8.0,17.95,28.65,50.65,104.6,156.6,260.7,364.4,468.0,507.0,520.0};
		double[] y = {5.28794,13.84,20.2,24.9,31.1,36.5,36.6,31.0,20.9,7.8,1.5,0.2};
		double[] t = {4.0,14.0,30.0,60.0,130.0,230.0,450.0,515.0};
		
		// 微分
		double[] dy = new double[n];
		double[] ddy = new double[n];
		dy[0] = 1.86548;
		dy[n-1] = -0.046115;

		// 插值運算
		int m = 8;
		double[] z = new double[m];
		double[] dz = new double[m];
		double[] ddz = new double[m];
		
		double yt = Interpolation.getValueSpline1(n, x, y, dy, ddy, m, t, z, dz, ddz);
		System.out.println("定積分值 = " + yt);
		System.out.println("\nt[i]\t\tz[i]\t\t\tdz[i]\t\t\tddz[i]\n");
		for (int i=0; i<m; ++i)
		{
			System.out.println(t[i] + "\t" + z[i] + "\t" + dz[i] + "\t" + ddz[i]);
		}
	}
}