r0 = 0.89999
n = 2.7
vc = [0, 0, 0]
function create(n, r0, vc)
    state = Dict{}()
    state["coord"] = []; state["vel"] = []; state["force"] = []; state["pe"] = []; state["size"] = []; state["base"] = []
    state["size"] = [5*n, 5*n, 5*n]
    vc = [0, 0, 0]
    for i in range(int(floor(n/r0)))
        for j in range(int(floor(n/r0)))
            for k in range(int(floor(n/r0)))
                push!(state["coord"],[r0 * (i + 0.5), r0 * (j + 0.5), r0 * (k + 0.5)])
                push!(state["base"],[r0 * (i + 0.5), r0 * (j + 0.5), r0 * (k + 0.5)])
                push!(state["vel"],[randint(-2, 2) / 2., randint(-2, 2) / 2., randint(-2, 2) / 2.])
                push!(state["force"],[0.0, 0.0, 0.0])
                push!(state["pe"], 0.0)
            end
        end
    end

    for i in range(int(floor(n/r0))^3)
        for idx in range(3)
            vc[idx] += state["vel"][i][idx]
        end
    end
    for i in range(3)
        vc[i] /= int(floor(n/r0))
    end
    for i in range(int(floor(n/r0))^3)
        for idx in range(3)
            state["vel"][i][idx] -= vc[idx]
        end
    end
    update_forces(state)
    return state
end

system = create(n, r0, vc)