r = 20;
r0 = 50;
x2r=0.5:2*r+0.5;
xr0=0.5:r0+0.5;
x11r0 = 0.5:1.1*r0+0.5;
x12r0 = 0.5:1.2*r0+0.5;
x13r0 = 0.5:1.3*r0+0.5;
x15r0 = 0.5:1.5*r0+0.5;
//
j2r = sqrt(%pi/(2*r)) * besselj(x2r, r);
j11r0 = sqrt(%pi/(2*r)) * besselj(x11r0, r);
j12r0 = sqrt(%pi/(2*r)) * besselj(x12r0, r);
j13r0 = sqrt(%pi/(2*r)) * besselj(x13r0, r);
fctj_r_10_r0 = j13r0(r:10:r0)
//
y2r = sqrt(%pi/(2*r0)) * bessely(x2r, r0);
y11r0 = sqrt(%pi/(2*r0)) * bessely(x11r0, r0);
y12r0 = sqrt(%pi/(2*r0)) * bessely(x12r0, r0);
y13r0 = sqrt(%pi/(2*r0)) * bessely(x13r0, r0);
y15r0 = sqrt(%pi/(2*r0)) * bessely(x15r0, r0);
fcty_r0_5_1_5r0 = y15r0(r0:5:1.5*r0)
//
//xbasc; plot2d(x2r, j2r), xtitle('Bessel');
//xbasc; plot2d(x11r0', [j11r0',y11r0'], xtitle('Bessel et Hankel'), leg="besselj(r)@bessely(2r)");
//xbasc; plot2d(x11r0, y11r0),xtitle('Hankel');
//xbasc; plot2d(x12r0, y12r0),xtitle('Hankel');
//xbasc; plot2d(x15r0, y15r0),xtitle('Hankel');
//
xbasc; plot2d(x12r0', [abs(j12r0'), abs(j12r0.*y12r0)'], logflag="nl", leg="besselj(r)@besselj(r)*hankel(r0)") , xtitle('Hankel*Bessel');
//