200909014110陈伟数值分析.docxVIP

华 北 科 技 学 院 上 机 报 告 系（部） 基础部 专业、班级 计算B091 姓名 陈伟 学号 200909014110 课程名称 数值分析 上机题目 数值积分、解线性方程组的直接法 任课教师 李慧 指导教师 李慧 成绩（优、良、中、及格、不及格） 华 北 科 技 学 院 基 础 部 1.编写复化柯特斯求积公式的程序，并计算例1和例2，观察n为多少时与6位有效数字。 #include stdio.h #include conio.h void main() { int i,n=2; float s; float f(float); float Simpson(float(*)(float),float,float,int); for(i=0;i=10;i++) { s=Simpson(f,0,1,n); printf(s(%d)=%f\n,n,s); n*=2; } } float Simpson(float (*f)(float),float a,float b, int n) { int k; float s,s1=0,s2,s3; float h=(b-a)/n; s3=f(a+h/2); s2=(f(a+h/4)+f(a+h*3/4)); for(k=1;k=n-1;k++) { s1+=f(a+k*h); s2+=(f(a+k*h+h/4)+f(a+k*h+h*3/4)); s3+=f(a+k*h+h/2); } s=h/90*(7*f(a)+14*s1+32*s2+12*s3+7*f(b)); return s; } float f(float x) { return 1/(1+x*x); } 结果 s(2)=0.785399 s(4)=0.785398 s(8)=0.785398 s(16)=0.785398 s(32)=0.785398 s(64)=0.785398 s(128)=0.785398 s(256)=0.785398 s(512)=0.785398 s(1024)=0.785398 s(2048)=0.785398 Press any key to continue 2.用龙贝格方法上机计算： #include stdio.h #include conio.h #include math.h float f(float x) { return sin(2*cos(x))*sin(x)*sin(x); } float Romberg(float a,float b,float (*f)(float),float epsilon) { int n=1,k; float h=b-a,x,temp; float T1,T2,S1,S2,C1,C2,R1,R2; T1=(b-a)/2*((*f)(a)+(*f)(b)); while(1) { temp=0; for(k=0;k=n-1;k++) { x=a+k*h+h/2; temp+=(*f)(x); } T2=(T1+temp*h)/2; if(fabs(T2-T1)epsilon) return T2; S2=T2+(T2-T1)/3.0; if(n==1){T1=T2;S1=S2;h/=2;n*=2;continue;} C2=S2+(S2-S1)/15; if(n==2){C1=C2;T1=T2;S1=S2;h/=2;n*=2;continue;} R2=C2+(C2-C1)/63; if(n==4){R1=R2;C1=C2;T1=T2;S1=S2;h/=2;n*=2;continue;} if(fabs(R2-R1)epsilon) return R2; R1=R2;C1=C2;T1=T2;S1=S2;h/=2;n*=2; } } void main() { float epsilon=5E-6; float t=3.1415926/2; printf(R=%f\n,Romberg(0,t,f,epsilon)); } 结果 R=0.507967 Press any key to continue 用列主元消去法解方程组: #include stdio.h