plot_solenoid(6)

f(2,10,2)
function []=f(R,h,N)
    full=2*pi*N;
    fi=0:0.1:2*full;
    x=R*cos(fi);
    y=R*sin(fi);
    z=h*fi/full;
    plot3(x,y,z);
end