// Całkowanie numeryczne - metoda Simpsona
// www.algorytm.org
// (c)2007 by Tomasz Lubinski

public class Metoda_Simpsona {

	//funkcja dla ktorej obliczamy calke
	private static double func(double x) {
		return x*x+3;
	}
	
	/**
	 * @param args
	 */
	public static void main(String[] args) {
		double xp, xk, dx, calka, s, x;
		int n;
		
		System.out.println("Podaj poczatek przedzialu calkowania");
		xp = Console.readDouble("");
		
		System.out.println("Podaj koniec przedzialu calkowania");
		xk = Console.readDouble("");

		System.out.println("Podaj dokladnosc calkowania");
		n = Console.readInt("");

		dx = (xk - xp) / (double)n;
		
		calka = 0;
        s = 0;
        for (int i=1; i<n; i++) {
                x = xp + i*dx;
                s += func(x - dx / 2);
                calka += func(x);
        }
        s += func(xk - dx / 2);
        calka = (dx/6) * (func(xp) + func(xk) + 2*calka + 4*s);
		
		System.out.println("Wartosc calki wynosi w przyblizeniu " + calka);

	}

}