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 _______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion