- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 4770 个
- 通用积分
- 3737.0293
- 学术水平
- 167 点
- 热心指数
- 195 点
- 信用等级
- 163 点
- 经验
- 54405 点
- 帖子
- 593
- 精华
- 0
- 在线时间
- 721 小时
- 注册时间
- 2008-4-15
- 最后登录
- 2026-1-5
|
表示实根?这个应该是初一的数学吧?
- get_root <- function(a,b,c){
- delta <- b^2-4*a*c
- if(delta < 0)
- return(print('no real root'))
-
- x1 <- (-b + sqrt(delta))/(2 * a)
- x2 <- (-b - sqrt(delta))/(2 * a)
-
- return(c(x1, x2))
- }
- get_prob <- function(n_sim){
- a <- runif(n_sim, 1, 5)
- b <- rnorm(n_sim, 3, 10)
- c <- rexp(n_sim, rate = 1) # mean is 1/rate
- return(sum(b^2-4*a*c >= 0)/n_sim)
- }
复制代码
|
|