function y = psinorm(K,m,n)
y = max(Psi(K,m,n),Psi(-K,m,n));
return


function y = psinorm0(K)
global Ainl Binl Ain1 Bin1 GG iif m n mylinopt Aing Bing



%y = Psi(K,m,n);
y = psipl(K,m,n);
% should really consider max over +/- K

return