Rick, /* X*b, where b is a column vector */ should be /* X*b`, where b` is a column vector */ ?