东南大学数值分析上机练习前三章 -.docVIP

数值分析上机练习 （以VC++6.0为操作平台） 第一章 20.（1）程序如下： #includeiostream.h #includemath.h #includeiomanip.h void main() { int N; double sum; float sum1=0, f; cout请输入N的初值：; cinN; for(float i=2;i=N;i++) { f=1/(i*i-1); sum1=sum1+f; } coutsum1=sum1endl; sum=0.5*(1.5‐1.0/N‐1.0/(N+1)); coutsum=sumendl; } 20．（2）程序如下： #includeiostream.h #includemath.h #includeiomanip.h void main() { int N; float sum2=0, f; double sum; cout请输入N的初值：; cinN; for(float i=N;i=2;i--) { f=1/(i*i-1); sum2=sum2+f; } coutsum2=sum2endl; sum=0.5*(1.5‐1.0/N‐1.0/(N+1)); coutsum=sumendl; } 20.(1)运行结果如下： 100 sum1=0.740049 sum=0.74005 Press any key to continue 10000 sum1=0.749852 sum=0.7499 Press any key to continue 1000000 sum1=0.749999 sum=0.749999 Press any key to continue 20.(2)运行结果如下： 100 sum2=0.74005 sum=0.74005 Press any key to continue 10000 sum2=0.7499 sum=0.7499 Press any key to continue 1000000 sum2=0.749999 sum=0.749999 Press any key to continue 结果分析： 当n=100时，sum1的6位数字中有5位有效数字，而sum2所有数字都是有效数字。 当n=10000时，sum1的6位数字中有4位有效数字，而sum2所有数字都是有效数字。 当n=1000000时sum1，sum2有相同位数的有效数字。 由运行结果可知：当精度较低时按从大到小的顺序累加的精度低于按从小到大的顺序累加结果的精度。至于当精度较高时，sum1，sum2会得出相同结果。 第二章 20.程序如下： #includeiostream.h #includemath.h #includeiomanip.h void main() { int count=0; double fl,f2,y; double x,xl,u=1.0; cout请输入x0的初值:; cinxl; cout请输入容许误差y的初值:; ciny; for(count=0;fabs(u)=y;count++) {x=xl; fl=x*x*x/3-x; f2=x*x-1; xl=x-fl/f2; u=xl-x; if(count10000) { cout迭代结果不收敛于0!endl; break; } } cout方程的一个根为:xlendl; } 运行结果如下： 请输入x0的初值：100 请输入容许误差y的初值:1e-6 方程的一个根为：1.73205 Press any key to continue 请输入x0的初值：0.7745966 请输入容许误差y的初值:1e-6 方程的一个根为：0 Press any key to continue 请输入x0的初值：0.7745967 请输入容许误差y的初值:1e-6 方程的一个根为：-1.73205 Press any key to continue 请输入x0的初值：-0.7745966 请输入容许误差y的初值:1e-6 方程的一个根为：0 Press any key to continue 请输入x0的初值：-0.7745967 请输入容许误差y的初值:1e-6 方程的一个根为：1.73205 Press any key to continue 请输入x0的初值：-100 请输入容许误差y的初值:1e-6 方程的一个根为：-1.73205 Press any key to continue 其中尽可能大的=0.7745966是经过多次试算得出的结论，并得出以下结