Что не так с этим тиражом?

Question

Если я удалю while и выберите одно условие в if, else if, else, я могу получить результат. Но когда я добавляю циркуляцию while в решение функции, происходит бесконечный цикл. Пожалуйста, выясните, где проблема?

#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

Прежде всего это c не C++ !!! Если вы возвращаете diff, вы каждый раз получаете значение меньше 0,001, как решение. Я думаю, вы должны вернуть x3. Если вы спросите fun(x3) == 0, то забава (x3) должна быть полностью равна нулю. Вы никогда не достигнете этого, в большинстве случаев с поплавками. Лучше сделайте так:

 funx3=fun(x3); 

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

Кроме того, это не будет работать во всех случаях. Чтобы исправить это использование, удваивается в вашей программе. Вы также получите больше десятичных знаков. Я не полностью уверен, что все десятичные разряды будут правы. Может быть, вы хотите нарезать одну eralier.

#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; 
}