---------------
#  Aufgabe 3  #
---------------

require(linprog)


# LOP Aufgabe 1:

c = c(1,4)
b = c(60,20,40)
v1 = c(5,1,0)
v2 = c(3,2,5)
A = cbind(v1,v2)

c
b
A

res = solveLP(c,b,A,maximum=TRUE)
names(res)

res$solution

# also x1 = x = 4 und x2 = y = 8, das passt.

res$opt               # 36, passt ebenfalls.

res$iter1             # Phase1-Iterationen, hier nicht 
                      # erforderlich -> 0

res$iter2             # 2 Simplex-Iterationen



# LOP Aufgabe 2:

# remove everything from before:
rm(list=ls())

c = c(3,2,-4)
b = c(5,6,2)
v1 = c(1,2,1)
v2 = c(4,4,1)
v3 = c(0,-2,-2)
A = cbind(v1,v2,v3)

c
b
A

res = solveLP(c,b,A,maximum=TRUE)

res$solution

# also x1 = x = 4, x2 = y = 0, x3 = z = 1, das passt.

res$opt                         # 8, passt ebenfalls.

res$iter1             # Phase1-Iterationen, hier nicht 
                      # erforderlich -> 0

res$iter2             # 2 Simplex-Iterationen