clear all
load data
procent=round(size(data, 1) * 0.85);
Xu=data(1:procent,1:2)';
Du=data(1:procent, 3)';
Xw=data((procent+1):end, 1:2)';
Dw=data((procent+1):end, 3)';

B1=(Du==1);
B2=(Du==2);
B3=(Du==3);
Xu1=Xu(:,B1);
Xu2=Xu(:,B2);
Xu3=Xu(:,B3);
plot(Xu1(1,:), Xu1(2,:),".r", Xu2(1,:), Xu2(2,:),".g", Xu3(1,:), Xu3(2,:),".b")