x=0:1:20;
niu=[0 1 2];J=besselj (niu,x);
plot(x,J) ;grid