In [1]:

    

include("kicked-rotor.jl")


    

    




Syntax: krawczyk2d(function, [Interval(lo, hi), Interval(lo, hi)], precision [default is 64])






    
Out[1]:





zeros4 (generic function with 1 method)

In [2]:

    

;pwd


    

    




/Users/dsanders/Dropbox/postdocs/nikolay_kryukov/interval-methods/draft

In [3]:

    

;cd ..


    

    




/Users/dsanders/Dropbox/postdocs/nikolay_kryukov/interval-methods

In [ ]:

    

include("kicked-rotor.jl")


    

In [10]:

    

import Base.zeros


    

    




Warning: import of Base.zeros into Main conflicts with an existing identifier; ignored.

In [7]:

    

using PyPlot


    

In [2]:

    

function plot_points(N, k=0.5)
    data = Base.zeros(N, 2)
    x0 = [2pi*rand(), 2pi*rand()]

    for i in 1:N
        output = T0(x0, k, i)
        data[i, 1:2] = output
    end
    
    data
end


    

    
Out[2]:





plot_points (generic function with 2 methods)

In [22]:

    

k = 0.5
for i in 1:500
    points = plot_points(200, k);
    plt.plot(points[:,1], points[:,2], "b.", markersize=1, alpha=0.3)
end

plt.ylim(2,4)
plt.xlabel(L"$\theta$")
plt.ylabel(L"$\phi$")

plt.savefig("standard_k$k.pdf")


    

In [4]:

    

using PyPlot


    

    




INFO: Loading help data...
Warning: using PyPlot.plot in module Main conflicts with an existing identifier.

In [35]:

    




    

    













    
Out[35]:





1-element Array{Any,1}:
 PyObject <matplotlib.lines.Line2D object at 0x11f4e05d0>

In [ ]: