Hola,
les mando las instrucciones para graficar las funciones de Bessel de primer orden en Octave-Matlab y en Mathematica
Octave-Matlab
Para graficar en el intervalo [0,100], creamos la variable x en ese intervalo, variando cada 0.001
x=0:0.001:100;
Generamos los primeros 3 (n=0,1,2) y los graficamos
y0=besselj(0,x);
y1=besselj(1,x);
y2=besselj(2,x);
Mathematica
Para graficar las primeras tres funciones de Bessel, usamos la siguiente instrucción:
Plot[{BesselJ[0,x],BesselJ[1,x],BesselJ[2,x]}, {x,0,100}]
Para Mathematica 6.0 o versiones superiores pueden utilizar el comando Manipulate para explorar como se ven cuando variamos el índice
Manipulate[Plot[BesselJ[n,x], {x,0,100}],{n,1,25,1}]
Con ello aparece un boton que pueden desplazar y observaran como varian las funciones de Bessel cuando n crece en pasos de uno de 1 a 25.
En el siguiente post les indicaré como pueden graficarlo utilizando la serie de potencia que define a las funciones. Con la información de los post anteriores es posible que lo puedan hacer.