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

实验报告学院（系）名称：计算机与通信工程

实验步骤：二分法：

迭代法：

牛顿法：

双点弦截法：

用二分法求方程x2-x-1=0的正根，要求准确到小数点后第一位#include

#include

#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)

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);

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

#include

#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)

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

#include

#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)//牛顿法方程

{

return exp(-x)-x;

}

double res3(double x0,double e,int n,double (*fun)(double x),double (*dfun)(double x)) //牛顿迭代法，奇异标志为返回INF，失败标志为返回INF，并输出Fail!

{

int k=1;

double x1;

if(fabs(dfun(x0)-0)

return INF;

x1=x0-fun(x0)/dfun(x0);

while(k!=n)

{

if(fabs(x1-x0)

return x1;

x0=x1;

x1=x0-fun(x0)/dfun(x0);

}

if(k==n)

printf("Fail!\n");

return INF;

}

int main()

{

printf("%.3f\n",res3(0.5,0.001,100,func2_2,dfunc2_2));

return 0;

}