def projection_matrix(w):
return np.identity(<NUM_LIT:3>) - np.dot(np.reshape(w, (<NUM_LIT:3>,<NUM_LIT:1>)), np.reshape(w, (<NUM_LIT:1>, <NUM_LIT:3>)))<EOL>
Return the projection matrix  of a direction w.