牛顿迭代法求解方程

最新推荐文章于 2024-01-13 14:55:07 发布

TA-MEI

最新推荐文章于 2024-01-13 14:55:07 发布

阅读量493

点赞数

分类专栏： C语言文章标签：算法 c语言

本文链接：https://blog.csdn.net/weixin_43826113/article/details/107383070

版权

C语言专栏收录该内容

2 篇文章 0 订阅

订阅专栏

C语言编程百例

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

int Function(double,double *,double *);
int Newton(double *,double,int);

int Function(x,f,dy)
double x;
double *f;
double *dy;
{
	*f = x*x*(x-1)-1;
	*dy = 3*x*x-2*x;
	return(1);
}

int Newton(x,eps,l)
double *x;
double eps;
int l;
{
	double f,dy,x1;
	Function(*x,&f,&dy);
A:	if(fabs(dy) == 0)
	{
		l = 0;
		return (0);
	}
	x1=*x-f/dy;
	Function(x1,&f,&dy);
	if(fabs(x1-*x)>=eps||fabs(f)>=eps)
	{
		l-=1;
		*x=x1;
		if(l==0)
			return(1);
		goto A;
	}
	*x = x1;
	return 1;
}

void main()
{
	double x,eps;
	int l;
	eps=1.e-6;
	x=1.5;
	l=60;
	if(!Newton(&x,eps,l))
	{
		printf("该函数不可以用牛顿跌代法求根!\n");
	}
	printf("利用牛顿跌代法求的的根为：\n");
	printf("x=%.10f\n",x);
}