###Ramsey Model,BK's Method
bk< - function(mA,mB,mD,nk);
{
n=nrow(mA);
B=solve(mA)%*%mB;## B:B~
D=solve(mA,mD);##D:D~
B11=B[1:nk,1:nk];

Lamda=diag(eigen(B)Svalues[n:1]);
Lamda1=Lamda[1:nk,1:nk];
Lamda2=Lamda[(nk+1):n,(nk+1):n];

P=eigen(B)$vectors[ ,n:1];
Pi= solve(P);

J=Pi%*%B;
M=Pi%*%D;
P11=Pi[1:nk,1:nk];
P22=Pi[(nk+1):n,(nk+1):n];
P12=Pi[1:nk,(nk+1):n]
P21=Pi[(nk+1):n,1:nk]


F=matrix(NA,n,nk);
G=matrix(NA,n,ncol(D));

F[(nk+1):n,]=-solve(P22)%*%P21;
G[(nk+1):n,]=-solve(P22)%*%solve(Lamda2)%*%M[(nk+1):n,];

F[1:nk,]=B11-B12%*%solve(P22)%*%P21;
G[1:nk,]=D[1:nk,]-B12%*%solve(P22)%*%solve(lamda2)%*%M[(nk+1):n,];
out<-NULL;
out % F=F;
out % G=G;
return(out);
}
alpha=0.7;beta=0.95;sigma=2;
Kbar=(beta*(1-alpha))^(1/alpha);##Kbar=0.1664
Ybar=(Kbar)^(1-alpha)-Kbar;##Cbar=0.4175

mA=matrix(c(Kbar/Ybar,alpha,0,sigma),2,2);
mB=matrix(c(1-alpha,0,-Cbar/Ybar,sigma)2,2);
nk=1;
bk(mA,mB,mD,nk=1);


请问为什么会出现
??? Error using ==> dynare at 114
DYNARE: preprocessing failed
以后在编写程序时如何避免这种问题的出现？？？
谢谢！！