数值计算方法编程作业（C语言版）
5星 · 超过95%的资源 需积分: 26 191 浏览量
2009-05-22
17:50:23
上传
评论 4
收藏 61KB DOC 举报 
1:第二章
(1)二分法求解非线性方程：
#include<stdio.h>
#include <math.h>
#define f(x) ((x*x-1)*x-1)
void main()
{ float a,b,x,eps;
int k=0;
printf("intput eps");/* 容许误差 */
scanf("%f",&eps);
printf("a,b=");
for(;;)
{scanf("%f, %f",&a ,&b);
if(f(a)*f(b)>=0) /* 判断是否符合二分法使用的条件 */
printf("二分法不可使用,请重新输入:");
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(" The root is x=%f, k=%d",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 <stdio.h>
#include <math.h>
float f(float x) /* 定义函数 f(x) */
剩余12页未读，继续阅读
资源评论
BelievenChanging2013-01-02很好，很好用。- GloriaD2011-12-20算法很全，适合做北航数值分析作业的同学们
