第一题可不可以这么做啊？

A = rcauchy(500)

B = cumsum(A)/(c(1:500))

plot(A~B)

恩，不错。不过我觉得用plot(B,A)可能符合题意

发第二题答案：

1)  n = 100000, m = 10000, k = n/m, pi = 0.05, eta = 0.97, theta = 0.95

> n = 100000

> m = 10000

> k = n/m

> Pi = 0.05

> Eta = 0.97

> Theta = 0.95

> p = 1 - ( Eta*Pi + (1-Theta)*(1-Pi) )

> N = 0

>

> for (i in 1:m){

+ if(rbinom(1,1,p^k) == 1){

+ N = N + 1

+ }

+ else{

+ N = N + k + 1

+ }

+ }

>

> ev.N = n*(1 + (1/k) - (p^k))

> ev.N

2)  n = 100000, m =   5000, k = n/m, pi = 0.05, eta = 0.97, theta = 0.95

> n = 100000

> m = 5000

> k = n/m

> Pi = 0.05

> Eta = 0.97

> Theta = 0.95

> p = 1 - ( Eta*Pi + (1-Theta)*(1-Pi) )

> N = 0

>

> for (i in 1:m){

+ if(rbinom(1,1,p^k) == 1){

+ N = N + 1

+ }

+ else{

+ N = N + k + 1

+ }

+ }

>

> ev.N = n*(1 + (1/k) - (p^k))

> ev.N