以y=a^3*x+b^2*x+c*x+d+sin(e)为例:
//用牛顿迭代法求方程的根
#include<stdio.h>
#include<math.h>
double a, b, c, d, e;
double f1(double x)//返回y的原函数
{
return a * pow(x, 3) + b * pow(x, 2) + c * x + d + sin(e);
}
double f2(double x)//返回y函数的一阶导函数
{
return 3 * a * pow(x, 2) + 2 * b * x + c;
}
void main()//主程序
{
scanf("%lf %lf %lf %lf %lf", &a, &b, &c, &d, &e);//从主程序输入系数
double x1 = 1, x;//求1附近的根
do
{
x = x1;
x1 = x - f1(x) / f2(x);//迭代关系式
} while (fabs(x1 - x) > 1e-8);//对迭代过程进行控制
printf("%.3f", x1);//从主程序输出实根
}
该方法一大特点是在外围定义了系数a,b,c,d,e,避免了在各个函数内重复定义的繁琐过程。