Pascal/Turbo Pascal - Ayuda por favor

Me estoy volviendo loca.Si alguien sabe como hacer el método de Newton-Raphson para la ecuación z^2+c (para números complejos) en Turbopascal por favor que me la mande o se ponga en contacto conmigo si quereis cobrar adelante estoy dispuesta a todo.

uses crt;
var
x1:real;
acc:integer;

function f(t:real):real;
begin
f:=t*sin(t)+cos(t);
end;

function fd(t:real):real;
begin
fd:=t*cos(t);
end;

begin
clrscr;
writeln('Valor Inicial?: ');
readln(x1);

if fd(x1)=0 then
writeln('Valor Inicial no válido')
else
begin
writeln('A cuántos decimales? ');
readln(acc);
while (abs(f(x1))>1/exp(acc*(ln(10)))) do
begin
x1:=x1-f(x1)/fd(x1);
end;
writeln('Solución: ',x1:5:acc+1);
end;
readkey;
end.

Este algoritmo no me pertenece pero olvidé de dónde lo saqué :(.

gracias pero lo necesito con números complejos.Igual tambien lo sabes.De todas formas muchas gracias.