11.07.2015 Views

Probleme de examen

Probleme de examen

Probleme de examen

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!