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