Question

我正在尝试编写一个算法来查找f（x）= x ^ 4 -4x +1的根我应该得到这个函数的4个根2实数和想象。我在c中写这个算法。但是不要，如果它写得很好，我应该为a和b输入什么样的初始猜测，因为每次我运行程序时它都会给我不同的数字这是我的代码，感谢您的帮助：

#include<stdio.h>
#include<math.h>
int computeroots(double, double, double, double);
int main()
{
    double a ,b, soln;
    double  epsilon = pow(10, -8);
    int MaxIter;
    printf("please enter an initial guess for a: \n");
    scanf("%ld",&a);
    printf("please enter a second guess for b: \n");
    scanf("%ld", &b);
    printf("please enter the maximum number of iteration: \n");
    scanf("%d", &MaxIter);
    soln = computeroots(a, b, epsilon, MaxIter);
    printf("Here is the solution: %d \n", soln);
    return 0;
}

int computeroots(double a, double b, double epsilon, double MaxIter)
{
    float FA = pow(a,4) - 4*a + 1;
    float FB = pow(b,4) - 4*b +1;
    float FP;
    int i =1;
    float p;
    if(FA * FB < 0)
    {
            while(i<MaxIter)
            {
                    p = a + (b-a)/2;
                    FP = pow(p,4) - 4*p +1;
                    if(FP == 0 || (b-a)/2 < epsilon)
                    {
                            return p;
                            break;
                    }
                    i++;
                    if(FA*FP > 0)
                    {
                            a = p;
                            FA = FP;
                    }
                    else
                    {
                            b = p;
                    }


            }
    }
}

Answer 1

我认为你必须找到两个猜测，这样一个猜测会让你产生负面结果而另一个会产生正面结果。维基百科有一个例子。迭代次数可以根据精度和初始猜测而大大改变。也不要将变量从double降为float，否则你将失去精度

c编程中的二分法

1 个答案: