matematicas

# Definir la matriz de coeficientes técnicos
A <- matrix(c(
  0.293, 0.000, 0.000,
  0.014, 0.207, 0.017,
  0.044, 0.010, 0.216
), nrow = 3, byrow = TRUE,
dimnames = list(c("agricultura", "manufactura", "energía"),
                c("agricultura", "manufactura", "energía")))

# a) Unidades de producción agrícola para obtener 1 unidad de producto agrícola
# Esto es simplemente el coeficiente técnico directo
a_a <- A["agricultura", "agricultura"]
cat("a) Unidades de producción agrícola para 1 unidad de producto agrícola:", a_a, "\n")

## a) Unidades de producción agrícola para 1 unidad de producto agrícola: 0.293

# b) Unidades de producción agrícola para 200,000 unidades de producto agrícola
unidades_b <- 200000 * a_a
cat("b) Unidades para 200,000 unidades agrícolas:", unidades_b, "\n")

## b) Unidades para 200,000 unidades agrícolas: 58600

# c) Unidades de producción agrícola para 50,000 unidades de energía
# Necesitamos calcular la demanda directa e indirecta usando la matriz inversa de Leontief
I <- diag(3) # Matriz identidad
L <- solve(I - A) # Matriz inversa de Leontief

# Para producir 50,000 unidades de energía, necesitamos:
demanda <- c(0, 0, 50000) # Vector de demanda final
produccion <- L %*% demanda

unidades_c <- produccion["agricultura"]
cat("c) Unidades agrícolas para 50,000 unidades de energía:", unidades_c, "\n")

## c) Unidades agrícolas para 50,000 unidades de energía: NA

# d) Unidades de energía para 50,000 unidades de productos agrícolas
demanda_d <- c(50000, 0, 0) # Vector de demanda final para agricultura
produccion_d <- L %*% demanda_d

unidades_d <- produccion_d["energía"]
cat("d) Unidades de energía para 50,000 unidades agrícolas:", unidades_d, "\n")

## d) Unidades de energía para 50,000 unidades agrícolas: NA

# Exportaciones
exportaciones <- c(agricultura = 13213, manufactura = 17597, energia = 1786)

# Calcular la producción total necesaria para satisfacer las exportaciones
produccion_export <- L %*% exportaciones

cat("\nProducción total necesaria para las exportaciones:\n")

## 
## Producción total necesaria para las exportaciones:

print(produccion_export)

##                  [,1]
## agricultura 18688.826
## manufactura 22597.858
## energía      3615.162