Was stimmt nicht mit dieser Zirkulation?

Question

Wenn ich die while entferne und eine Bedingung in der if, else if, else auswähle, kann ich ein Ergebnis erhalten. Aber wenn ich die while Zirkulation in die Funktion auflöse, passiert eine Endlosschleife. Bitte finde heraus, wo das Problem liegt.

#include<stdio.h> 
#include<math.h> 
int a, b, c, d; 
float fun(float x); 
float solve(void); 

int main() 
{ 
    printf("Put in the coefficient of the equation:\n"); 
    scanf("%d", &a); 
    scanf("%d", &b); 
    scanf("%d", &c); 
    scanf("%d", &d); 

    float ans = solve(); 

    printf("A solve for the equation：x=%.3f", ans); 

    return 0; 
} 

float fun(float x) 
{ 
    float value = a * x * x * x + b * x * x + c * x + d; 

    return value; 
} 

float solve(void) 
{ 
    float x1 = -100, x2 = 100, x3; 
    float diff = fabs(fun(x1) - fun(x2)); 
    while (diff > 0.001) 
    { 
     x3 = (x1 * fun(x2) - x2 * fun(x1))/(fun(x2) - fun(x1)); 

     if (fun(x3) == 0) 
     { 
      x1 = x3; 
      x2 = x3; 
     } 
     else if ((fun(x3) * fun(x1)) < 0) 
     { 
      x2 = x3; 
     } 
     else 
     { 
      x1 = x3; 
     } 

     diff = fabs(fun(x1) - fun(x2)); 
    } 

    return diff; 
} 

Answer 1

Vor allem das ist nicht C++ !!! Wenn Sie diff zurückgeben, erhalten Sie immer einen Wert, der kleiner als 0,001 ist. Ich denke, du solltest x3 zurückgeben. Wenn Sie fragen fun(x3) == 0 Spaß (x3) muss vollständig gleich Null sein. Das erreichen Sie in den meisten Fällen niemals mit Floats. Besser geht es so:

 funx3=fun(x3); 

     if (funx3 == 0.0 || funx3>0.0 &&funx3<0.001 || funx3<0.0 && funx3>-0.001) 
     { 
      x1 = x3; 
      x2 = x3; 
     } 

Auch das wird nicht in allen Fällen funktionieren. Um das zu beheben, verdoppelt sich das in Ihrem gesamten Programm. Sie erhalten auch mehr Nachkommastellen. Ich bin nicht ganz sicher, ob alle Nachkommastellen stimmen. Vielleicht hast du einen Geralier geschnitten.

#include<stdio.h> 
#include<math.h> 
int a, b, c, d; 
double fun(float x); 
double solve(void); 

int main() 
{ 
    printf("Put in the coefficient of the equation:\n"); 
    scanf("%d", &a); 
    scanf("%d", &b); 
    scanf("%d", &c); 
    scanf("%d", &d); 

    float ans = solve(); 

    printf("A solve for the equation：x=%.8f", ans); 

    return 0; 
} 

double fun(float x) 
{ 
    double value = a * x * x * x + b * x * x + c * x + d; 

    return value; 
} 

double solve(void) 
{ 
    double x1 = -100, x2 =100, x3; 
    double diff = fabs(fun(x1) - fun(x2)); 
    double funx3; 
    while (diff > 0.0000001) 
    { 
     x3 = (x1 * fun(x2) - x2 * fun(x1))/(fun(x2) - fun(x1)); 

     funx3=fun(x3); 

     if (funx3 == 0.0 || funx3>0.0 &&funx3<0.0000001 || funx3<0.0 && funx3>-0.0000001) 
     { 
      x1 = x3; 
      x2 = x3; 
     } 
     else if ((fun(x3) * fun(x1)) < 0.0) 
     { 
      x2 = x3; 
     } 
     else 
     { 
      x1 = x3; 
     } 

     diff = fabs(fun(x1) - fun(x2)); 
    } 

    return x3; 
}