Continuamos explorando las posibilidades de hacer plots en julia. Hoy el paquete Gaston, que implementa un interfaz a gnuplot, que tiene que ser instalado en el sistema separadamente.

julia>   Pkg.add("Gaston")

julia> using Gaston

julia> gaston_demo()
WARNING: contains(collection, item) is deprecated, use in(item, collection) instead
 in contains at reduce.jl:238

Se abre rápidamente muchos ventanas de gnuplot con plots.  Mostramos unos ejemplos, tomados del manual de Gaston.

julia> t = 0:0.0001:.15
0.0:9.999999999999999e-5:0.15

julia> length(t)
1501

julia> carrier = cos(2pi*t*200) ;

julia> modulator = 0.7+0.5*cos(2pi*t*15) ;

julia> am = carrier.*modulator ;

julia> plot(t,am,"color","blue","legend","AM DSB-SC","linewidth",1.5,
       t,modulator,"color","black","legend","Envelope",
       t,-modulator,"color","black","title","AM DSB-SC example",
       "xlabel","Time (s)","ylabel","Amplitude",
       "box","horizontal top left")
0

julia> set_filename("ejemplo_plot_2.png")
"ejemplo_plot_2.png"

julia> printfigure("png")
1

Y ahora el plot esta escrito a fichero, y aparece como abajo:

Continuamos con otro ejemplo tomado del manual, un histograma de datos simulados de une densidad jicuadrada con dos grados de libertad (el man ual lo llama "densidad  Rayleigh"):

julia> y = sqrt( randn(10000).^2 + randn(10000).^2 )  ;

julia> histogram(y,"bins",25,"norm",1,"color","blue","title","Rayleigh density (25 bins)")

1

julia> set_filename("ejemplo_plot_3.png")

"ejemplo_plot_3.png"

julia> printfigure("png")

1