求矩阵用改进的Cholesky分解法解方程组用改进的Cholesky分解法解方程组请帮忙给出正确答案

正确答案：由改进的Cholesky分解(3．1．7)式对i＝1234计算得i＝1时 d₁＝a₁₁＝4i＝2时 c₂₁＝a₂₁＝－2 l₂₁＝c₂₁／d₁＝－0．5 d₂＝a₂₂－c₂₁l₂₁＝5i＝3时 c₃₂＝a₃₂－c₃₁l₂₁＝－1．5 l₃₁＝c₃₁／d₁＝0．25l₃₂＝c₃₂／d₂＝－0．3 d₃＝a₃₃－c₃₁l₃₁－c₃₂l₃₂＝5．3i=4时 c₄₁＝a₄₁＝0c₄₂＝a₄₂－c₄₁l₂₁＝1c₄₃＝a₄₃－c₄₁l₃₁－c₄₂l₃₂＝－1．7 l₄₁＝c₄₁／d₁＝0l₄₂＝c₄₂／d₂＝0．2l₄₃＝c₄₃／d₃＝－0．320754717 d₄＝a₄₄－c₄₁l₄₁－c₄₂l₄₂－c₄₃l₄₃＝3．254716981 D＝diag(455．33．254716981) 由Ly＝b(b＝(3333)^T)解得 y₁＝b₁＝3 y₂＝b₂－l₂₁y₁＝4．5 y₃＝b₃－l₃₁y₁－l₃₂y₂＝3．6 y₄＝b₄－l₄₁y₁－l₄₂y₂－l₄₃y₃＝3．254716981 再由L^Tz＝D^-1y解得 χ₄＝y₄／d₄＝1 χ₃＝y₃／d₃－l₄₃χ₄＝1 χ₂＝y₂／d₂－l₃₂χ₃－l₄₂l₄＝1 χ₁＝y₁／d₁－l₂₁χ₂－l₃₁χ₃－l₄₁l₄＝1 故解得原方程组的解χ＝(1111)^T。
由改进的Cholesky分解(3．1．7)式,对i＝1,2,3,4计算得i＝1时,d1＝a11＝4i＝2时,c21＝a21＝－2l21＝c21／d1＝－0．5d2＝a22－c21l21＝5i＝3时,c32＝a32－c31l21＝－1．5l31＝c31／d1＝0．25,l32＝c32／d2＝－0．3d3＝a33－c31l31－c32l32＝5．3i=4时,c41＝a41＝0,c42＝a42－c41l21＝1,c43＝a43－c41l31－c42l32＝－1．7l41＝c41／d1＝0,l42＝c42／d2＝0．2,l43＝c43／d3＝－0．320754717d4＝a44－c41l41－c42l42－c43l43＝3．254716981D＝diag(4,5,5．3,3．254716981)由Ly＝b(b＝(3,3,3,3)T)解得y1＝b1＝3y2＝b2－l21y1＝4．5y3＝b3－l31y1－l32y2＝3．6y4＝b4－l41y1－l42y2－l43y3＝3．254716981再由LTz＝D-1y解得χ4＝y4／d4＝1χ3＝y3／d3－l43χ4＝1χ2＝y2／d2－l32χ3－l42l4＝1χ1＝y1／d1－l21χ2－l31χ3－l41l4＝1故解得原方程组的解χ＝(1,1,1,1)T。