% coquillage tubulaire

lambdar=1.011;
lambdaz=lambdar;
n=120;
dt=4*pi/180;
r0=0.01;
z0=0;
helico(lambdar,lambdaz,n,dt,r0,z0)
axis equal
view(-120,35)
colormap(bone)
pause

% coquillage turbiné
lambdar=1.007;
lambdaz=lambdar;
n=200;
dt=10*pi/180;
r0=0.1;
z0=1.9;
rc=linspace(0,1,n);
helico(lambdar,lambdaz,n,dt,r0,z0,rc)
colormap copper
view(-120,0)
pause

% coquillage tubulaire
lambdar=1.011;
lambdaz=lambdar;
n=200;
dt=10*pi/180;
r0=0.001;
z0=0.01;
rc=linspace(0,0.02,n);
helico(lambdar,lambdaz,n,dt,r0,z0,rc)
colormap copper
view(-120,0)
pause
