Newton法、一般迭代法Steffensen法、弦截法C语言代码.doc

一、#includemath.h #includestdio.h double f(double x) { return (3*x*x-exp(x)); } double f1(double x) { return (6*x-exp(x)); } void main() { double x1=1,x; do { x=x1; x1=x-f(x)/f1(x); printf(x=%.9lf\n,x1); }while(fabs(x1-x)0.000005); } 说明：f 原函数，f1f的导函数x1为初始值，通过x=%.9lf“控制输入输出格式 、#include stdio.h #include math.h int main() { double x=1,x1; while(1) { x1=pow(3*x+1,0.2); printf(x=%.6lf\n,x1); if(fabs(x1-x)0.000005 ) break; x=x1; } return 0; } 说明：x1初始值x1=pow(3*x+1,0.2);为迭代格式，0.000005误差通过x=%.6lf“控制输入输出格式 、#includestdio.h #includemath.h #define phi(x) pow(3*(x)+1,0.2); void main() { double x,x0,del,y,z; printf(x0=); scanf(%lf,x0); printf(\ndel=:); scanf(%lf,del); while(1) { y=phi(x0); z=phi(y); x=x0-(y-x0)*(y-x0)/(z-2*y+x0); printf(\n%.6lf,x); if(fabs(x-x0)del) break; x0=x; } } 说明：x0初始值，pow(3*(x)+1,0.2);φ（x）格式del为允许误差通过x=%.6lf“控制输入输出格式、#includemath.h #includestdio.h double f(double x) { //计算f(x)的值 return pow(2,x)+pow(3,x)-pow(4,x); } double point(double x1,double x2) { //计算与x轴交点的x值 printf(x=%.5f\n,(x1*f(x2)-x2*f(x1))/(f(x2)-f(x1))); return (x1*f(x2)-x2*f(x1))/(f(x2)-f(x1)); } int main() { //输入两个数x1,x2 double x1,x2,x; do { printf(输入两个数x1,x2:); scanf(%lf%lf,x1,x2); }while (f(x1)*f(x2)= 0); // 当输入两个数大于0为真时，继续重新输入 //关键循环步骤： do { x=point(x1,x2);//得到交点的值 if(f(x)*f(x1)0) x1=x;//新的x1 else x2=x; }while (fabs(f(x)) 0.000005); }