You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
12 SEMINAR 6. EXAMINARE10: %-------------------------------------------------------------------------11:12: [ n,m ]=size(A);13: if n∼=m14: error(’Matricea nu este patratica!’);15: end16: for j=2:n17: if A(1,j)∼=018: error(’Matricea nu este inferior triunghiulara!’);19: end20: end21: for i=2:n22: for j=i+1:n23: if A(i,j)∼=024: error(’Matricea nu este inferior triunghiulara!’);25: end26: if A(j,i)∼=027: error(’Matricea are elemente nenule sub diagonala care nu se afla pe prima coloana!’);28: end29: end30: end31:32: X(1,1)=1/A(1,1);33: for i=2:n34: X(i,1)=-A(i,1)/(A(1,1)*A(i,i));35: end36: for j=2:n37: for i=1:n38: X(i,j)=0;39: end40: X(j,j)=1/A(j,j);41: end1: function [ L,U ]=Crout MH(H)2: %-------------------------------------------------------------------------3: % Seminar 6 Problema 2 punctul a.4: % Algoritmul calculeaza factorizarea Crout fara pivotare in cazul unei5: % matrice patratice inferior Hessenberg cu toate submatricele li<strong>de</strong>r principal6: % nesingulare.7: % Apelul: [ U,L ] = Crout MH(H)8: %9: % Valentin Buta, mai, 200610: %-------------------------------------------------------------------------11:12: [ n,m ]=size(H);13: if n∼=m14: error(’Matricea nu este patratica!’);15: end16: for i=1:n-217: for j=i+2:n