控制结构与语句.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 例 用迭代公式求a的平方根，迭代公式如下： #includestdio.h #include math.h #includestdlib.h void main() { double x0,x1,a; printf(输入数据a:\n); scanf(%lf,a); if (a0) { printf(data error\n); exit(0); } x0=a/2; //给出一个任意值作为迭代的初始值x0 do { x1=(x0+a/x0)/2; x0=x1; //本次计算结果x1作为下次迭代式中的x0值 } while(fabs(x1*x1-a)1e-5); printf(a的平方根=%.5lf\n,x1); } x y f(xn+1) f(xn) 例 编程，求∫12(x2+x+1)dx ⑴用矩形公式求函数f(x)在[a,b]区间的定积分。 所求定积分的近似值： #include stdio.h void main() { float a,b,h,x,s=0;int n,i; scanf(“%f%f%d”,a,b,n); h=(b-a)/n; for(i=0;in;i++) { x=a+i*h; s+=x*x+x+1; } printf(“%.10f\n”,s*h); } f(x) x0 x1 x2 x3 x4 x5 x6 x7 x8 h 例 编程，求∫12(x2+x+1)dx ⑵用梯形公式求函数f(x)在[a,b]区间的定积分。 所求定积分的近似值： n-1 ∫abf(x)dx=?[f(a+i?h)+f(a+(i+1)?h)]?h/2 i=0 #include stdio.h void main() { float a,b,h,x1,x2,s=0;int n,i; scanf(“%f%f%d”,a,b,n); h=(b-a)/n; for(i=0;in;i++) { x1=a+i*h; x2=x1+h; s+=(x1*x1+x1+1+x2*x2+x2+1)*h/2; } printf(“%.10f\n”,s); } f(x) x0 x1 x2 x3 x4 x5 x6 x7 x8 h 例 编程，求∫12(x2+x+1)dx ⑶求数值积分的自适应方法 ①先m等份[a,b]区间，求的积分近似值s1 ②使m=2*m后，再m等份[a,b]区间，求的积分近似值s2 ③如果|s1-s2|≥ε，则执行①，否则求的积分值s2 #include stdio.h #include math.h void main() { float a,b,h,x,s1,s2,eps;int m,i; scanf(“%f%f%f”,a,b,eps);m=4; do { s1=0;h=(b-a)/m;x=a; for(i=0;im;i++) {s1+=(x*x+x+1+(x+h)*(x+h)+(x+h)+1)*h/2;x+=h;} s2=0;m*=2;h=(b-a)/m;x=a; for(i=0;im;i++) {s2+=(x*x+x+1+(x+h)*(x+h)+(x+h)+1)*h/2;x+=h;} }while(fabs(s1-s2)=eps); printf(“%.10f\n”,s1); } f(x) x0 x1 x2 x3 x4 x5 x6 x7 x8 h * * * * * * * * * * * * * * * * * * * * * * * * * * * * 例 编程，输入一个自然数，判断该数是否为素数，如果是素