非常感谢帮忙的人!m<(_ _)>m
程序如下:
u=runif(1000)
v=runif(1000)
if v<=u y==0
else y==1
s=sum(y*exp(-((1/u)+(1/v)-2))*(1/(u*v)^2))/1000
s
显示结果如下:(就是错误的地方不知道如何改……)
> u=runif(1000)
> v=runif(1000)
>
> if v<=u y==0
错误: 意外的符号在"if v"里
>
> else y==1
错误: 意外的'else'在"else"里
> s=sum(y*exp(-((1/u)+(1/v)-2))*(1/(u*v)^2))/1000
> s
[1] 0.503547


雷达卡




京公网安备 11010802022788号







