A = createA(n)
    b = createB(n)
    c = createC(n)
    m = Model(solver = s)

    @variable(m, x[1:n] >= 0)

    @objective(m, Min, sum(c[i]*x[i] for i=1:n))

    @constraint(m, A*x .>= b)

    print(m)
    println("Objective value: ", getobjectivevalue(m))
    println("x = ", getvalue(x))

    return solve(m)