x <- seq(from=0,to=15,len=100)
plot(x,bessel_J0(x),xlim=c(0,15),ylim=c(-0.8,1),type="l",xaxt="n",yaxt="n",bty="n",xlab="",ylab="",main="Figure 9.1, p359")
axis(1,pos=0)
axis(2,pos=0)
lines(x,bessel_Y0(x))
lines(x,bessel_J1(x),lty=2)
lines(x,bessel_Y1(x),lty=2)
x <- seq(from=0,to=13,len=100)
matplot(x,t(bessel_jl_array(3,x)),col="black",type="l",xaxt="n",yaxt="n",bty="n",xlab="",ylab="",main="Figure 10.1, p438")
axis(1,pos=0)
axis(2,pos=0)
x <- seq(from=0 ,to=5,by=0.1)
cbind(x, bessel_J0(x),bessel_J1(x),bessel_Jn(2,x)) #table 9.1, p390
cbind(x, bessel_Y0(x),bessel_Y1(x),bessel_Yn(2,x)) #table 9.2, p391
t(bessel_Jn_array(3,9,x*2)) #table 9.2, p398
cbind( x, #table 9.8, p416
exp(-x)*bessel_I0 (x),
exp(-x)*bessel_I1 (x),
x^(-2)*bessel_In(2,x)
)
cbind( x, #table 9.8, p417
exp(x)*bessel_K0 (x),
exp(x)*bessel_K1 (x),
x^(2)*bessel_Kn(2,x)
)
cbind(x, #table 10.1 , p457
bessel_j0(x),
bessel_j1(x),
bessel_j2(x),
bessel_y0(x),
bessel_y1(x),
bessel_y2(x)
)
bessel_Jnu(0.5,1.123)
bessel_Ynu(0.5,1.123)
bessel_Inu(0.5,1.123)
bessel_Inu_scaled(0.5,1.123)
bessel_Knu(0.5,1.123)
bessel_lnKnu(0.5,1.123)
bessel_Knu_scaled(0.5,1.123)
x <- 1:20
cbind(x,bessel_zero_J0(x), bessel_zero_J1(x), bessel_zero_Jnu(n=3,x))
#Table 9.5, p409, odd columns
bessel_sequence_Jnu(nu=0.5,v=c(0.4,0.44,0.55))
Run the code above in your browser using DataCamp Workspace