function Y=rpsi(a,b)
% rpsi - Version of PSI that returns the same shape array as supplied
%
% rpsi :: [[Size]->real] -> [[Size]->real].
if nargin==1
	Y=reshape(psi(a),size(a));
else
	Y=reshape(psi(a,b),size(b));
end

