P = [1 2; 2 2; 1 1]';
p1 = [2; 1.5];
Tc = [1 2 3];
T = ind2vec(Tc);
spread = 1;
net = newpnn(P,T,spread);

a = sim(net,p1);
ac = vec2ind(a);
hold on
plot(p1(1),p1(2),'.','markersize',30,'color',[1 0 0])
text(p1(1)+0.1,p1(2),sprintf('class %g',ac))
hold off
title('Classifying a new vector.')
xlabel('P(1,:) and p(1)')
ylabel('P(2,:) and p(2)')
pause;