Составить алгоритм, определяющий, является ли конечный неориентированный граф гамильтоновым (теорема Дирака)

program HamiltonCycle;
 
type
  TMatrixAdjacent = array of array of integer;
  TArray = array of integer;
 
{
Степень вершины (англ. degree, также валентность, англ. valency)
в теории графов — количество рёбер графа G, инцидентных вершине x.
При подсчёте степени ребро-петля учитывается дважды.
}
  {для неориентированного графа}
  function deg(v: integer; const M: TMatrixAdjacent; NVertex: integer): integer;
  var
    j: integer;
  begin
    deg := 0;
    for j := 0 to NVertex - 1 do
      if M[v, j] <> 0 then
        Inc(deg);
    if M[v, v] <> 0 then
      Inc(deg);
  end;
 
{
   Необходимое условие существования гамильтонова пути в неориентированном графе:
   если неориентированный граф G содержит гамильтонов цикл, тогда в нём не существует
   ни одной вершины x(i) с локальной степенью p(x(i)) < 2.
  }
  function HamiltonianProperty(const M: TMatrixAdjacent; NVertex: integer): boolean;
  var
    x, y: integer;
    Res:  boolean;
  begin
    Res := (NVertex >= 3);
    if Res then
    begin
      for x := 0 to NVertex - 1 do
        Res := Res and (deg(x, M, NVertex) >= 2);
    end;
    HamiltonianProperty := Res;
  end;
 
{
  условие Дирака существования гамильтонова пути:
  пусть p — число вершин в данном графе и p>=3;
  если степень каждой вершины не меньше, чем p/2, то данный граф — гамильтонов.
  Это достаточное условие не является необходимым.
  }
  function Dirac(const M: TMatrixAdjacent; NVertex: integer): boolean;
  var
    Res: boolean;
    v: integer;
  begin
    Res := (NVertex >= 3);
    if Res then
    begin
      for v := 0 to NVertex - 1 do
        Res := Res and (deg(v, M, NVertex) >= (NVertex div 2));
    end;
    Dirac := Res;
  end;
 
var
  MA: TMatrixAdjacent;
  NVertex: integer;
begin
...........заполнение матрицы смежности
  if HamiltonianProperty(MA, NVertex) then
    writeln('По необходимому условию, граф может быть гамильтоновым.')
  else
    writeln('По необходимому условию, граф не - гамильтонов');
  Write('Проверка графа по теореме Дирака.');
  if Dirac(MA, NVertex) then
    writeln('Условие Дирака выполнилось. Граф - гамильтонов.')
  else
    writeln('Условие Дирака не выполнилось. Граф может быть гамильтоновым.');
end.