数值作业：牛顿迭代法

数值作业：牛顿迭代法

记录学习过程
写于2020.11.14

#include<stdio.h>  //调用标准函数库，使输入输出函数可用
#include<math.h>  //调用数学函数库
#define N 100
#define eps 1e-6  //定义全局变量eps（根的容许误差），delta（函数绝对值的容许误差）
#define eta 1e-8
double a, b, c, d;  //定义全局变量a,b,c,d

double f1(double x)  //定义double型子函数f1，传入x值进行运算（写在最前，则后续被调用时无需声明）
{
	return a * x * x * x + b * x * x + c * x + d;  //将运算和返回过程合并书写
}
double f2(double x)  //定义f1的导数为f2
{
	return (3 * x * x + 2 * b * x + c);  
}

void main()  //定义void型主函数，无返回参数
{
 
	double x0, x1, z;  //定义双精度型数组x[N]

	printf("请输入三元一次方程组的系数:a,b,c,d\n");
	scanf_s("%lf%lf%lf%lf", &a, &b, &c, &d);  //输入三元一次方程组的系数
loop:printf("请输入初值x(0)=\n");
	scanf_s("%lf", &x0);  //输入初值x0
	int t = 0;  //定义并赋值整型数t

	do    //do...whlie结构，先循环一次，再判断条件是否满足；条件为真，则继续循环
	{
		x1 = x0 - f1(x0) / f2(x0);  //牛顿迭代法公式
		if (t++ > N || fabs(f2(x1)) < eps)  //判断，当迭代次数t大于与迭代次数上限N 或 迭代值x1代入函数的导数的绝对值小于根的容许误差时执行其后语句
		{
			printf("\n牛顿迭代法发散\n");
			break;  //跳出do，while循环
		}
		z = fabs(x1) < 1 ? x1 - x0 : (x1 - x0) / x1;  //？为条件运算符，当满足fabs(x1)<1时，将x1-x0 的值赋给z;否则将(x1 - x0) / x1赋给z （追求迭代最快）
		x0 = x1;  //迭代，将x1的值赋给xo
		printf("x(%d)=%f\t", t, x0);  //输出第t次牛顿迭代法结果x(t)
		if (t % 3 == 0) printf("\n");  //隔三列换行
	} while (fabs(z) > eps && fabs(f1(x1)) > eta);  //&&逻辑与，全真则真；z的绝对值小于根的容许误差或函数值的绝对值小于其容许误差eta，终止循环
	printf("方程的近似根为%lf",x1 );

chong:printf("\n\n是否继续输入初值：\n");  //选择是否继续输入初值进行计算
	printf("1.继续输入初值\t 2.结束程序\n");
	printf("请选择：");
	int w;
	scanf_s("%d", &w);
	switch (w) //使用选择及结构语句switch，根据w的值选择方案
	{
	case 1:goto loop; break;  //方案1，跳转到loop标识符
	case 2:printf("程序结束！！！\n"); break;  //方案2结束整个程序
	default:    //当所输入w无法对应case时,运行其后语句
	{
		printf("输入错误，请重新输入！！！\n");
		goto chong;  //goto语句和标识符chong相搭配，程序无条件转移至chong处运行
	}
	}
}

结果图