plot[x*besselJ(1,x)/besselJ(0,x)] between 0 and 15