Calcul matriceal elementar

1.7. ANEXA B 2133: end34: end35: end1: function [ c ]=GAXPY1(A,b,c)2:3: %-------------------------------------------------------------------------4: % Algoritmul 1.65: % Functia actualizeaza vectorul (coloana) c cu produsul unei matrice6: % A cu un vector (coloana) b. Algoritmul este cunoscut sub numele de GAXPY7: % pe linii.8: % Apelul: [ c ]=GAXPY1(A,b)9: %10: % Buta Valentin, aprilie, 200611: %-------------------------------------------------------------------------12:13: [ m,n ]=size(A);14: [ mb,nc ]=size(b);15: if mb∼=116: error(’b nu este vector coloana’);17: end18: if nc∼=n19: error(’A si b nu pot fi inmultite’);20: end21: [ mc,nc ]=size(c);22: if mc∼=123: error(’c nu este vector coloana’);24: end25: for i=1:m26: for j=1:n27: c(i)=c(i)+A(i,j)*b(j);28: end29: end1: function [ c ]=GAXPY2(A,b,c)2:3: %-------------------------------------------------------------------------4: % Algoritmul 1.75: % Functia actualizeaza vectorul (coloana) c cu produsul unei matrice6: % A cu un vector (coloana) b. Algoritmul este cunoscut sub numele de GAXPY7: % pe coloane.8: % Apelul: [ c ]=GAXPY2(A,b)9: %10: % Buta Valentin, aprilie, 200611: %-------------------------------------------------------------------------12:13: [ m,n ]=size(A);14: [ mb,nc ]=size(b);