[Numpy-discussion] matrix default to column vector?

[Alan G Isaac]

On 6/4/2009 12:08 PM Olivier Verdier apparently wrote:
> I really don't see any advantage of matrices over arrays for teaching. I 
> prefer to teach linear algebra with arrays. 


beta = (X.T*X).I * X.T * Y

beta = np.dot(np.dot(la.inv(np.dot(X.T,X)),X.T),Y)

I rest my case.

I would have to switch back to GAUSS for teaching
if NumPy discarded matrices.

Cheers,
Alan Isaac