Как реализовать метод Р. Ковэю в Паскале?

[Korlet]

Для случая, когда m представляется степенью двойки, интересный квадратичный метод предложил Р. Ковэю:
X0 mod 4=2, Xn+1=Xn(Xn+1) mod 2e, n>=0.
Данный метод гарантирует достаточно большой период.
---
Как написать это в Паскале? вот я пробовал, но не выходит:
begin
for var j:=1 to n-1 do
x[j+1]:=(x[j]*(x[j+1])) mod t;
writeln('С использованием метода Ковэю:');
for var j:=1 to n-1 do
writeln(x[j]);
end;

t предлагается ввести пользователю заранее как любую степень двойки.