工程常用算法作业一-45组

《工程常用算法作业一-45组》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、题目：已知方程y=x3+2x2-ex，请分別用二分法、牛顿法和割线法求此方程的根。谋差限取:割线法xOpd*e£=10七。注意先确定出方程的有根区间。(初步确定有根区间为[0,1])二、程序流程图>X2=xTW(xl・)G/CKxl>ffM£)：l)《工程常用算法》综合实践作业一班级学号姓名主要「•作说明自评成绩编写程序排版，及总结绘制流程图作业名称：非线性方程求根完成日期：2015年05月30日开始'开始二分法丨

2、3-b«FXO=xlXI虫°愉lUx2fft三、完整的程序及简要的注释/7“%*4**1^%%*4**1^*4**4*1-^!•匸a

3、%*1**1^%%%*4**1^//“杠宀木容枉厅一半押杓杓和:押和押种//•上•土7厶•丄•*■•士•丄■・丄・乂■7■XVv「T"•丄■•土•丄■•土7•丄■•丄■•丄•*4*//不“*"水*不*****杆*那一分法水"和和不"*/#include#includevoidmain(){doublea=O,b=l,x,fx,fa;〃初始区间为[0,1]intk=0;printfC二分法：初始区间为：[0,l]“)；while(((b-a)/2)>pow(10,-12))〃终止条件{x=(a+b)/2;f

4、x=x*x*x+2*x*x-exp(x);fa=a*a*a+2*a*a-exp(a);if(fa*fx<0)b=x;〃变换区间elsea=x;k++;printf("k=%-3dx=%-.13ff(x)=%.13f",k,x,fx);〃输出次数，x值和函数值/「r

5、—)I//*x**7*%**x**7**4*%•*x**y^#includev

6、oidmain(){doublex0=0,x1,fx,fd,fx1,e=pow(10,-12);〃初始点为x0=0intk=0;printf(“牛顿法:初始点为x0=0H);while(1){fx=x0*x0*x0+2*x0*x0・exp(x0);〃计算x0处的函数值fd=3*x0*x0+4*x0-exp(x0);〃计算xO处的导数值printf("k=%-3dx=%・.13ff(x)=%-.13iAnH,k,x2,fx2);〃控制格式输出，x和f(x)小数点后显示13位if(fabs(fx2)

7、0=xl;〃交换处标xl=x2;getcharO;〃保持命令窗口始终打开四、程序运行结果E:PROGRAMFILEC语言'MICROSOFTVISUALSTUDIOMYPROJE

8、de#ineludevoidmain(){doublex0=1,x1=2,x2,fx0,fx1,fx2,fd,e二pow(10厂12);//初始点为xO=Lxl=2intk=();printf(“割线法:初始点:x0=l,xl=2n);while(1){fx0=x0*x0*x0+2*x0*x0-exp(x0);〃计算xO处前数值fxl=xl*xl*xl+2*xl*xl-exp(x1);〃计算xl处函数值fd=(fx1-fx0)/(x1-xO);〃用割线斜率代替切线斜率初始区间为:[0,13k=lX0.50

9、00000000000k=2X0■7500000000000k=3X0・8750000000000k=4X0.9375000000000k=5X0・9062500000000k=6X0.9218750000000k=7X-0.9296875000000k=8X0.9335937500000k=9X0・9316406250000k=10X=0.9306640625000k=llX0.9301757812500k=12X0.9299316406250k=13X0.9300537109375k=14X0.9301147460938k=15X0.9300

10、842285156k=16X0.9300994873047k=17X0.9300918579102k=18x=0.93009567260