% LE TUBE HELICOIDAL

% Ce script crée une hélice en volume de A tours, de rayon de courbure B, de rayon du

% tube C, de hauteur D centrée en (E;F;G).

A = 3;

B = 4;

C = 1;

D = 16;

E = 0;

F = 0;

G = 0;

[t,p] = meshgrid(0:pi/25:A*2*pi,0:pi/25:2*pi);

x = (B+C*cos(p)).*cos(t);

y = (B+C*cos(p)).*sin(t);

z = D*t/(A*2*pi)+C*sin(p);

surf(x+E,y+F,z+G)

colormap jet

shading interp

axis image

% NOTE: Les équations sont les mêmes que pour le tore, à l'élévation et à l'azimut près.