function map=fuzzy(rho,alpha)F1=4;CAT=0;wij = ones(1,F1);wijold=1;load points;ns=size(points);np=ns(1);niter=1while 1>0, for point=1:1:np    a(1)=points(point,1);    a(2)=1-points(point,1);    a(3)=points(point,2);    a(4)=1-points(point,2);        for j=1:CAT+1      T(j) = norm(min(a,wij(j,:)),1)/(alpha+norm(wij(j,:),1));    end    while 1>0,          [Tmax,Jmax]=max(T);        if norm(min(a,wij(Jmax,:)),1) >= rho*norm(a,1)           map(point)=Jmax;           break;        end        T(Jmax)=0;    end    wij(Jmax,:)=min(a,wij(Jmax,:));    if Jmax==CAT+1      CAT=CAT+1;      wij = aug(wij,1);   endendif size(wij)==size(wijold)  if wij == wijold        break;  endendwijold=wij;ss=sprintf('niter=%d CAT=%d',niter,CAT)niter = niter+1;endsave wij wij -asciisave map map -ascii