fork download
  1. N<-25
  2. Q<-10
  3. J<-10
  4. V<-numeric()
  5. for(i in 1:Q){
  6. v<-rep(i/(Q+1),J)
  7. V<-cbind(V,v)
  8. }
  9. P<-rep(1/Q,Q)
  10. lf <- 1.00-0.04*Q
  11. G<-diag(rep(lf,Q))
  12. for(i in 2:(Q-1)){
  13. G[c(i-1,i+1),i]<-(1-lf)/2
  14. }
  15. G[c(1,2),1]<-c(lf/(lf+((1-lf)/2)),((1-lf)/2)/(lf+((1-lf)/2)))
  16. G[c(Q-1,Q),Q]<-c(((1-lf)/2)/(lf+((1-lf)/2)),lf/(lf+((1-lf)/2)))
  17.  
  18. q<-1
  19. tmp<-rep(V[, q], N)
  20. Vq <- matrix(tmp, nrow = N, ncol = J, byrow = TRUE)
  21. #print(tmp)
  22.  
  23.  
  24. c<-matrix(rep(1:9,9),nrow=3, ncol=3, byrow=TRUE)
  25. #print(c[,3])
  26. d<-matrix(rep(c[,3],N),nrow=N,ncol=3, byrow=TRUE)
  27. #print(P)
  28.  
  29. v=0.3
  30. u=0
  31. z=1
  32. print(v^u*(1-v)^(1-u))
  33. print(0.3^0)
  34. print(0.7^1)
  35. #4*9=36
  36.  
Success #stdin #stdout 0.25s 39032KB
stdin
Standard input is empty
stdout
[1] 0.7
[1] 1
[1] 0.7