រូបមន្ត៖ចំណូល=តម្លៃផលិតផល*បរិមាណផលិតផល
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 <- 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
#3.1 រកបរិមាណដែលធ្វើអោយប្រាក់ចំណេញអតិបរមាq*
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
#3.2 គណនាប្រាក់ចំណេញអតិបរមា
R <- function(q) {(100 - 2 * q) * q}
C <- function(q) {50 + 20 * q + q^2}
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] 13.33333