Debugging

You may check your formula by running ``sch.m'' and then typing W'*O(W), which should now be the $4 \times 4$ identity matrix (to within machine precision).