from Proj.m (Projection for matrix) by rajchandar k
Projection file .

out=proj(v)
function out=proj(v)
out=(1/dot(v,v))*v*v';

Contact us at files@mathworks.com