Минимизация логической функции Опции

[Tenshi]

День добрый, Уважаемые программисты
Прошу помощи в данном вопросе: "составить программу минимизации логической функции произвольной длины".
На данном этапе мне нужна теория и желательно алгоритм действий. У кого есть ссылки на источники или знания помогите

[Tenshi]

Что бы разобраться в принципе работы, просьба объясните (если можно комменты в код), что и какая функция и процедура здесь делает на этом примере (с данного форума), с логическими функциями еще не работал до этого момента (а будет еще курсовая). И еще вопрос: здесь происходит минимизация или просто вывод в таблицу? Заранее благодарю.

 program Minimization;
  function fromdec(n,m : longint): string;
	var
	  s: string;
	const
	  digit: string[16] = '0123456789ABCDEFQ';
  begin
	s := '';
	repeat
	s := digit[(n mod m) + 1] + s;
	n := n div m;
	until n = 0;
	while length(s) < 2 do s := '0' + s;
	fromdec := s;
  end;

  const
	s: string = 'A*/B+/A*/C*/D+B*/C*D';
  type
	tmatrix = array[0 .. 3, 0 .. 3] of 0 .. 1;
  var
	mx: tmatrix;
	p, i: integer;
	mask, sub_s: string;

  function matches(s, mask: string): boolean;
	var i: integer;
	begin
	  matches := true;
	  for i := 1 to length(s) do begin
		  if (mask[i] = 'X') or (s[i] = mask[i]) then {}
		  else matches := false

	end;
  end;

  procedure print_mx(const mx: tmatrix);
	 var i, j: integer;
	 begin
		for i := 0 to 3 do begin
		for j := 0 to 3 do write(mx[i, j]:3);
		writeln;
	 end;
  end;

  procedure fill_mx(var mx: tmatrix;
			first, second: string);
	 var i, j: integer;
	 begin
		for i := 0 to 3 do
		for j := 0 to 3 do
		if matches(fromdec(i, 2), first) and matches(fromdec(j, 2), second)
		then mx[i, j] := 1;
	 end;

begin
   write ('function: '); readln(s);
   s := s + '+';
   repeat
	 p := pos('+', s);
	 if p > 0 then begin
		sub_s := copy(s, 1, p-1);
		mask := '';
		for i := 1 to 4 do mask := mask + 'X';
		for i := 1 to length(sub_s) do begin
			if sub_s[i] in ['A' .. 'D'] then
			if (i > 1) and (sub_s[i - 1] = '/') then
			mask[ord(sub_s[i]) - ord('A') + 1] := '0'
			else
			mask[ord(sub_s[i]) - ord('A') + 1] := '1';
		end;
	 fill_mx(mx, copy(mask, 1, 2), copy(mask, 3, 2));
	 delete(s, 1, p);
	 end;
   until p = 0;
   print_mx(mx);
   readln;
end.