from scipy.optimize import linprog
obj = [-2, -8, -4]

lhs_ineq = [[ 2, 1, 1],  
            [1,3, 0], 
            [ -1, 2, -3]]  

rhs_ineq = [20,  
            30,  
             15]  


bnd = [(0, float("inf")),  
       (0, float("inf"))]  
opt=linprog(c=obj,A_ub=lhs_ineq,b_ub=rhs_ineq, bounds=bnd,method="highs")
print(opt)