Nullstellenberechnung durch Intervallhalbierung

Lyb schrieb:

Dann 2,5 * 2,5 -2 gibt 4,25

Wo in deinem Programm berechnest du das?

Lyb

fvonx1 = a * a -2;
fvonx2 = b * b -2;

Lyb schrieb:

fvonx1 = a * a -2;
fvonx2 = b * b -2;

a = 1; b = 4;
Wo berechnest du f(4,25) ?

Sorry.
a = 1; b = 4;
Wo berechnest du f(2.5) = 4.25?

Lyb

Erste Intervallhalbierung:

(1 + 4) / 2 = 2,5!

Diese 2,5 musst du nun wieder in die Gleich
x * x -2 einsetzen, um zu testen ob ein Vorzeichen-
wechsel statt fand.

2,5 * 2,5 -2 = 4,25

Lyb schrieb:

Diese 2,5 musst du nun wieder in die Gleich
x * x -2 einsetzen, um zu testen ob ein Vorzeichen-
wechsel statt fand.

2,5 * 2,5 -2 = 4,25

Das ist mir alles klar.

Nur in deinem Programm brechnest du das nicht!
Du musst doch prüfen ob der Vorzeichenwechsel im Bereich a,m oder m,b stattfindet.
Du berechnest aber nirgends den Funktionswert an der Stelle m.
Das fvonm = *m * *m -2; fehlt.

Lyb

Absolut richtig Dirk, ich danke dir!
Aber zusätzlich stimme auch noch die Prüfung im
if-statement nicht:

#include <stdio.h>
#include <stdlib.h>
/************************************************************************/
int bisection(double a, double b, double *m)
{
	int iterations = 0;
	double fvonx1 = 0;
	double fvonx2 = 0;
	double fvonm = 0;

	while (iterations < 100)
	{
		fvonx1 = a * a -2;
		fvonx2 = b * b -2;
		*m = (a + b) / 2;
		fvonm = *m * *m -2;

		if (fvonm < fvonx2 && fvonm > 0) {
			b = *m;
		}
		else {
			a = *m;
		}

		iterations++;

		if (iterations >= 100) {
			return 1;
		}
	}
	return 0;
}
/************************************************************************/
int main(void)
{
	double a = 1;
	double b = 4;
	double m = 0;

	if (bisection(a, b, &m) == 1) {
		printf("Die Berechnung wurde nach zu vielen Durchläufen beendet!\n");
	} else {
		printf("Die Nullstelle der Funktione lautet: %f\n", m);
	}
	system("PAUSE");
/* fertig! */
   exit(0);
}
/************************************************************************/

So funktioniert der Code
Also zumindest die Berechnung, das drum rum macht noch nicht
viel Sinn, weil er jedes mal abbrechen würde wegen zu vieler
Iterationen

Danke!

Da Abbruchkriterium aus deiner ersten Version war auch besser.

while (b - a > 1.0e-10)

Mach evtl. erstmal 1.0e-5.

Lyb

Bei dieser Bedingung geht er kein einziges mal in die Schleife rein.

Es geht: http://codepad.org/51EbPkks

Beim letzten mal waren a und b auch Zeiger.