R 向量上的循环单根
下面的代码正在求解q1=10的方程,最终得到根q2=170R 向量上的循环单根,r,loops,vector,R,Loops,Vector,下面的代码正在求解q1=10的方程,最终得到根q2=170 q1 <- 10 fun <- function(q2) 100-q1-0.5*q2-5 root <- uniroot(fun, c(0, 200))$root q1uniroot提供了一个..参数,可用于向函数传递其他命名参数。因此,您可以将q1作为参数添加到fun中,并将其作为命名参数传递给sapply中的uniroot: fun <- function(q1, q2) 100-q1-0.5*q2-5 s
q1 <- 10
fun <- function(q2) 100-q1-0.5*q2-5
root <- uniroot(fun, c(0, 200))$root
q1uniroot
提供了一个..
参数,可用于向函数传递其他命名参数。因此,您可以将q1
作为参数添加到fun
中,并将其作为命名参数传递给sapply
中的uniroot
:
fun <- function(q1, q2) 100-q1-0.5*q2-5
sapply(seq(1, 200, 0.05), function(q1) uniroot(fun, c(-1000, 200), q1=q1)$root)
# [1] 188.0 187.9 187.8 187.7 187.6 187.5 187.4 187.3 187.2 187.1 187.0 186.9 186.8 186.7
# [15] 186.6 186.5 186.4 186.3 186.2 186.1 186.0 185.9 185.8 185.7 185.6 185.5 185.4 185.3
# ...
有趣
fun <- function(q1, q2) 100-q1-0.5*q2-5
sapply(seq(1, 200, 0.05), function(q1) uniroot(fun, c(-1000, 200), q1=q1)$root)
# [1] 188.0 187.9 187.8 187.7 187.6 187.5 187.4 187.3 187.2 187.1 187.0 186.9 186.8 186.7
# [15] 186.6 186.5 186.4 186.3 186.2 186.1 186.0 185.9 185.8 185.7 185.6 185.5 185.4 185.3
# ...