您的位置:首页 > 其它

3元一次方程(牛顿迭代法求方程的根)

2014-07-15 17:03 246 查看
牛顿迭代方求方程根的公式原理请自行谷歌或百度相关资料,具体代码如下:

// test1.cpp : 定义控制台应用程序的入口点。
//

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

#define t_3
//#define t_2

int _tmain(int argc, _TCHAR* argv[])
{
float solution(float a, float b, float c, float d);
float solution1(float a,float b, float c);
float a,b,c,d;
#ifdef t_3
printf("input a,b,c,d:");
scanf("%f,%f,%f,%f", &a,&b,&c,&d);
printf("The solution is: x=%10.7f\n", solution(a,b,c,d));
#else if t_2
printf("input a,b,c:");
scanf("%f,%f,%f", &a,&b,&c);
printf("The solution is: x=%10.7f\n", solution1(a,b,c));
#endif
return 0;
}

float solution(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-5);
return x;
}

float solution1(float a,float b, float c)
{
float x=1,x0,f,f1;
do
{
x0=x;
f=(a*x0+b)*x0+c;
f1=(3*a)*x0+b;
x=x0-f/f1;
}while(fabs(x-x0)>=1e-5);
return x;
}


1.求解精度可根据自身情况进行调节,目前精度为1e-5

2.该代码在进行2元一次方程求解时,未考虑负值解。具体2元一次方程求解直接参考http://baike.baidu.com/view/397767.htm?fr=aladdin 求解即可。


内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: