Actually I had this script beacuse an accident, when I was trying to make an equipotential line in order to modeling a geoelectrical method. but I found that my fail in this script can make it looks like a multilevel helix :D

phi = 3.14;
k=1;
for m = n:n:r;
    i(k)=m;
    k=k+1;    
end
V = ((I*rho)/(4*phi))*(1./i);
a = [0:1:360];
for p=1:length(i);
    x(:,p) = i(p)*cos(a.*pi/180);
    y(:,p) = i(p)*sin(a.*pi/180);
    z(:,p) = linspace(0,100,length(a));
end
figure 
for nl=1:length(i);
    plot3(x(:,:,nl),y(:,:,nl),z(:,:,nl))
    %text(x(1,nl),y(1,nl),[num2str(V(nl),2)])
    %title(‘Nilai equipotential pada tiap lintasan’)
    hold on
    grid on
end

phi = 3.14;
k=1;for m = n:n:r;    i(k)=m;    k=k+1;    end
V = ((I*rho)/(4*phi))*(1./i);
a = [0:1:360];
for p=1:length(i);    x(:,p) = i(p)*cos(a.*pi/180);    y(:,p) = i(p)*sin(a.*pi/180);    z(:,p) = linspace(0,100,length(a));end
figure for nl=1:length(i);    plot3(x(:,:,nl),y(:,:,nl),z(:,:,nl))  

 %text(x(1,nl),y(1,nl),[num2str(V(nl),2)])    

%title(‘Nilai equipotential pada tiap lintasan’)    hold on    grid on

end