# Pkg.add("Decimals") # Not needed because this is on juliabox

using Decimals

# Defining Variables
a = Decimal(0.15)
b = Decimal(0.2)
c = Decimal(1)
d = Decimal(-2)
n = Decimal(-0.5)
t = 0
phi = Decimal(0.5)
values = []
column = Dict()

# Creating the first column
while n < Decimal(0.5)
    column[string(n)] = a + b*n^1 + c*n^2 + d*n^3
    n = n + Decimal(0.02)
end

# Add the first column to the array of values
push!(values, column)

# Starting an infinite loop
while true
    n = Decimal(-0.5)
    column = Dict()

    # Creating the n'th column
    while n < Decimal(0.5)
        index = string(n)
        previousIndex = string(n-Decimal(0.02))
        nextIndex = string(n+Decimal(0.02))
        if index == "0.5"
            column[index] = Decimal(0.25)
        else
            column[index] = values[end][index] + (values[end][previousIndex]*phi - 2*values[end][n]*phi + values[end][nextIndex]*phi)
        end
        n = n+Decimal(0.02)
    end

    if values[end]["0.48"] == values[end]["0.46"]
      break
    end

    push!(values, column)
    t = t + 4
end

println("done!")