In [1]:

    

using PyPlot


    

In [2]:

    

function generate_data(n)
    epsilon_values=Array(Float64,n)
    for i=1:n
        epsilon_values[i] = randn()
    end
    return epsilon_values
end

ts_length = 100
data = generate_data(ts_length)
plot(data, "b-")


    

    













    
Out[2]:





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

In [3]:

    

ts_length = 100
data = randn(ts_length)
plot(data, "b-")


    

    













    
Out[3]:





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

In [4]:

    

Pkg.add("Distributions")


    

    




INFO: No packages to install, update or remove
INFO: Package database updated
INFO: METADATA is out-of-date — you may not have the latest version of Distributions
INFO: Use `Pkg.update()` to get the latest versions of your packages

In [5]:

    

Pkg.update()


    

    




INFO: Updating METADATA...
INFO: Computing changes...
INFO: No packages to install, update or remove

In [8]:

    

using PyPlot
using Distributions

function plot_histogram(distribution, n)
    epsilon_values = rand(distribution, n)  # n draws from distribution
    plt[:hist](epsilon_values)
end

lp = Laplace()
plot_histogram(lp, 500)


    

    













    
Out[8]:





([10.0,33.0,81.0,177.0,127.0,44.0,18.0,7.0,1.0,2.0],[-3.7234031778113614,-2.7363658537428535,-1.749328529674346,-0.7622912056058384,0.22474611846266956,1.2117834425311775,2.1988207665996846,3.1858580906681926,4.172895414736701,5.1599327388052085,6.1469700628737165],Any[PyObject <matplotlib.patches.Rectangle object at 0x000000001B685588>,PyObject <matplotlib.patches.Rectangle object at 0x000000001B9A1198>,PyObject <matplotlib.patches.Rectangle object at 0x000000001B998EF0>,PyObject <matplotlib.patches.Rectangle object at 0x000000001B99DBE0>,PyObject <matplotlib.patches.Rectangle object at 0x000000001B9A6D68>,PyObject <matplotlib.patches.Rectangle object at 0x000000001B9AA860>,PyObject <matplotlib.patches.Rectangle object at 0x000000001B9AAFD0>,PyObject <matplotlib.patches.Rectangle object at 0x000000001B9B2E10>,PyObject <matplotlib.patches.Rectangle object at 0x000000001B998240>,PyObject <matplotlib.patches.Rectangle object at 0x000000001B9B7F98>])

In [7]:

    

rand(3)


    

    
Out[7]:





3-element Array{Float64,1}:
 0.916872
 0.396554
 0.871333

In [9]:

    

factorial(6)


    

    
Out[9]:





720

In [11]:

    

function factorial2(n)
    k=1
    for i in 1:n
        k=k*i
    end
    return k
end
factorial2(6)


    

    
Out[11]:





720

In [14]:

    

function binomial_rv(n, p)
    count = 0
    U = rand(n)
    for i in 1:n
        if U[i] < p
            count = count + 1    # Or count += 1
        end
    end
    return count
end

for j in 1:20
    b = binomial_rv(15, 0.8)
    print("$b, ")
end


    

    




12, 12, 13, 10, 12, 14, 13, 13, 12, 13, 14, 12, 13, 10, 9, 13, 10, 13, 14, 9, 

In [16]:

    

n=100000
count=0
for i in 1:n
    u,v=rand(2)
    d=sqrt((u-0.5)^2+(v-0.5)^2)
    if d<0.5
        count+=1
    end
end
area_estimate=count/n
print(area_estimate*4)


    

    




3.14024

In [ ]: