% discretization

function [u,ul1,ul2,umin,umax]=disc(uz)

global T my sn

 umax=max(uz);
 umin=min(uz);
 ur=umax-umin;

% shifting times(sn) and shifting internals(siy,siu)
 
  sn=6;                                   % shift numbers sn=4
  siu=ur/sn;                              % shift interval of y
  
for s=1:sn-2                                % shifting
    for i=1:T
      ul1(s)=umin+siu*s;
      ul2(s)=umin+siu*(s+1);
      if uz(i)<=ul1(s)
         u(s,i)=1;
      elseif uz(i)<=ul2(s)
         u(s,i)=2;
      else
         u(s,i)=3;
      end
    end
end