N<-25
Q<-10
J<-10
V<-numeric()
for(i in 1:Q){
v<-rep(i/(Q+1),J)
V<-cbind(V,v)
}
P<-rep(1/Q,Q)
lf <- 1.00-0.04*Q
G<-diag(rep(lf,Q))
for(i in 2:(Q-1)){
G[c(i-1,i+1),i]<-(1-lf)/2
}
G[c(1,2),1]<-c(lf/(lf+((1-lf)/2)),((1-lf)/2)/(lf+((1-lf)/2)))
G[c(Q-1,Q),Q]<-c(((1-lf)/2)/(lf+((1-lf)/2)),lf/(lf+((1-lf)/2)))
q<-1
tmp<-rep(V[, q], N)
Vq <- matrix(tmp, nrow = N, ncol = J, byrow = TRUE)
#print(tmp)
c<-matrix(rep(1:9,9),nrow=3, ncol=3, byrow=TRUE)
#print(c[,3])
d<-matrix(rep(c[,3],N),nrow=N,ncol=3, byrow=TRUE)
#print(P)
v=0.3
u=0
z=1
print(v^u*(1-v)^(1-u))
print(0.3^0)
print(0.7^1)
#4*9=36
TjwtMjUKUTwtMTAKSjwtMTAKVjwtbnVtZXJpYygpCmZvcihpIGluIDE6USl7Cgl2PC1yZXAoaS8oUSsxKSxKKQoJVjwtY2JpbmQoVix2KQp9ClA8LXJlcCgxL1EsUSkKbGYgPC0gMS4wMC0wLjA0KlEKRzwtZGlhZyhyZXAobGYsUSkpCmZvcihpIGluIDI6KFEtMSkpewpHW2MoaS0xLGkrMSksaV08LSgxLWxmKS8yCn0KR1tjKDEsMiksMV08LWMobGYvKGxmKygoMS1sZikvMikpLCgoMS1sZikvMikvKGxmKygoMS1sZikvMikpKQpHW2MoUS0xLFEpLFFdPC1jKCgoMS1sZikvMikvKGxmKygoMS1sZikvMikpLGxmLyhsZisoKDEtbGYpLzIpKSkKCnE8LTEKdG1wPC1yZXAoVlssIHFdLCBOKQpWcSA8LSBtYXRyaXgodG1wLCBucm93ID0gTiwgbmNvbCA9IEosIGJ5cm93ID0gVFJVRSkKI3ByaW50KHRtcCkKCgpjPC1tYXRyaXgocmVwKDE6OSw5KSxucm93PTMsIG5jb2w9MywgYnlyb3c9VFJVRSkKI3ByaW50KGNbLDNdKQpkPC1tYXRyaXgocmVwKGNbLDNdLE4pLG5yb3c9TixuY29sPTMsIGJ5cm93PVRSVUUpCiNwcmludChQKQoKdj0wLjMKdT0wCno9MQpwcmludCh2XnUqKDEtdileKDEtdSkpCnByaW50KDAuM14wKQpwcmludCgwLjdeMSkKIzQqOT0zNgo=