数值分析答案.docVIP

数值分析答案 本人利用课余时间做了一下数值分析的题目，只有2-4题，第1题设计数组太复杂，用C编不出来，谁会用MATLAB的同学请做第1题。 声明： 答案仅供参考，如有考试出现雷同，与本人无关！ ————刘洪波 第二题： #includeiostream #includemath.h using namespace std; double f(double x) { double y; y=1-x-sin(x); return (y); } void main() { float a=0,b=1; float A,D=b-a; while(D1e-4) { A=(a+b)/2; if(f((a+b)/2)*f(b)0) a=A; else if(f((a+b)/2)*f(b)0) b=A; D=b-a; } cout方程的根为:\nA\n; } 结果： 第三题： #includecmath #includeiostream using namespace std; double f(double x) { return x==0 ? 1 : 1/(1+x); } double S( double a, double b, int n) { double h = (b - a) / (2.0 * n); double x, F0, F1, F2; F0 = f(a) + f(b); F1 = 0; F2 = 0; for (int j=1; j=2*n-1; ++j) { x = a + j * h; if (j % 2 == 0) F2 += f(x); else F1 += f(x); } return h*(F0 +2*F2 +4*F1)/3.0 ; } void main() { int i = 1; double xx=0, sx=1, r1, r2; r1 = S(xx, sx, i); do { r2 = r1; r1 = S(xx, sx, ++i); } while (fabs(r1-r2) 1e-6); cout 积分结果： r1 endl; } 结果： 第四题： #includestdio.h #define Z y*y double f(double x,double y) { return Z; } void R_K(double y0,double x0,double h,double b) { double y1;double R[4];int i; for (i=0;i=4;i++) { R[0]=f(x0,y0); R[1]=f(x0+h/2.0,y0+h*R[0]/2.0); R[2]=f(x0+h/2.0,y0+h*R[1]/2.0); R[3]=f(x0+h,y0+h*R[2]); y1=y0+h*(R[0]+2*R[1]+2*R[2]+R[3])/6.0; y0=y1;x0=x0+h; printf(%d\t%f\t%f\n,i,x0,y1); } } void main() { R_K(-1.0,1.0,0.1,1.5); } 结果：