solve函数c语言_数值计算方法编程作业(C语言版)

1:

第二章

(1)

二分法求解非线性方程：

#include

#include 

#define f(x) ((x*x-1)*x-1)

void main()

{

float a,b,x,eps;

int k=0;

printf("intput eps\n");/*

容许误差

*/

scanf("%f",&eps);

printf("a,b=\n");

for(;;)

{scanf("%f, %f",&a ,&b);

if(f(a)*f(b)>=0) /*

判断是否符合二分法使用的条件

*/

printf("

二分法不可使用

,

请重新输入

:\n");

else break;

}

do

{

x=(a+b)/2;

k++;

if(f(a)*f(x)<0) /*

如果

f(a)*f(x)<0

，则根在区间的左半部分

*/

b=x;

else if(f(a)*f(x)>0)

/*

否则根在区间的右半部分

*/

a=x;

else break;

}while(fabs(b-a)>eps);/*

判断是否达到精度要求

,

若没有达到

,

继续循环

*/

x=(a+b)/2; /*

取最后的小区间中点作为根的近似值

*/

printf("\n The root is x=%f, k=%d\n",x,k);

}

运行结果：

intput eps

0.00001

a,b=

2,-5

The root is x=1.324721, k=20

Press any key to continue

总结：本题关键在于两个端点的取值和误差的判断，此程序较容易。二分法收敛速度较快，

但缺点是只能求解单根。

(2)

牛顿法求解非线性方程：

#include 

#include 

float f(float x)

/*

定义函数

f(x) */