clear
T=100000;
eps=0.9;
x0=10;
G=rand(1,T,"normal");
//G=grand(1,T,'nor',5,1); 
d=1000;
//G=floor(d*rand(1,T,'def'));
U=grand(1,T,"def");
X=[];
X=[X,x0];
for s=1:T



             if U(s)>eps  then X=[X,G(s)];
                               else X=[X, X(floor(1+s*grand(1,1,"def")))]; 
             end

end

xset("window",1);
     xbasc();
plot2d(1:T+1,X)

xset("window",2);
     xbasc();
histplot(100,X')