/* Bisektionsverfahren 2.0 */

#include <stdio.h>
#include <stdlib.h>

double funktion (double x)            // hier wird das Polynom als Funktion definiert
{
        double erg = x*x*x*x*x+x+1;  // f(x) = x^5 + x + 1
        return erg;
}

int main (void)
{
        double a;                          // Beginn
        printf("\nGeben Sie einen Wert fuer a ein:\n");
        scanf("%lf",&a);
        double b;                          // Ende
        printf("\nGeben Sie einen Wert fuer b ein:\n");
        scanf("%lf",&b);
        double x = (a + b) / 2;         // Zwischenergebnisse und auch das Endergebnis
        double e;                       // die Genauigkeit Epsilon
        printf("\nGeben Sie einen Wert fuer e (z.B.:10e-6) ein:\n");
        scanf("%lf",&e);
        int summe = 0;
        int k = 1;

        while( (b - a) > e )
        {
                if(funktion(x) * funktion(b) < 0)
                {
                        a = x;
                }
                else
                {
                        b = x;
                }
                x = (a + b) / 2;
                summe = summe + k;
        }

        printf("a: %lf\nb: %lf\nMittelwert: %lf\n", a, b, x); // Ausgabe des Ergebnisses
        printf("\nAnzahl der durchlaufenen While-Schleifen:%i\n\n",summe);
        system("PAUSE");              // verhindert das Schliessen des Konsolen-Fensters
        return 0;
}