using Plots;
gr()
function helicoid(r)
n = 100
θ = range(-2pi, 2π; length = n)
u = range(-10, 10; length = n)
x = r*u*cos.(θ)'
y = r*u*sin.(θ)'
z = ones(size(u))*θ'
return x, y, z
end
my_cg = cgrad([:red,:blue])
surface(helicoid(3), label="helicoid")