axis([0 4 0 10])
x=0:.1:4;
for k = 0:2
    y= besseli(k,x);
    hold on;   
    if k == 0
        plot(x,y,'k');
    end
    if k == 1
        plot(x,y,'k:');
    end
    if k == 2
        plot(x,y,'k--');
    end
end
