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

using Decimals

# Defining Variables
Da = Decimal(3e-5)
Db = Decimal(4e-5)
dx = Decimal(0.003)
dt = Decimal(1)
n = Decimal(-0.5)
t = 0

phia = Decimal((Da*dt)/(dx*dx))
phib = Decimal((Db*dt)/(dx*dx))

cs1 = decimal(0)
cs2 = decimal(1)

values = []
column = Dict()

# Creating the first column
# while n <= Decimal(0)
#     if n == Decimal(0)
#         column[string(n)] = .36
#     else
#         column[string(n)] = 0
#         n = n + Decimal(0.02)
#     end
# end
 println("a")
while n >= Decimal(0)
    if n == Decimal(0)
        column[string(n)] = .96
    else
        column[string(n)] = 1
        n = n + Decimal(0.02)
    end
end

println("a")

# 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(Decimal(n-Decimal(0.02)))
        if previousIndex == "-0.00"
            previousIndex = "0"
        end
        nextIndex = string((n+Decimal(0.02)))
        if index == "-0.5" || index == "0.5"
            column[index] = a + b*n^1 + c*n^2 + d*n^3
        else
            column[index] = Decimal(values[end][index] + (values[end][previousIndex]*phi - 2*values[end][index]*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(length(values))
println(values[1]["-0.48"])
println("done!")