- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 2 个
- 通用积分
- 0
- 学术水平
- 3 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 931 点
- 帖子
- 76
- 精华
- 0
- 在线时间
- 184 小时
- 注册时间
- 2010-9-15
- 最后登录
- 2022-6-22
硕士生
还不是VIP/贵宾
- 威望
- 0 级
- 论坛币
 - 2 个
- 通用积分
- 0
- 学术水平
- 3 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 931 点
- 帖子
- 76
- 精华
- 0
- 在线时间
- 184 小时
- 注册时间
- 2010-9-15
- 最后登录
- 2022-6-22
|
30论坛币
|
急求R 指令, 画下面三个概率密度函数的图像:1. f(x)=4/3 * (1-x^3) x属于区间(0,1)
2. f(x)=3/7 * x^2 x 属于区间 [1,2]
3. f(x)=1/10 x属于区间 [-2,8]
- > a<-0
- > b<-1
- > myfunc<-function(x,a,b){
- + out<-rep(0,length(x))
- + out[x>a & x<b]<-(4/3)*(1-x[x>a & x<b]^3)
- + return(out)
- + }
- > plot(x, myfunc(x,a,b), type = "l", col="blue", lwd=4)[code]
复制代码 错误于plot(x, myfunc(x, a, b), type = "l", col = "blue", lwd = 4) :
找不到对象'x'
这是老师给的第一题的code,但运行出来有错误,哪位高人指点一下
|
最佳答案
楚韵荆风 查看完整内容
a=0;b=1
x=seq(a,b,by=0.01)
fun1=function(x){ 4/3*(1-x^3)}
plot(x, fun1(x), type = "l", col="blue", lwd=4)
a=1;b=2
x=seq(a,b,by=0.01)
fun2=function(x){ 3/7*x^2}
plot(x, fun2(x), type = "l", col="red", lwd=4)
a=-2;b=8
x=seq(a,b,by=0.01)
fun3=function(x){ n=length(x);re=rep(0,n);for(i in 1:n){ re=0.1} }
plot(x, fun3(x), type = "l", col="green", lwd=4)
回帖推荐
楚韵荆风 发表于7楼
查看完整内容
或者用:
fun1=function(x){ 4/3*(1-x^3)}
fun2=function(x){ 3/7*x^2}
fun3=function(x){ n=length(x);re=rep(0,n);for(i in 1:n){ re=0.1} }
final.fun=function(fun=c("fun1","fun2","fun3"),a,b,eps=0.01)
{ x=seq(a,b,by=eps)
plot(x,fun(x), type = "l", col="red", lwd=4)
}
final.fun(fun1,0,1,eps=0.01)
楚韵荆风 发表于2楼
查看完整内容
a=0;b=1
x=seq(a,b,by=0.01)
fun1=function(x){ 4/3*(1-x^3)}
plot(x, fun1(x), type = "l", col="blue", lwd=4)
a=1;b=2
x=seq(a,b,by=0.01)
fun2=function(x){ 3/7*x^2}
plot(x, fun2(x), type = "l", col="red", lwd=4)
a=-2;b=8
x=seq(a,b,by=0.01)
fun3=function(x){ n=length(x);re=rep(0,n);for(i in 1:n){ re=0.1} }
plot(x, fun3(x), type = "l", col="green", lwd=4)
|