Tugas Pemodelan No 6

function [flag] = pdtest(P)

	% Mencari nilai N, yaitu besarnya matrix P
	% Matrix P harus merupakan rectangular matrix
	[row, column] = size(P);
	if row ~= column
		error('Matrix harus rectangular');
	end
	N = row;

	% Menguji sifat positive definite matrix P dengan Sylvester's criterion
	% Yaitu dengan menguji determinan submatrix 1x1, 2x2, ... NxN
	order = 0;
	for iterDet = 1:N
		if det(P(1:iterDet, 1:iterDet)) > 0
			order = order + 1;
		end
	end

	% Apabila jumlah determinan submatrix yang positif sama dengan N,
	% maka matriks P adalah matriks positive definite 
	if order == N
		flag = 1;
	else
		flag = 0;
	end
end