function helico(lambdar,lambdaz,n,dt,r0,z0,rc)

% t varie de 0 jusqu'à l'ouverture du coquillage
t=zeros(1,n);
r=zeros(size(t));	% rayon de l'hélico-spirale
z=zeros(size(t));	% déplacement le long des z
r(1)=r0;
z(1)=z0;

% calcule des coordonnées de l'hélico-spirale
for i=2:n
	t(i)=t(i-1)+dt;
	r(i)=r(i-1)*lambdar;
	z(i)=z(i-1)*lambdaz;
end

% courbe génératrice (un cercle si absent)
if nargin==6
	rc=linspace(r(1),max(r),n);
end

% représentation
% plot3(r.*cos(t),r.*sin(t),z)
% return
ver3(r.*cos(t),r.*sin(t),z,rc,24)
axis off
rotate3d on
