Question

c语言 用迭代法求解

键盘输入一个实数a，输出它的立方根。请自行设计算法编写程序，不允许调用pow函数。程序用3_18.c命名并上交。
【提示】
假设a的立方根为x，则x3=a，问题转化求一元高次方程x3-a=0的解

Answer 1

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

int main()
{
double a;
double x2=0,x1=1.0;
printf("请输入被开方数：");
scanf("%lf",&a);
for(x1=1.0;fabs(x2-x1)>1e-10;)
{
x1=x2;
x2=(2*x1+a/(x1*x1*3.0));
}
printf("结果为：%f",x2);
return 0;
}

Answer 2

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

int main()
{
double a;
double temp,xk2,xk1=0;
printf("请输入一实数：");
scanf("%lf",&a);
do{
xk2=(2*xk1+a/(xk1*xk1))/3.0; //Newton 法
temp=xk1;
xk1=xk2;
} while(fabs(xk2-temp)<1e-6)
return 0;
}

追问

怎么运行结果不对呢

追答

刚刚少了个分号

#include
#include

int main()
{
double a;
double temp,xk2,xk1=1.0;
printf("请输入一实数：");
scanf("%lf",&a);
do{
xk2=(2*xk1+a/(xk1*xk1))/3.0; //Newton 法
temp=xk1;
xk1=xk2;
}while(fabs(xk2-temp)<1e-6);
printf("所得结果：%f",xk2);
return 0;
}

Answer 3

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

double fun(double a, double x)
{
return x*x*x-a;
}

int main()
{
double a,i,j;
scanf("%lf",&a);
if(a==0) printf("%f\n",0.0);
else if(a>0)
{
i=0;j=1;
while(fun(a,j)<0)
{i=j;j*=2;}
}
else if(a<0)
{
i=-1;j=0;
while(fun(a,i)>0)
{j=i;i*=2;}
}
if(a!=0)
{
while(fabs(fun(a,(i+j)/2))>0.00001) // 二分法
{
if(fun(a,(i+j)/2)>0) j = (i+j)/2;
else if(fun(a,(i+j)/2)<0) i = (i+j)/2;
else
{
i = j = (i+j)/2;
break;
}
}
i = (i+j)/2;
printf("%f\n",i);
}
}