循环结构程序设计之科学计算.docVIP

循环结构程序设计之科学计算 1 编写程序，用矩形法求一元函数f(x)=x^3+2x^2-12x-6（其中^表示幂运算），在区间[3.5，5]上的积分近似值s，保留2位小数，小区间数n=25，此参数不能改动，否则影响答案。 参考答案： 方法1，左85.31； 方法2，右90.69 答案： 1 /* #####矩形法(rectangle method) 编写程序，用矩形法求一元函数f(x)=x^3+2x^2-12x-6 （其中^表示幂运算），在区间[3.5，5]上的积分近似值S， 保留2位小数 （小区间数n=25，此参数不能改动，否则影响答案）。 参考答案： 方法1，左85.31 方法2，右90.69 */ #include stdio.h #include math.h main() { int n,k; double a,b,d,f0,s=0,s1,x; n=25; a=3.5; b=5; d=(b-a)/n; x=a; f0=pow(x,3)+2*pow(x,2)-12*x-6; for(k=1;k=n;k++) { //方法1，第1个矩形从函数曲线的左边开始计算 //只能使用两种方法之一,以下两行为第一种方法 s1=f0*d; s=s+s1; x=x+d; f0=pow(x,3)+2*pow(x,2)-12*x-6; //方法2,第1个矩形从函数曲线的右边开始计算 /* s1=f0*d; s=s+s1; */ } printf(s=%.2lf\n,s); } 2 编写程序，用矩形法求一元函数f(x)=ln(x+1)+x/2（其中ln为自然对数），在区间[1，5]上的积分近似值S，小区间数n=15，保留2位小数。 参考答案： 方法1，左10.95 方法2，右11.78 答案： 2 /* #####矩形法(rectangle method) 编写程序，用矩形法求一元函数f(x)=ln(x+1)+x/2 （其中ln为自然对数），在区间[1，5]上的积分近似值S， 保留2位小数（小区间数n=15，此参数不能改动，否则影响答案）。 参考答案： 方法1，左10.95 方法2，右11.78 */ #include stdio.h #include math.h main() { int n,k; double a,b,d,f0,s=0,s1,x; n=15; a=1; b=5; d=(b-a)/n; x=a; f0=log(x+1)+x/2; for(k=1;k=n;k++) { //方法1，第1个矩形从函数曲线的左边开始计算 //只能使用两种方法之一,以下两行为第一种方法 s1=f0*d; s=s+s1; x=x+d; f0=log(x+1)+x/2; //方法2,第1个矩形从函数曲线的右边开始计算 /* s1=f0*d; s=s+s1;*/ } printf(s=%.2lf\n,s); } 3 编写程序，用梯形法求一元函数f(x)=sin(2x)+3x，在区间[0，3.1416]上的积分近似值S，小区间数n=10,要求结果保留3位小数。) 参考答案: s=14.804 答案： 3 /* #####梯形法(keystoning)##### 编写程序，用梯形法求一元函数f(x)=sin(2x)+3x， 在区间[0，3.1416]上的积分近似值S，要求结果保留3位小数 （小区间数n=10）。 s=14.804 */ #includestdio.h #includemath.h main() { int i,n=10; double a=0,b=3.1416,d,t1,t2,s=0,x; d=(b-a)/n; for(i=1;i=n;i++) { x=a+(i-1)*d; t1=sin(2*x)+3*x;