在我尝试练习 Julia 时,我编写了一个绘制分岔图的程序。我的代码如下:
function bifur(x0,y0,a=1.3,b=0.4,n=1000,m=10000)
i,x,y=1,x0,y0
while i < n && abs(x) < m
x,y = a - x^2 + y, b * x
i += 1
end
if abs(x) < m
return x
else
return 1000
end
end
la = Float64[];
lx = Float64[];
for a=0:400
for j = 1:1000
x0 = rand()
y0 = rand()
x = bifur(x0,y0,a/100)
if x != 1000
push!(la,a/100)
push!(lx,x)
end
end
end
using Gadfly
myplot = Gadfly.plot( x=la, y=lx , Scale.x_discrete, Scale.y_continuous, Geom.point)
draw(PNG("myplot.png",10inch,8inch),myplot)
我得到的输出是这张图片:
为了使我的绘图看起来更像这样: 我需要能够将点大小设置为小至一个像素。然后通过增加迭代长度,我应该能够得到一个更好的分岔图。有谁知道如何在 Julia 的牛虻图中设置点大小?