clf; clear
format short e
a = [ -1 -1 -1 1 1 1 2 2 2]
b = [-.1 .1 .2 -.1 .1 .2 -.1 .1 .2]

t = linspace(1,50,1000)
for i = 1:9
    for j = 1:1000
    x(i,j) = a(i).*cos(t(j));
    y(i,j) = a(i).*sin(t(j));
    z(i,j) = b(i).*t(j);
    end
end
for i = 1:9
subplot(3,3,i)
plot3(x(i,:), y(i,:), z(i,:))
grid on
axis([-1.5 1.5 -1.5 1.5 -10 15])
if i == 1
title(sprintf('a = %0.1f b = %0.1f (cl176)', a(i), b(i)))
end
if i > 1
title(sprintf('a = %0.1f b = %0.1f', a(i), b(i)))
end
end

print -deps test