Heron-Verfahren mit "C"

08.05.2013, 23:06

Marshal

Heron-Verfahren mit "C"

Meine Frage:
Hi Leute.

Ich soll ein Programm schreiben, das mittels des Heronverfahrens $[latex] x_{n+1}=\frac{1}{2}\left(x_n+\frac{a}{x_n}\right) [/latex]$ die Wurzel einer nicht-negativen vom Nutzer eingegebenen double Zahl [latex] a [/latex]

berechnet.

Im Quellcode habe ich den Startwert [latex] x_0=a [/latex]

gesetzt.

code:

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:

#include <stdio.h>
#include <math.h>

int main(void)
{
    double a, x_alt, x_neu;

    /* Radikant inizialisieren */

    printf("Radikand a eingeben:\n\n");
    printf("a = ");
    scanf("%lf", &a);

    /* falls Radikand negativ */

    if( a < 0. )
    {
        printf("\n\nFehler: Radikand ist negativ\n");
        return 0;
    }

    /* falls Radikand gueltig */

    for(x_alt = a; fabs( x_neu - sqrt(a) ) >= 1e-2; x_alt == x_neu)
        x_neu = 0.5*(x_alt + a/x_alt);
   
    printf("\nWurzel %f = %f\n", a, x_neu);
    return 0;
}

Meine Ideen:
Allerdings kommt keine Ausgabe, wenn ich eine nicht-negative double-Zahl eingebe... Ich vermute eine Endlosschleife...

Sieht jemand den oder die Fehler?

09.05.2013, 00:20

Auf diesen Beitrag antworten »

Airblader

code:
1:	for(x_alt = a; fabs( x_neu - sqrt(a) ) >= 1e-2; x_alt == x_neu)

Hier verbirgt sich der Fehler: Du machst hinten einen Vergleich, willst aber eigentlich eine Zuweisung machen.

Im Übrigen ist die Abbruchbedingung gemogelt: Wenn du $[latex]\sqrt{a}[/latex]$ berechnen sollst, dann darfst du diesen Wert natürlich nicht in der Abbruchbedingung stehen haben. Du musst dir hier was anderes einfallen lassen.

Neue Frage »

Antworten »

Heron-Verfahren mit "C"

Verwandte Themen