function L=CholeskyFelbontas(A) % Elkészíti egy pozitív definit szimmetrikus mátrix Cholesky-felbontását, % azaz azon L alsó háromszög mátrixot, melyre L*L^T=A. if A(1,1)<0 error('A mátrix nem pozitív definit'); end L=zeros(size(A)); L(1,1)=sqrt(1,1); for i=1:size(A,1) for j=1:i-1 L(i,j)=(A(i,j)-dot(L(i,1:j-1),L(j,1:j-1)))/L(j,j); s=A(i,i)-dot(L(i,1:i-1),L(i,1:i-1)); if s<0 error('A mátrix nem pozitív definit'); else L(i,i)=sqrt(s); end end end
Not enough input arguments. Error in CholeskyFelbontas (line 5) if A(1,1)<0