Inteligencia Artificial - ALGORITMO CORDIC

Hola:
tengo que implementar La función sigmoide en cordic .he pensado en programar el algoritmo cordic para seno y coseno hiperbólicos para sacar la exp(-x). el tema es que el código funcióna solo el en el rango de convergencia y cuando lo expando no . a ver si alguien me puede echar una mano . el código esta implementado en scilab :

function [xn,yn,zn,k]= cos(a)
k=1;
k1=1;
k2=1;
y(1)=0;
z(1)=a;
l=4;
x1=0;
y1=0;
z1=0;
for j=-5:-1:0
k1 =k1*(sqrt(1-(1-(2^(2* j)/4))));
end
for j=1:1:100
k2=k2*(sqrt(1-(1/(2^(2* j)))));
end
k=k1*k2;

x(1)=1/k;
for i=-5:-1:0 // iteraciones adicionales
if z(i)<0
d=-1;
else
d=1;
end
x(i+1)=x(i)+y(i)*(d*(1-(2^i/4)));
y(i+1)=y(i)+x(i)*(d*(1-(2^i/4)));
z(i+1)=z(i)-d*atanh((1-(2^i/4)));
x1=x(i+1);
y1=y(i+1);
z1=z(i+1);

end


i=1;
x(1)=x1;
y(1)=y1;
z(1)=z1;
for i=1:1:100 // iteraciones normales
if z(i)<0
d=-1;
else
d=1;
end
if i+1==l
x(i+1)=x(i)+y(i)*(d/2^(i));
y(i+1)=y(i)+x(i)*(d/2^(i));
z(i+1)=z(i)-d*atanh(1/2^(i));
l=3*l+1;
end


x(i+1)=x(i)+y(i)*(d/2^(i));
y(i+1)=y(i)+x(i)*(d/2^(i));
z(i+1)=z(i)-d*atanh(1/2^(i));


end

gracias :)


if z==0.0001


break;
end

plot(x,y,'-',x,y,'b-')
n=i;
xn=x(i+1);
yn=y(i+1);
zn=z(i+1);
endfunction


gracias