function QM_PIBplot(Psi,x,A,L,parity)
% QM_PIBplot(Psi,x,A);

if parity
    signPsi = 1;
    PsiColor = 'r-';
else
    signPsi = -1;
    PsiColor = 'g-';
end

top = 1.25*A;
x = x(:);
x = [flipud(-x(2:end)) ;x];
Psi = Psi(:);
Psi = [flipud(signPsi*Psi(3:end)) ;Psi(2:end)];


clf
set(gca,'nextplot','add')
plot([-L L],[0 0],'k-',[-L -L],[0 top],'k-',[L L],[0 top],'k-',...
    [-1.5*L -L],[top top],'k-',[1.5*L L],[top top],'k-','linewidth',2)
axis equal
set(gca,'ylim',[-top top*1.1])

plot(x,Psi,PsiColor,'linewidth',4)


return