package internewton;

//Autor : Rafał‚ Sikora
import java.util.Scanner;

public class InterNewton {

    public static void wczytajDane(double[] tablica)
    {
        Scanner scanner = new Scanner(System.in);
        for (int i = 0; i < 5; i++)
        {
            System.out.println("Podaj " + (i + 1) + " element : ");
            tablica[i] = scanner.nextInt();
        }
    }
    public static double Newton_InterpolationIR(double[] tablicax, double[] tablicay, double podstawowyX) {
        int index1, index2, temp;
        double pomocnicza, pomocnicza2, wynik = 0;
        for (index2 = 0; index2 <= tablicax.length - 1; index2++)
        {
            for (index1 = tablicay.length - 1; index1 > index2; index1--)
            {
                temp = index1 - index2 - 1;
                tablicay[index1] = (tablicay[index1] - tablicay[index1 - 1])
                        / (tablicax[index1] - tablicax[temp]);
            }
        }
        for (index1 = tablicax.length - 1; index1 >= 0; index1--)
        {
            pomocnicza = 1;
            for (index2 = 0; index2 <= index1 - 1; index2++)
            {
                pomocnicza2 = podstawowyX - tablicax[index2];
                pomocnicza = pomocnicza * pomocnicza2;
            }
            pomocnicza = pomocnicza * tablicay[index2];
            wynik = wynik + pomocnicza;
        }
        return wynik;
    }
    public static void main(String[] args) {
        double podstawowyX = 3;
        // Dane do programu : x= -4 -2 0 2 4   y= 537 37 1 45 553
        double[] mojeX = new double[5];
        double[] mojeY = new double[5];
        System.out.println("Elementy tablicy X: ");
        wczytajDane(mojeX);
        System.out.println("Elementy tablicy Y: ");
        wczytajDane(mojeY);
        System.out.println(Newton_InterpolationIR(mojeX, mojeY, podstawowyX));
    }
}