C语言程序设计（第三版） 教学课件 作者 恰汗.合孜尔 第6章 函数与编译预处理.pptVIP

例6.35 用牛顿迭代法求方程-3x3+4x2-5x+6=0在1.0附近的实根，精度要求为 。 牛顿迭代公式如下所示： 只要相邻两次迭代值误差的绝对值小于等于给定的精度要求，即| xk-xk-1|≤ 满足，我们就认为最后一次的迭代值为方程的近似解。实现牛顿迭代法的基本步骤如下： (1) 给出初始近似根x0及精度 ； (2) 计算 (3) 若| xk-xk-1|≤ ，转向(4)；否则x0=x1，转向(2)； (4) 输出满足精度的根x1，结束。 6.8 程序举例 图6-5 牛顿迭代法的流程图和N-S图描述 6.8 程序举例 当f(x)=-3x3+4x2-5x+6时，有f‘(x)=-9x2+8x-5 ，设初值x0=1.0。 #include stdio.h #include math.h #define eps 0.00001 /* 宏定义容许误差eps为0.00001 */ float f(float x) /* 定义函数f(x) */ { return(((-3*x+4)*x-5)*x+6); } float f1(float x) /* 定义函数f(x)的导数f1(x) */ { return((-9*x+8)*x-5); } void main() { float x0,x1=1.0; do { x0=x1; /* 准备下一次迭代的初值 */ x1=x0-f(x0)/f1(x0); /* 牛顿迭代 */ }while(fabs(x1-x0)eps); /* 若满足精度，输出近似根 */ printf(x=%f\n,x1); } 程序运行结果： x=1.265328 6.8 程序举例 例 6.36 使用梯形求积公式求下列定积分的值。 梯形积分法的基本思想是，将区间[a,b]分成n个相等的小区间，则每个小区间的长度为h=(b-a)/n，对每个小区间均实施如下的梯形求积： 将这些小梯形的求积值加起来，可以得到如下梯形求积公式: 6.8 程序举例 基于上述分析，实现梯形积分法的基本步骤： (1)输入区间[a,b]的端点a,b值以及分割数n； (2)将区间[a,b]等分成n个小区间,每一个小区间的长度h=(b-a)/n； (3)计算每一个等分点的函数值与yi=f(a+ih)(i=0,1,2,…,n)； (4)计算: (5) 输出s的值，结束。 6.8 程序举例 图6-6 梯形积分法的N-S图描述 6.8 程序举例 #include stdio.h #define N 16 /* 宏定义等分数N为16 */ /* 定义函数func，计算被积函数的值 */ float func(float x) { float y; y=4.0/(1+x*x); return(y); } /* 定义函数gedianzhi，计算等分点的函数值 */ void gedianzhi(float y[],float a,float h) { int i; for(i=0;i=N;i++) y[i]=func(a+i*h); } 6.8 程序举例 /* 定义函数trapeze，计算梯形面积 */ float trapeze(float y[],float h) { float s; int i; s=(y[0]+y[N])/2.0; for(i=1;iN;i++) s+=y[i]; return(s*h); } 6.8 程序举例 void main() { float a,b,h,s,f[N+1]; printf(input a,b=); scanf(%f,%f,a,b); h=(b-a)/N; gedianzhi(f,a,h); /* 调用函数gedianzhi，计算等分点的函数值 */ s=trapeze(f,h); /* 调用函数trapeze，计算梯形面积 */ printf(s=%f\n