clear all
 s=20;
a=13;
b=1;
m=32;
for I=1:m
       n1=a*s+b;
u1=mod(n1,m);
v1=u1/m;
D1=fix(10*v1)+1;
  if D1==1 
           PD1=0.11;
       elseif D1==2 
           PD1=0.12;
       elseif D1==3 
           PD1=0.09;
       elseif D1==4
           PD1=0.08;
       elseif D1==5 
           PD1=0.12;
       elseif D1==6
           PD1=0.10;
       elseif D1==7
           PD1=0.09;
       elseif D1==8
           PD1=0.09;
       elseif D1==9
           PD1=0.10;
       elseif D1==10
           PD1=0.10;
       end
       Q=PD1/0.12; %revisar como se calculaba Q
s=u1;
n2=a*s+b;
u2=mod(n1,m);
v2=u2/m;
if v2<Q
    D2= D1;
    s=u2;
          D2
      else
      I==I+1;
        end
    end
