em có đoạn function này nhưng không chạy được bác nào báo lỗi giùm em với:
Code:
function T = phantich(a)
% phan tich ma tran cholesky
n = size(A, 1);
for i=1:n
for j=1:n
if a(i,j)> a(j,i)
disp('ma tran khong doi xung ')
a= input('nhap lai ma tran a = ')
end
for j = 1:n
p = A(j, j) - dot(A(j, 1:j-1),A(j, 1:j-1));
if p < 0.0
disp('Ma tran xac dinh duong,vui long nhap lai')
a= input('nhap lai ma tran a = ')
end
for i=1:length(a)
if ( det( a(1:i, 1:i) ) <= 0 )
disp('ma tran khong xac dinh duong')
a= input('nhap lai ma tran a = ')
end
T=chol(a) %ham phan tich cholesky.đây cũng là chỗ nó báo sai
;
Comment