学工程计算作业.docVIP

学工程计算作业 姓名：刘欣 班级：统计1201 学号：201247100125 学院：数计学院 1,解方程x^3-3x-1=0在【1,2】的近似根(精确到10^-3)(二分法）.#include stdio.h #include math.h main() { double x1,x2,x,e,y; x1=1.0;x2=2.0; do { e=x1-x2; x=(x1+x2)/2; y=pow(x,3.0)-3*x-1.0; if(y!=0) { if(y0) x2=x; else x1=x; } else printf(%0.20f\n,x); } while(fabs(e)1.0e-10); printf(近似根是：%0.20f\n,x); return x; } 说明：运行程序后，显示：“近似根是：1.87938524156925270000”，结果位数与e的截断精度有关。 2，解方程x*e^x-1=0.(牛顿迭代） #includestdio.h #includemath.h int main() { double f(double x); double x,y,d; x=1; y=0; d=0.000000001; while(fabs(x-y)=d) { x=y; y=f(x); } printf(%.3f\n,y); printf(%f\n,exp(y)*y); return 0; } double f(double x) { return x-(x*exp(x)-1)/(exp(x)+x*exp(x)); } 说明：结果显示：“0书中表2.9给出了迭代四次后的结果为0结果表明算法收敛，有四位有效数字。 3，解方程x*e^x-1=0.(简单迭代） #includestdio.h #includemath.h main() { ; double x0,x1,d; x0; x1=0.0;//迭代值。 d=0.000000001;//误差控制 while(fabs(x0 -x1)=d) { x0=x1; x1=exp(-x0); } printf(%.3f\n,x1); printf(%f\n,exp(x1)*x1); return x1; } 说明：结果显示与3相同。 4,二分法搜索x^3-1.8x^2+0.15x+0.65=0的有根区间。 #includestdio.h #includemath.h main() { double a,x,y,k; a=-2; k=0.5; do { x=a+k; y=x*x*x-1.8*x*x+0.15*x+0.65; if(y!=0) { if(y0) printf(区间[%f,%f]\n,a,x); else a=x;t=t+1;printF(循环次数是%d\n,t); } else printf(精确值是%f\n,x); }while(y0); } 说明：结果显示；【0.5,1】，与用逐步搜索法笔算结果一致，循环次数为3次。 matlab语言实现。 f=inline(x^3-1.8x^2+0.15x+0.65) %把f=x^3-1.8x^2+0.15x+0.65设定为内嵌函数 x1=1.0,x2=2.0 %给定区间[a,b]即[x1,x2]一个初值，这里为[1.0,2.0]. err=1e-3; %给出精度为1E-3。 k=0; %k为一个计数变量，为了统计符合精度时，二分的循环次数。给初值0. while abs(x1-x2)err; %二分循环的结束条件，最后一次的解变量与倒数第二次的解变量 %的差的绝对值，比精度err小时，则循环结束，否则，循环继续。 x=(x1+x2)/2; %给方程的解变量赋一个初值，为区间的中点。 k=k+1; %循环一次k自加1. format long %对从此行一下的变量进行长精度格式化。 x %这里的x的作用是，显示每次循环时二分之后的当前x的值。 if(f(x)0) %若为负数即f(x)0,则把（（x1+x2）/2）赋给x1,即新区间的左端点。 x1=x; else if(f(x)0)