1. អនុគមន៍ថ្លៃដើម ចំណូល និងប្រាក់ចំណេញ

1.1 រកអនុគមន៍ចំណូលR(q)

រូបមន្ត៖ចំណូល=តម្លៃផលិតផល*បរិមាណផលិតផល

R<-makeFun((50-0.5*q)*q~q)
R

## function (q) 
## (50 - 0.5 * q) * q

P<-makeFun((((50-0.5*q)*q)-(100+10*q+0.5*q^2)~q))
P

## function (q) 
## ((50 - 0.5 * q) * q) - (100 + 10 * q + 0.5 * q^2)

## function (q) 
## 100 + 10 * q + 0.5 * q^2

2. រកដេរីវេ C,R,

C<- makeFun((100 + 10 * q + 0.5 * q^2) ~ q)
MC=D(C(q)~q)
MC

## function (q) 
## q + 10

MR= D(R(q)~q)
MR

## function (q) 
## 50 - q

MP= D(P(q)~q)
MP

## function (q) 
## 2 * (20 - q)

សង់ក្រាប MC,MR &MP

MC <- function(q) {
  10 + q
}

MR <- function(q) {
  20 - q
}

# បង្កើតចន្លោះ q
q_vals <- seq(0, 20, by = 0.1)

# គណនាតម្លៃ MC(q) និង MR(q)
MC_vals <- MC(q_vals)
MR_vals <- MR(q_vals)

# គូសក្រាប
P<- plot(q_vals, MC_vals, type = "l", col = "blue", lwd = 2, xlab = "q", ylab = "Cost/Revenue", main = "Graphs of MC and MR")
lines(q_vals, MR_vals, col = "red", lwd = 2)
legend("topright", legend = c("MC(q)", "MR(q)"), col = c("blue", "red"), lwd = 2)

P

## NULL

library(Deriv)
R <- function(q) {
  (100 - 2 * q) * q
}
profit <- function(q) {
  R(q) - C(q)
}
profit_prime <- Deriv(profit)
q_optimal <- uniroot(function(q) profit_prime(q), c(0, 50))$root
q_optimal

## [1] 18