熟悉用二分法,迭代法,牛頓法和弦截法求解非線性方程

1、 實(shí)驗(yàn)報(bào)告學(xué)院（系）名稱：計(jì)算機(jī)與通信工程姓名 趙云鵬學(xué)號(hào)20091931專業(yè)計(jì)算機(jī)科學(xué)與技術(shù)班級(jí)09-1實(shí)驗(yàn)項(xiàng)目實(shí)驗(yàn)一 方程求根課程名稱數(shù)值計(jì)算方法課程代碼實(shí)驗(yàn)時(shí)間2011年5月26日實(shí)驗(yàn)地點(diǎn)#7-215批改意見：實(shí)驗(yàn)?zāi)康模菏煜び枚址?，迭代法，牛頓法和弦截法求解非線性方程。實(shí)驗(yàn)環(huán)境：硬件環(huán)境：IBMPC或兼容機(jī)軟件環(huán)境：Windows操作系統(tǒng)編程語言：C語言實(shí)驗(yàn)內(nèi)容：1、用二分法求方程x2-x-1=0的正根，要求準(zhǔn)確到小數(shù)點(diǎn)后第一位2用迭代法和牛頓法求解方程x=e-x在x=0.5附近的一個(gè)根，要求精確到小數(shù)點(diǎn)后三位3用雙點(diǎn)弦截法求方程x3+3x2-x-9=0在區(qū)間1,2內(nèi)的一個(gè)實(shí)根，精確到

2、五位有效數(shù)字。成績(jī)教師簽字：實(shí)驗(yàn)步驟：二分法：迭代法：牛頓法：雙點(diǎn)弦截法：n 用二分法求方程x2-x-1=0的正根，要求準(zhǔn)確到小數(shù)點(diǎn)后第一位#include <stdio.h>#include <math.h>#define ADJUST1 0.01#define ADJUST2 0.001#define EX 0.000001#define INF 999999999.99double func1(double x)/二分法求的方程return (x*x-x-1);double func2_1(double x)/迭代法的方程return exp(-x);double

3、 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*fa<0)b=(a+b)/2;fa=fun1(a);fb=fun1(b);fmid=fun1(a+b)/2);else if(fmid*fb<0)a=(a+b)/2;fa=fun1(a);fb=fun1(b);

4、fmid=fun1(a+b)/2);return (a+b)/2;int main()printf("%.2fn",res1(0,100,func1);printf("%.2fn",func2_1(-1);return 0;n 用迭代法和牛頓法求解方程x=e-x在x=0.5附近的一個(gè)根，要求精確到小數(shù)點(diǎn)后三位#include <stdio.h>#include <math.h>#define ADJUST1 0.01#define ADJUST2 0.001#define EX 0.000001#define INF 999999

5、999.99double func2_1(double x)/迭代法的方程return exp(-x);double res2(double x0,double e,int n,double (*fun)(double x)/迭代法,迭代失敗標(biāo)志，輸出Fail!，并返回INFint 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

6、n",res2(0.5,0.001,100,func2_1);/q2.1printf("%.3fn",func2_1(0.567);/for testreturn 0;用雙點(diǎn)弦截法求方程x3+3x2-x-9=0在區(qū)間1,2內(nèi)的一個(gè)實(shí)根，精確到五位有效數(shù)字#include <stdio.h>#include <math.h>#define ADJUST1 0.01#define ADJUST2 0.001#define EX 0.000001#define INF 999999999.99double dfunc2_2(double x)/牛

7、頓法方程導(dǎo)數(shù)return (0-exp(-x)-1);double func2_2(double x)/牛頓法方程return exp(-x)-x;double res3(double x0,double e,int n,double (*fun)(double x),double (*dfun)(double x) /牛頓迭代法 ，奇異標(biāo)志為返回INF，失敗標(biāo)志為返回INF，并輸出Fail!int k=1;double x1;if(fabs(dfun(x0)-0)<EX)return INF;x1=x0-fun(x0)/dfun(x0);while(k!=n)if(fabs(x1-x0)<e)return x1;x0=x1;x1=x0-fun(x0)/dfun(x0);if(k=n)printf