sundries

 1using JuMP, Clp
 2
 3model = Model(Clp.Optimizer)
 4
 5@variable(model, 0 <= d <= 100)
 6@variable(model, 0 <= a <= 100)
 7@variable(model, j)
 8@variable(model, k)
 9@variable(model, l)
10
11@constraints(model, begin
12        d + a == 65
13        d == 3 * (a - k)
14        d - k == (a + j) / 2
15        a + j == 3 * (d - l)
16        d - l == 3 * (a - l);
17end)
18JuMP.optimize!(model)
19
20println("d=$(value(d)) a=$(value(a))")