Universidade Federal da Paraíba

Curso de Graduação em Economia

Análise Microeconômica I

options(OutDec=",")

# capital fixo - hectares de terra
x1 <- 10
# Criar sequências de 10 números ordenados - trabalhadores
x2 <- seq(1,10,length=10)

# Definir função de produção
y <- function(K,L) {K - L + 6*L^2 -0.5*L^3}

# Definir função de produtividade marginal:
pmg <- function(K,L) {-1 + 12*L -1.5*L^2}

# Definir função de produtividade média:
pme <- function(K,L) {(K - L + 6*L^2 -0.5*L^3)/L}


# desenhar gráficos:

par(oma=c(0,1.8,0,0.8))
plot(x2,pmg(x1,x2),type="o",bty="l",cex.lab=0.8,cex.axis=0.8,col="brown",lwd=2,las=1,
     ylab=expression(f(K,L)),xlab=expression(L))
grid(col=grey(0.80))
points(x2,pme(x1,x2),type="o",col="blue")
par(fig=c(0, 1, 0, 1), oma=c(0, 0, 0, 0),mar=c(0, 0, 0, 0), new=TRUE)
plot(0,axes=F,xlab="",ylab="",col="transparent")
legend("topright",c(expression(PMgL(K,L) == -1 + 12*L -1.5*L^2),expression(PMeL(K,L) == (K - L + 6*L^2 -0.5*L^3)/L)),
       bty="n", pch =c(1,1),col=c("brown","blue"),cex=0.9)

# resultados da simulação:
data.frame(K=x1,L=round(x2,1),Y=round(y(x1,x2),1),PMgL=round(pmg(x1,x2),1),PMeL=round(pme(x1,x2),1))
##     K  L     Y  PMgL PMeL
## 1  10  1  14,5   9,5 14,5
## 2  10  2  28,0  17,0 14,0
## 3  10  3  47,5  21,5 15,8
## 4  10  4  70,0  23,0 17,5
## 5  10  5  92,5  21,5 18,5
## 6  10  6 112,0  17,0 18,7
## 7  10  7 125,5   9,5 17,9
## 8  10  8 130,0  -1,0 16,2
## 9  10  9 122,5 -14,5 13,6
## 10 10 10 100,0 -31,0 10,0