AXIS([0 20 -1 1])
x=0:0.01:20;
for k = 0:2
    		y= BESSELJ(k,x);
    		hold on;  
    	    plot(x,y,'k');
end
