问题为x^0.5=X;x为非负实数。
一、运用牛顿迭代法
两边平方移项得 X^2-x=0,可以用牛顿迭代法求解二元一次方程的根
#include<stdio.h>
#include<math.h>
double func(double x,double c) //要求解的函数方程
{
return x*x-c;
}
double func1(double x) //导函数
{
return 2*x;
}
int Newton(double *x, double precision, int maxcyc,double c)//maxcyc 最大迭代次数,c被开方数
{
double x1, x0;
int k;
x0 = *x;
for (k = 0; k < maxcyc; k++)
{
if (func1(x0) == 0.0)
{
printf("迭代过程中导数为0\n");
return 0;
}
x1 = x0 - func(x0,c) / func1(x0);
if (fabs(x1 - x0) < precision || fabs(func(x1,c)) < precision)//达到设定的精度
{
*x = x1;
return 1;
}
else
x0 = x1; //准备下次迭代
}
printf("迭代次数超过预期!\n"); //迭代次数达到,仍没有达到精度
return 0;
}
int main()
{
double x, precision,c;
int maxcyc;
printf("输入被开方数");
scanf("%lf",&c);
printf("输入初始迭代值x0");
scanf("%lf", &x);
printf("输入最大迭代次数:");
scanf("%d", &maxcyc);
printf("请输入迭代要求的精度:");
scanf("%lf", &precision);
if (Newton(&x, precision, maxcyc,c) == 1)
printf("该值附近的根为:%lf\n", x);
else
printf("迭代失败!\n");
while(1);
return 0;
}
本文介绍了如何使用牛顿迭代法来实现开根运算,详细解析了针对方程X^2 - x = 0的迭代求解过程。
1万+

被折叠的 条评论
为什么被折叠?



