function res = psi1(x)
for k=1:length(x)
    if x(k)>=0 res(k)=psi(x(k));
    else res(k) = -psi (-x(k));
    end
end
        