资源描述：
《【精品】数值分析编程》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、习题一：已知A=r12.384122」15237-1.0610741.112336-0.1135840.7187191.7423823.067813一2.031743、2」1523719.141823-3.125432-1.0123452.1897361.836742-0.7841651.1123483.123124-1.061074-3.12543215.5679143.1238482.031454-3.741856-1.0567810.336993-1.0101031.112336-1.0123453.
2、12384827.1084374.101011-3.7418562.101023-0.71828-0.037585-0.1135842.1897362.0314544.10101119.8979180.431637-3.1112232.1213141.7843170.7187191.5638491.836742-3.7418560.4316379.789365-0.103458-1.1034560.2384171.742382-0.784165-1.0567812.101023-3.111223-0.103
3、45814.71384653.123789-2.2134743.0678131.1123840.336993-0.718282.121314-1.1034563.12378930.7193344.446782厂2.0317433.123124-1.010103-0.0375851.7843170.238417-2.2134744.44678240.00001丿b=(2.1874369,33.992318,-25.173417,0.84671695,1.784317,-86.612343,1.1101230,
4、4.719345,-5.6784392)；1.用Household变换，把A化为三对角阵B。(并打印E)。算法：(1)(1)(1)令4()=A,atj二伽，已知A—,?!一］二(“”•)。且(r)1/(2)Sr=(£血)2)%i=r+l(r)(3)%=Sr+ar+lrSrr(r)(r)(r)(r)V(4)ur=0,0,0…,0,ar+Lr+sign(ar+.r)sr,ar+u，…,arhr(6)qr=yr-krur(7)Ar=Ar_{-(qriirT+乞％)“=1,2,・・・/一2C源程序：#include
5、#include#defineNmainO{floata[N][N]={{12.38412,2.115237厂1.061074,1.112336厂0.113584,0.718719,1.742382,3.067813,-2.031743},{2.115237,19.141823,-3.125432-1.012345,2.189736,1.563849,-0.784165,1.112348,3.123124},{-1.061074厂3.125432,15.567914,
6、3.12384&2.031454,1.836742,-1.056781,0.336993厂1.010103},{1.112336厂1.012345,3.123848,27.108437,4.101011,-3.741856,2.101023,-0.71828,-0.037585},{-0.113584,2.189736,2.031454,4.101011,19.89791&0.431637,-3.111223,2.121314,1.784317},{0.718719,1.563849,1.836742厂3.
7、741856,0.431637,9.789365厂0.103458,-1.103456,0.238417},{1.742382,-0.784165,-1.056781,2.101023,-3.111223,-0.103458,14.7138465,3.123789,-2.213474},{3.067813,1.112348,0.336993,-0.7182&2.121314,-1.103456,3.123789,30.719334,4.446782},{-2.031743,3.123124,-1.01010
8、3,-0.037585,1.784317,0.238417,-2.213474,4.446782,40.00001}}；inti,j,r,n=N；/*分别定义变量*/floatm=0；floats,p,k=0；floatu[N],y[N],q[N]；for（r=l；r<=n-2；r++）/*求S的值*/{for（i=r+l；i〈=N；i++）m=m+a[i~l][r-1]*a[i-l][r-1]；s=sqrt（m
