《计算方法与实》书后实验题.docVIP

计算方法与实习 实习报告 2010-2011学年第三学期 院系：xxx 学号：xxx 姓名：xxx 日期：xxx 目录 实习题1—舍入误差与数值稳定………………………………3 实习题2—方程求根……………………………………………7 实习题3—线性方程组数值解法………………………………8 实习题4—插值法………………………………………………16 实习题5—曲线拟合……………………………………………17 实习题6—数值积分……………………………………………21 实习题1 P215/1.用两种不同的顺序计算,分析其误差的变化。 思路分析 用一个循环语句，对从1到10000进行叠加，两种不同顺序指从1叠加到10000和从10000叠加到1，每隔一定的叠加次数就比较一次误差。 用C++语言编程 （1）从1叠加到10000源代码如下： #includeiostream #includecmath #include iomanip using namespace std; int main(){ double N=10000,i=0; int a; double n=0,S=1.644834; for(i=1;i=N;i++){ n+=1/(i*i); a=i; if(a%500==0)coutsetprecision(7)迭代a次时，和为：S=n\t误差为：i=fabs(S-n)endlendl; } return 0; } 运行结果如下： 迭代500次时，和为：S=1.642936 误差为：i=0.001897934 迭代1000次时，和为：S=1.643935 误差为：i=0.0008994333 迭代1500次时，和为：S=1.644268 误差为：i=0.0005663776 迭代2000次时，和为：S=1.644434 误差为：i=0.0003998082 迭代2500次时，和为：S=1.644534 误差为：i=0.0002998532 迭代3000次时，和为：S=1.644601 误差为：i=0.0002332109 迭代3500次时，和为：S=1.644648 误差为：i=0.0001856066 迭代4000次时，和为：S=1.644684 误差为：i=0.0001499019 迭代4500次时，和为：S=1.644712 误差为：i=0.0001221307 迭代5000次时，和为：S=1.644734 误差为：i=9.991315e-005 迭代5500次时，和为：S=1.644752 误差为：i=8.173481e-005 迭代6000次时，和为：S=1.644767 误差为：i=6.658593e-005 迭代6500次时，和为：S=1.64478 误差为：i=5.376747e-005 迭代7000次时，和为：S=1.644791 误差为：i=4.278009e-005 迭代7500次时，和为：S=1.644801 误差为：i=3.32576e-005 迭代8000次时，和为：S=1.644809 误差为：i=2.492534e-005 迭代8500次时，和为：S=1.644816 误差为：i=1.757329e-005 迭代9000次时，和为：S=1.644823 误差为：i=1.103809e-005 迭代9500次时，和为：S=1.644829 误差为：i=5.19077e-006 迭代10000次时，和为：S=1.644834 误差为：i=7.184807e-008 Press any key to continue （2）从10000叠加到1的源代码如下： #includeiostream #includecmath #include iomanip using namespace std; int main(){ double N=10000,i=0; int a; double n=0,S=1.644834; for(i=N;i=1;i--){ n+=1/(i*i); a=i; if(a==N)cout以下为n大于200时每隔500所积累的和：endlendl; if(a%500==0a=200)coutsetprecision(7)n=a时，和为：S=n\t误差为：i=fabs(S-n)endlendl; if(a==200)coutendl以下为n小于200时每隔20所积累的和endlendl; if(a200a%10==0||a=10)coutsetprecision(7)n=a时，和为：S=n\t误差为：i=fabs(S-n)endlendl; } return 0; } 运行结果如下：