help с задачкой

[SYL@R]

Даны натуральные числа р и q. Получить все делители числа q, взаимно простые с р.

[eoln]

var
  p, q: word;
  i, j: word;
  b: boolean;
begin
  p := 8;
  q := 60;
  for i := 2 to q do
  begin
    if q mod i = 0 then //если найден делитель числа q
    begin
      b := true;
      write(i:3); //выводим его на экран
      for j := 2 to i do
        if ((i mod j = 0) //пытаемся найти общий делитель
          and (i <> j)) //не равный 1 и не равный самому числу
          and (p mod j = 0)
          then b := false; //если нашли такой делитель
      if b then writeln('v':2) else writeln //ставим галочку если удовл условию задачи
    end
  end;
  readln
end.

[puporev]

Числа называются взаимно простыми, если их наибольший общий делитель равен 1. Поэтому используем функцию NOD(m,n);

uses crt;
var p,q,i:longint;
    k:byte;
function NOD(m,n:longint):longint;
begin
while m<>n do
if m>n then m:=m-n else n:=n-m;
NOD:=m;
end;
begin
clrscr;
write('p=');readln(p);
write('q=');readln(q);
for i:=2 to trunc(sqrt(q))do //перебирам до корня квадратного из q.
if (q mod i=0)and(NOD(i,p)=1)then
   begin
     writeln(i);
     k:=k+1;//считаем делители, подходящие по условию
   end;
if k=0 then writeln('Takih delitelej net!');//если к=0, выводим сообщение
readln
end.