P  = [-3 -2 -2 0 0 0 0 2 2 3 ; 0 1 -1 2 1 -1 -2 1 -1 0];
Tc = [ 1  1  1 2 2 2 2 1 1 1 ];
i=[4 5 6 7];
j=[1 2 3 8 9 10];
plotvec(P(:,i),Tc(:,i),'square');
hold on 
plotvec(P(:,j),Tc(:,j),'+');
title('Input Vectors');
m=['P1 ';'P2 ';'P3 ';'P4 ';'P5 ';'P6 ';'P7 ';'P8 ';'P9 ';'P10'];
x=P(1,:)';
y=P(2,:)';
text(x+0.1,y+0.1,m);
axis([min(P(1,:))-1 max(P(1,:))+1 min(P(2,:))-1 max(P(2,:))+1]);
hold off
T = ind2vec(Tc);
targets = full(T)
net = newlvq(minmax(P),4,[0.6 0.4],0.1);
net.IW{1,1}
net.LW{2,1}
Y = sim(net,P);
Y = vec2ind(Y)