Calcul matriceal elementar

24 SEMINAR 1. CALCUL MATRICEAL ELEMENTAR21: end22: end23: for k=1:n24: for i=1:m25: for j=1:p26: C(i,j)=C(i,j)+A(i,k)*B(k,j);27: end28: end29: end1: function [ A ]=inmultire LU(L,U)2:3: %-------------------------------------------------------------------------4: % Algoritmul 1.115: % Functia calculeaza eficient produsul a doua matrice patratice L si U,6: % prima fiind inferior triunghiulara in timp ce a doua este superior7: % triunghiulara. Rezultatul este memorat in matricea A.8: % Apelul: [ A ]=inmultire LU(L,U)9: %10: % Buta Valentin, aprilie, 200611: %-------------------------------------------------------------------------12:13: [ n,x ]=size(L);14: [ m,y ]=size(U);15: if n∼=x16: error(’Matricea L nu este patratica’);17: end18: if m∼=y19: error(’Matricea U nu este patratica’);20: end21: if n∼=m22: error(’Matricele nu pot fi inmultite!’)23: end24: for i=1:n25: for j=1:i26: A(i,j)=0;27: for k=1:j28: A(i,j)=A(i,j)+L(i,k)*U(k,j);29: end30: end31: for j=i+1:n32: A(i,j)=0;33: for k=1:i34: A(i,j)=A(i,j)+L(i,k)*U(k,j);35: end36: end37: end