c语言迭代计算求方程的解,计算机辅助分析与设计(c语言二分法牛顿迭代法求方程解).doc...

计算机辅助分析与设计(c语言二分法牛顿迭代法求方程解)

《计算机辅助分析与设计》实验报告

实验一：

编程作业1：

用二分法求方程 的一个正的近似解(精确到0.1)。

程序：

#include

#include

#define F(x) ((x)*(x)*(x)-3*(x)+1) /*宏定义F(x)为方程式的左边*/

void main()

{int k=0;

float a,b,c,m;

printf("input a,b,c:");

scanf("%f,%f,%f",&a,&b,&c); /*输入a,b的值以及精确度c的值*/

while(F(a)*F(b)>0) /*输入的a,b的值直到满足F(a)F(b)0为止*/

{printf("input a,b again:");

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

for(;F(a)*F(b)<=0;k++) /*二分法开始求方程的根*/

{if(F(a)*F(b)==0) /*若F(a)F(b)=0,则方程的解在a,b两者中，并找出其根*/

if(F(a)==0)

printf("A root of equation is:%f k=%d\n",a,k);

else {printf("A root of equation is;%f k=%d\n",b,k);break;}

else m=(a+b)/2; /*取a,b的中点*/

if(fabs(a-b)

{ printf("A root of equation is:%f k=%d\n",m,k); break;}

else if(F(a)*F(m)>0) /*找出二分后区间两端点值异号的区间，并赋给[a,b]*/

a=m;

else b=m;

}

}

程序运行结果如下图

编程作业2：

用牛顿迭代法求方程 在附近的根。

#include

#include

float f(float x) /*定义函数f返回的值(其中f(x)为方程式的左边)*/

{

float y;

y=(x*x*x-x-1)/(3*x*x-1);

return(y);

}

void main()

{

float x1,x2,x0,c;

printf("input x0,c:");

scanf("%f,%f",&x0,&c); /*输入x0及精确度c的值*/

x2=x0;

do /*牛顿迭代法开始求方程的根*/

{x1=x2;

x2=x1-f(x1);}

while(fabs(f(x1))>c); /*若不满足精度要求执行循环，反之则得到方程的根就是x2*/

printf("A root of equation is:%f\n",x2);

}

程序运行结果如下图：

实验二：

编程作业1：

用梯形法计算定积分 ( 取 n=100)。

程序：

#include

#include

float f(float x) /*定义函数f返回[1/(1+)]*/

{

float y;

y=1/(x*x+1);

return(y);

}

void main()

{

float a,b,h,F=0;

int n,k;

printf("Input a,b,n:");

scanf("%f,%f,%d",&a,&b,&n); /*输入积分上下限a,b及n的值*/

F=(f(a)+f(b))/2.0;

h=fabs(b-a)/n; /*计算每一个小梯形的高度h*/

for(k=1;k

F+=f(a+k*h);

F=F*h;

printf("The integral of 1/(1+x*x) is:%f\n",F);

getch();

}

程序运行结果如下图：

编程作业2：

用抛物线法计算定积分 ( 取 n=100)。

程序:

#include

float f(float x) /*定义函数f返回[1/(1+)]*/

{

float y;

y=1/(1+x*x);

return(y);

}

void main()

{

float a,b,h,F;

int