8.2 Programs
145
%eta=zeros(N,NR);
%for j=1:NR
%for i=1:N-1
%Pat1at = 0.5*(1. + exp(-i*dt/tau)); % Probability to stay in "-1"
%Pat2at1 = Pat1at;
%Pat1bt = 0.5*(1. - exp(-i*dt/tau)); %Transition Probability 1->0
%Pat2bt1 = Pat1bt;
%Pbt2bt1 = 0.5*(1. + exp(-i*dt/tau)); %Probability to stay in "1"
%Pbt1bt = Pbt2bt1;
%Pbt1at = 0.5*(1. - exp(-i*dt/tau)); %Transition Probability 0->1
%Pbt2at1 = Pbt1at;
% R=unifrnd(0,1);
% if(Pat1at > R)
% eta(i+1,j)=-1.;
% else
% eta(i+1,j)=1.;
% end
%if(Pat1at < R || Pat1at > R)
% R1=unifrnd(0,1);
% elseif(Pat2bt1 > R1 || Pat2at1 > R1) % || significa or
% eta(i+2,j)=-1.;
%else
% eta(i+2,j)=1.;
%end
%end
%end
%=========================================================
DU=zeros(N-1,NR);
Dd1U=zeros(N-1,NR);
Dd2U=zeros(N-1,NR);
Dd3U=zeros(N-1,NR);
DDeff=zeros(N-1,NR);
%Potential slightly tilded to the right.? C1=0. Denisov et al potential
DU(1,1)=(cos(kL*xij) - C1*0.25*cos(2*kL*xij));
Dd1U(1,1)=-sin(xij) + C1* 0.5*sin(2.*xij);
Dd2U(1,1)=-cos(kL*xij) + C1*cos(2.*kL*xij);
Dd3U(1,1)=sin(kL*xij) - C1*2.*sin(2.*kL*xij);
DDeff(1,1)=1./(1-Dlambda*Dbeta*Dd2U(1,1));
%=========================================================
% THERMAL QUANTUM NOISE
%=========================================================
tR=tc/ts;
INVTR=1/tR;
A1=6.50618; A2=3.89653; A3=3.01154;
t1=8.33507e-5; D1=A1*t1;