jueves, 22 de abril de 2010

Funciones de Bessel (Graficar)




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);
plot(x,y,x,y1,x,y2)



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.

No hay comentarios:

Publicar un comentario