Ejercicio 3.

Nombre: Danna Valentina Angarita Quintero Código: 1151984

Unidades de materia prima/unidad

de solución
Disponibilidad
Solución A Solución B
diaria (unidades)
MI 0.5 0.5 150
MII 0.6 0.4 145
Utilidad por tonelada
8 10
(miles)

Para formular el modelo matemático se definen las variables de decisión:

x 1 : unidades diarias fabricadas de solución A

x 2 : unidades diarias fabricadas de solución B

La función objetivo es

Maximizar U =8 x1 +10 x 2 (en miles)

Sujeta a las restricciones

0.5 x 1+ 0.5 x 2 ≤150 (1)

0.6 x 1+ 0.4 x 2 ≤ 145 (2)

30 ≤ x 1 ≤150

x 1 ≥ 30 (3)

x 1 ≤ 15 0 (4)

x 2 ≥ 40 (5)

x 2 ≤ 20 0 (6)

x 1 , x 2 ≥ 0 (Condición de no negatividad)
Solución con R

a)
> coef<-c(8,10)

> A<-matrix(c(0.5,0.6,1,1,0,0,0.5,0.4,0,0,1,1),ncol=2)

> b<-c(150,145,30,150,40,200)

> dir<-c("<=","<=",">=","<=",">=","<=")
> solucion <- solveLP(coef, b, A, maximum=TRUE, dir)

> summary(solucion)

Results of Linear Programming / Linear Optimization

Objective function (Maximum): 2800

Solution

opt

1 100

2 200

b)

> library(lpSolve)

> coef<-c(8,10)

> A<-matrix(c(0.5,0.6,1,1,0,0,0.5,0.4,0,0,1,1),ncol=2)

> b<-c(150,145,30,150,40,200)

> dir <- c("<=","<=",">=","<=",">=","<=")

> solucion <- lp('max', coef, A, dir, b)

> solucion$objval

[1] 2800

> solucion$solution

[1] 100 200

> solucion <- lp('max', coef, A, dir, b,compute.sens=TRUE)

> solucion$duals

[1] 16 0 0 0 0 2 0 0

Restricción Dual
1 16
2 0
3 0
4 0
5 0
6 2

 Un aumento de 1 unidad en la restricción 1 ocasiona un incremento de la función

objetivo en 16.

 Un aumento de 1 unidad en la restricción 2 no provoca alteración alguna sobre la

función objetivo

 Un aumento de 1 unidad en la restricción 3 no provoca alteración alguna sobre la

función objetivo
  Un aumento de 1 unidad en la restricción 4 no provoca alteración alguna sobre la
función objetivo

 Un aumento de 1 unidad en la restricción 5 no provoca alteración alguna sobre la

función objetivo

 Un aumento de 1 unidad en la restricción 6 ocasiona un incremento de la función

objetivo en 2.

Hallamos los intervalos de factibilidad para cada una de las restricciones:

> solucion$duals.from

[1] 1.15e+02 -1.00e+30 -1.00e+30 -1.00e+30 -1.00e+30 1.75e+02 -1.00e+30

[8] -1.00e+30

> solucion$duals.to

[1] 1.541667e+02 1.000000e+30 1.000000e+30 1.000000e+30 1.000000e+30

[6] 2.700000e+02 1.000000e+30 1.000000e+30

Restricción Dual Actual Mínima Máxima

1 16 150 115 154
2 0 145 -Inf Inf
3 0 30 -Inf Inf
4 0 150 -Inf Inf
5 0 40 -Inf Inf
6 2 200 175 270

Hallamos los intervalos de factibilidad para los coeficientes de la función objetivo:

> solucion$sens.coef.from

[1] 0 8

> solucion$sens.coef.to

[1] 1e+01 1e+30

Variable Actual Mínima Máxima

x 1(A) 8 0 10
x 2 (B) 10 8 Inf