clc
clear
close all
randn('state',9);
rand('state',9);

noise=0.1;
dim=2;
num=900;
num_pos=500;
num_neg=400;
variance=0.004;

radii = ones(num,1)+0.5+0.05*randn(num_pos+num_neg,1);%+0.22*randn(num_pos+num_neg,1);
phi = rand(num_pos+num_neg,1).*pi;
x=zeros(dim,num);

for i=1:num_pos
    x(2,i)=radii(i)*cos(phi(i))*0.8-.3;
    x(1,i)=radii(i)*sin(phi(i))-.5;
    y(i,1)=1;
end
for i=num_pos+1:num_neg+num_pos
    x(2,i)=1+radii(i)*cos(phi(i));
    x(1,i)=-radii(i)*sin(phi(i))+0.5;
    y(i,1)=2;
end

x = x + sqrt(variance)*randn(dim,num);

dat = [x(1,:)' x(2,:)'];
dat(:,1) = dat(:,1)-min(dat(:,1));
dat(:,1) = dat(:,1)/max(dat(:,1));
dat(:,2) = dat(:,2)-min(dat(:,2));
dat(:,2) = dat(:,2)/max(dat(:,2));
           
U = DifFUZZY(dat, 50);