紫面金刚王顺:c语言函数编程

来源:百度文库 编辑:查人人中国名人网 时间:2024/03/29 21:03:09
编写一 函数实现用牛顿迭代法求方程ax的3次方+bx的平方+cx+d=0在x=1附近的一个实根.主函数完成各系数值的输入及所求得的根值的输出.

#include <stdio.h>
#include <math.h>
float solut(float a,float b,float c, float d)
{
float x=1,x0,f,f1;
do
{
x0=x;
f=((a*x0+b)*x0+c)*x0+d;
f1=(3*a*x0+2*b)*x0+c;
x=x0-f/f1;
}
while(fabs(x-x0)>-1e-3);
return(x);
}
main()
{
float a,b,c,d;
printf("\ninput a,b,c,d:");
scanf("%f,%f,%f",&a,&b,&c,&d);
printf("equation is:%5.2fx^3+%5.2fx^2+%5.2fx+%5.2f=0",a,b,c,d);
printf(\nx=%10.7\n",solut(a,b,c,d);
}

看看以下这个程序吧
/*----------用牛顿迭代法求解方程的根-----------*/
#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#define LOWLIMIT 1.E-6
#define MAX 100
#define infinity 1.E300
double f(double x)
{
return x*x+3*x-3;/*依你的需要自定义修改函数*/
}
double diff(double x)
{
return 2*x+3; /*定义函数的导数*/
}
finite(double x)
{
if(x > infinity)return 0; /*此函数的功能是防止上溢*/
else return 1;
}
/*-------牛顿法迭代算法--------*/
void newton()
{
int k;
double x;
double delta_x;
printf("Enter x0,the initial guess:");
scanf("%lf",&x);
for(k=0;k<MAX;k++){
delta_x=-f(x)/diff(x);
if(!finite(delta_x))break;
x+=delta_x;
if(fabs(delta_x) < LOWLIMIT)break;
}
/*分析结果*/
if(k==MAX)
printf("\t No convergence in %d iterations.\n",MAX); /*与X轴无交点*/
else if(!finite(delta_x))
printf("\tOverflow on delta_x.\n");/*溢出*/
else
printf("\tProcess converged.\n"); /*与X轴有交点,成功求出根*/
printf("\tx=%e,delta_x=%e.\n",x,delta_x);
printf("f(x)=%e,f'(x)=%e.\n",f(x),diff(x));

}

void main()
{
int do_again;
puts("\n\n-------Solve f(x)=x^2+3x-3 by Newton's method.--------");
do{newton();
printf("\n\nEnter 1 to continue ,0 to exit:");
scanf("%d",&do_again);
}while(do_again);
}

dadalili不愧是专家,真热心。