EjrALGpseudoinv01

(%i2)

A : matrix([ 1, 1],[ 1, 0],[ - 1, 0]) ;
AtA : transpose( A). A ;

\begin{matrix} (A) & (\begin{matrix} 1 & 1 \\ 1 & 0 \\ - 1 & 0 \end{matrix}) \end{matrix}

\begin{matrix} (AtA) & (\begin{matrix} 3 & 1 \\ 1 & 1 \end{matrix}) \end{matrix}

(%i3)

X : addcol( AtA, ident( 2)) ;

\begin{matrix} (X) & (\begin{matrix} 3 & 1 & 1 & 0 \\ 1 & 1 & 0 & 1 \end{matrix}) \end{matrix}

(%i4)

X : rowop( X, 1, 2, 2) ;

\begin{matrix} 0 errores, 0 advertencias \end{matrix}

\begin{matrix} (X) & (\begin{matrix} 1 & - 1 & 1 & - 2 \\ 1 & 1 & 0 & 1 \end{matrix}) \end{matrix}

(%i5)

X : rowop( X, 2, 1, 1) ;

\begin{matrix} (X) & (\begin{matrix} 1 & - 1 & 1 & - 2 \\ 0 & 2 & - 1 & 3 \end{matrix}) \end{matrix}

(%i6)

X : rowop( X, 2, 2, 1 / 2) ;

\begin{matrix} (X) & (\begin{matrix} 1 & - 1 & 1 & - 2 \\ 0 & 1 & - \frac{1}{2} & \frac{3}{2} \end{matrix}) \end{matrix}

(%i7)

X : rowop( X, 1, 2, - 1) ;

\begin{matrix} (X) & (\begin{matrix} 1 & 0 & \frac{1}{2} & - \frac{1}{2} \\ 0 & 1 & - \frac{1}{2} & \frac{3}{2} \end{matrix}) \end{matrix}

(%i8)

invA : submatrix( X, 1, 2) ;

\begin{matrix} (invA) & (\begin{matrix} \frac{1}{2} & - \frac{1}{2} \\ - \frac{1}{2} & \frac{3}{2} \end{matrix}) \end{matrix}

(%i9)

L : invA. transpose( A) ;

\begin{matrix} (L) & (\begin{matrix} 0 & \frac{1}{2} & - \frac{1}{2} \\ 1 & - \frac{1}{2} & \frac{1}{2} \end{matrix}) \end{matrix}

(%i10)

L. A ;

\begin{matrix} (%o10) & (\begin{matrix} 1 & 0 \\ 0 & 1 \end{matrix}) \end{matrix}

Created with wxMaxima.