function y=P5plus(x);
% Extracts the positive chiral component
x=x(:);
ndim=max(size(x));
x=reshape(x,ndim/4,4);
%
y(:,1)=zeros(ndim/4,1);
y(:,2)=zeros(ndim/4,1);
y(:,3)=x(:,3);
y(:,4)=x(:,4);
%
x=reshape(x,ndim,1);
y=reshape(y,ndim,1);
% Copyright (C) 2006-2007 Artan Borici.
% This program is a free software licensed under the terms of the GNU General Public License