熟悉用二分法,迭代法,牛顿法和弦截法求解非线性方程。.docVIP

实验报告 学院（系）名称：计算机与通信工程 姓名 赵云鹏 学号专业 计算机科学与技术 班级 09-1 实验项目 实验一 方程求根 课程名称 数值计算方法 课程代码 实验时间 2011年5月26日 实验地点 #7-215 批改意见：实验目的：熟悉用二分法，迭代法，牛顿法和弦截法求解非线性方程。 实验环境：硬件环境：IBM－PC或兼容机 软件环境：Windows操作系统 编程语言：C语言 实验内容： 1、用二分法求方程x2-x-1=0的正根，要求准确到小数点后第一位 2用迭代法和牛顿法求解方程x=e-x在x=0.5附近的一个根，要求精确到小数点后三位 3用双点弦截法求方程x3+3x2-x-9=0在区间[1,2]内的一个实根，精确到五位有效数字。 成绩 教师签字： 实验步骤：二分法： 迭代法： 牛顿法： 双点弦截法： 用二分法求方程x2-x-1=0的正根，要求准确到小数点后第一位 #include stdio.h #include math.h #define ADJUST1 0.01 #define ADJUST2 0.001 #define EX 0.000001 #define INF 999999999.99 double func1(double x)//二分法求的方程 { return (x*x-x-1); } double func2_1(double x)//迭代法的方程 { return exp(-x); } double res1(double a,double b,double (*fun1)(double x))//二分法 { double fa=fun1(a); double fb=fun1(b); double fmid=fun1((a+b)/2); while(fabs(b-a)=ADJUST1) { if(fabs(fmid-0a)EX) return fmid; else if(fmid*fa0) { b=(a+b)/2; fa=fun1(a); fb=fun1(b); fmid=fun1((a+b)/2); } else if(fmid*fb0) { a=(a+b)/2; fa=fun1(a); fb=fun1(b); fmid=fun1((a+b)/2); } } return (a+b)/2; } int main() { printf(%.2f\n,res1(0,100,func1)); printf(%.2f\n,func2_1(-1)); return 0; } 用迭代法和牛顿法求解方程x=e-x在x=0.5附近的一个根，要求精确到小数点后三位 #include stdio.h #include math.h #define ADJUST1 0.01 #define ADJUST2 0.001 #define EX 0.000001 #define INF 999999999.99 double func2_1(double x)//迭代法的方程 { return exp(-x); } double res2(double x0,double e,int n,double (*fun)(double x))//迭代法,迭代失败标志，输出Fail!，并返回INF { int k=1; double x1; x1=fun(x0); while(k!=n) { if(fabs(x1-x0)e) return x1; x0=x1; x1=fun(x0); k++; } if(k==n) printf(Fail!\n); return INF; } int main() { printf(%.3f\n,res2(0.5,0.001,100,func2_1));//q2.1 printf(%.3f\n,func2_1(0.567));//for test return 0; } 用双点弦截法求方程x3+3x2-x-9=0在区间[1,2]内的一个实根，精确到五位有效数字 #include stdio.h #include math.h #define ADJUST1 0.01 #define ADJUST2 0.001 #define EX 0.000001 #define INF 999999999.99 double dfunc2_2(double x)//牛顿法方程导数 { return (0-exp(-x)-1); } double func2_2(double x)//牛顿