浙江大学数值分析C语言编程习题.doc

C语言编程****题第二章****题2--40x2+9=0的所有实根。#include<>#include<>#defineN10000doubleA,B,C;doublef(doublex){return(A*x*x*x*x+B*x*x+C);}voidBM(doublea,doubleb,doubleeps1,doubleeps2){intk;doublex,xe; doublevaluea=f(a); doublevalueb=f(b); if(valuea>0&&valueb>0||valuea<0&&valueb<0)return; printf("Findingrootintherange:[%.3lf,%.3lf]\n",a,b);for(k=1;k<=N;k++){x=(a+b)/2;xe=(b-a)/2;if(fabs(xe)<eps2||fabs(f(x))<eps1){printf("Thexvalueis:%g\n",x);printf("f(x)=%g\n\n",f(x));return;}if(f(a)*f(x)<0)b=x;elsea=x;}printf("Noconvergence!\n");}intmain(){doublea,b,eps1,eps2,step,start;printf("PleaseinputA,B,C:\n");scanf("%lf%lf%lf",&A,&B,&C);printf("Pleaseinputa,b,step,eps1,eps2:\n");scanf("%lf%lf%lf%lf%lf",&a,&b,&step,&eps1,&eps2); for(start=a;(start+step)<=b;start+=step){ doubleleft=start; doubleright=start+step; BM(left,right,eps1,eps2); }return0;}运行:PleaseinputA,B,C:6-409Pleaseinputa,b,step,eps1,eps2:-101011e-51e-5Findingrootintherange:[-,-]Thexvalueis:-(x)=-:[-,]Thexvalueis:-(x)=:[,]Thexvalueis:(x)=:[,]Thexvalueis:(x)=-(fabs(xe)<eps2||fabs(f(x))<eps1)改为if(fabs(xe)<eps2&&fabs(f(x))<eps1)会提高精度,对同一题运行结果:Findingrootintherange:[-,-]Thexvalueis:-(x)=--007Findingrootintherange:[-,]Thexvalueis:-(x)=--006Findingrootintherange:[,]Thexvalueis:(x)=--006Findingrootintherange:[,]Thexvalueis:(x)=--007****题2-=(弧度)附近的根。#include<>#include<>#defineN100#(doublex0,doubleeps){intk;doublex;doublex1,x2;for(k=1;k<=N;k++){x1=sin(x0)/cos(x0);x2=sin(x1)/cos(x1);x=x0-(x1-x0)*(x1-x0)/(x2-2*x1+x0);if(fabs(x-x0)<eps){printf("Thexvalueis:%g\n",x);return;}x0=x;}printf("Noconvegence!\n");}intmain(){doubleeps,x0;printf("Pleaseinputeps,x0:\n");scanf("%lf%lf",&eps,&x0)