{VERSION 5 0 "Linux" "5.0" } {USTYLETAB {PSTYLE "Ordered List 1" -1 200 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 1 } {PSTYLE "Ordered List 2" -1 201 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 36 2 0 2 2 -1 1 }{PSTYLE "Ordered List 3" -1 202 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 72 2 0 2 2 -1 1 }{PSTYLE "Ordered List 4" -1 203 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 108 2 0 2 2 -1 1 }{PSTYLE "Ordered List 5" -1 204 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 144 2 0 2 2 -1 1 }{PSTYLE "Author" -1 19 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 8 8 2 0 2 0 2 2 -1 1 }{PSTYLE "Warning" -1 7 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 1 0 0 1 } 1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Annotation Title" -1 205 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 1 0 0 1 }3 1 0 0 12 12 2 0 2 0 2 2 -1 1 }{PSTYLE "Fixed Width" -1 17 1 {CSTYLE "" -1 -1 "C ourier" 1 10 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Maple Plot" -1 13 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Line Printe d Output" -1 6 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Help" -1 10 1 {CSTYLE "" -1 -1 "Courier" 1 9 0 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Text Output" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 } {PSTYLE "Diagnostic" -1 9 1 {CSTYLE "" -1 -1 "Courier" 1 10 64 128 64 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Maple Output" -1 11 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Dash Item" -1 16 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 1 }{PSTYLE "Hyperlin kError" -1 206 1 {CSTYLE "" -1 -1 "Courier New" 1 12 255 0 255 1 2 2 1 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Error" -1 8 1 {CSTYLE "" -1 -1 "Courier" 1 10 255 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Title" -1 18 1 {CSTYLE "" -1 -1 "Tim es" 1 18 0 0 0 1 2 1 1 2 2 2 1 0 0 1 }3 1 0 0 12 12 2 0 2 0 2 2 -1 1 } {PSTYLE "Heading 1" -1 3 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 1 0 0 1 }1 1 0 0 8 4 2 0 2 0 2 2 -1 1 }{PSTYLE "Bullet Item" -1 15 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 1 }{PSTYLE "Heading 4" -1 20 1 {CSTYLE "" -1 -1 "Times" 1 10 0 0 0 1 1 1 2 2 2 2 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Heading 3" -1 5 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 1 1 2 2 2 2 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Heading 2" -1 4 1 {CSTYLE "" -1 -1 "Times" 1 14 0 0 0 1 2 1 2 2 2 2 1 0 0 1 }1 1 0 0 8 2 2 0 2 0 2 2 -1 1 }{PSTYLE "HyperlinkWarning" -1 207 1 {CSTYLE "" -1 -1 "Courier New" 1 12 0 0 255 1 2 2 1 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Left Justified Maple Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "List Item" -1 14 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 1 }{CSTYLE "An notation Text" -1 200 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 } {CSTYLE "Help Maple Name" -1 35 "Times" 1 12 104 64 92 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Bold" -1 5 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Menus" -1 36 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Italic" -1 3 "Times" 1 12 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Normal" -1 30 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Plot Text" -1 28 "Times" 1 8 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Nonterminal" -1 24 "Courier" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Heading" -1 26 "Times" 1 14 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Italic" -1 42 "Times" 1 12 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Variable" -1 25 "Courier" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Italic Bold" -1 40 "Times" 1 12 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "Default" -1 38 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Underlined Italic" -1 43 "Times" 1 12 0 0 0 1 1 2 1 2 2 2 0 0 0 1 }{CSTYLE "Maple Input" -1 0 "Courier" 1 12 255 0 0 1 2 1 2 2 1 2 0 0 0 1 }{CSTYLE "2D Output" -1 20 "Times" 1 12 0 0 255 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Bold Small" -1 10 " Times" 1 1 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Small" -1 7 " Times" 1 1 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Page Number" -1 33 "T imes" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Output Labels" -1 29 "Times" 1 8 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Plot Title" -1 27 "T imes" 1 10 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Emphasized" -1 22 "Times" 1 12 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Symbol 2 " -1 16 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Maple Comme nt" -1 21 "Courier" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Caption Reference" -1 201 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE " Maple Input Placeholder" -1 202 "Courier New" 1 12 200 0 200 1 2 1 2 2 1 2 0 0 0 1 }{CSTYLE "Code" -1 203 "Courier New" 1 12 255 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Italic Small205" -1 204 "Times" 1 1 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Inert Output" -1 205 "Times" 1 12 144 144 144 1 2 2 2 2 1 2 0 0 0 1 }{CSTYLE "2D Math Italic Small" -1 206 "Times" 1 1 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Comment" -1 18 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Underlin ed Bold" -1 41 "Times" 1 12 0 0 0 1 1 1 2 2 2 2 0 0 0 1 }{CSTYLE "Copy right" -1 34 "Times" 1 10 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Hyperl ink" -1 17 "Times" 1 12 0 128 128 1 2 2 1 2 2 2 0 0 0 1 }{CSTYLE "Help Underlined" -1 44 "Times" 1 12 0 0 0 1 2 2 1 2 2 2 0 0 0 1 }{CSTYLE " Prompt" -1 1 "Courier" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D M ath" -1 2 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Input" -1 19 "Times" 1 12 255 0 0 1 2 2 2 2 1 2 0 0 0 1 }{CSTYLE "Header and Footer" -1 207 "Times" 1 10 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Tex t" -1 208 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Note s" -1 37 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Bold" -1 39 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Equation Lab el" -1 209 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "LaTeX" -1 32 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Popup" -1 31 "Times" 1 12 0 128 128 1 1 2 1 2 2 2 0 0 0 1 }{CSTYLE "Atomic Variable " -1 210 "Times" 1 12 175 0 175 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "Dictio nary Hyperlink" -1 45 "Times" 1 12 147 0 15 1 2 2 1 2 2 2 0 0 0 1 } {CSTYLE "Caption Text" -1 211 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Fixed" -1 23 "Courier" 1 10 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 91 "read(\"/afs/eos.ncs u.edu/users/k/kaltofen/www/courses/LinAlgebra/Maple/LsqPkg/InitPkg.mpl \");" }}{PARA 11 "" 1 "" {XPPMATH 20 "6%Qjn/afs/eos.ncsu.edu/users/k/k altofen/www/courses/LinAlgebra/Maple6\"Q7/usr/local/maple18/libF$Q\".F $" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "with(LsqPkg);" }} {PARA 11 "" 1 "" {XPPMATH 20 "7#I-Gram_SchmidtG6\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "A := <<1,1,1,1>|<1,2,3,4>|<1,4,9,16>>;" }} {PARA 11 "" 1 "" {XPPMATH 20 "-I'MatrixG6$%*protectedGI(_syslibG6\"6#/ I$%idGF'\"59-V=eH%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "w ith(LinearAlgebra);" }}{PARA 11 "" 1 "" {XPPMATH 20 "7^sI#&xG6\"I$AddG F$I(AdjointGF$I3BackwardSubstituteGF$I+BandMatrixGF$I&BasisGF$I-Bezout MatrixGF$I/BidiagonalFormGF$I-BilinearFormGF$I%CAREGF$I5Characteristic MatrixGF$I9CharacteristicPolynomialGF$I'ColumnGF$I0ColumnDimensionGF$I 0ColumnOperationGF$I,ColumnSpaceGF$I0CompanionMatrixGF$I5CompressedSpa rseFormGF$I0ConditionNumberGF$I/ConstantMatrixGF$I/ConstantVectorGF$I% CopyGF$I2CreatePermutationGF$I-CrossProductGF$I%DAREGF$I-DeleteColumnG F$I*DeleteRowGF$I,DeterminantGF$I)DiagonalGF$I/DiagonalMatrixGF$I*Dime nsionGF$I+DimensionsGF$I+DotProductGF$I6EigenConditionNumbersGF$I,Eige nvaluesGF$I-EigenvectorsGF$I&EqualGF$I2ForwardSubstituteGF$I.Frobenius FormGF$I9FromCompressedSparseFormGF$I.FromSplitFormGF$I4GaussianElimin ationGF$I2GenerateEquationsGF$I/GenerateMatrixGF$I(GenericGF$I2GetResu ltDataTypeGF$I/GetResultShapeGF$I5GivensRotationMatrixGF$I,GramSchmidt GF$I-HankelMatrixGF$I,HermiteFormGF$I3HermitianTransposeGF$I/Hessenber gFormGF$I.HilbertMatrixGF$I2HouseholderMatrixGF$I/IdentityMatrixGF$I2I ntersectionBasisGF$I+IsDefiniteGF$I-IsOrthogonalGF$I*IsSimilarGF$I*IsU nitaryGF$I2JordanBlockMatrixGF$I+JordanFormGF$I1KroneckerProductGF$I(L A_MainGF$I0LUDecompositionGF$I-LeastSquaresGF$I,LinearSolveGF$I.Lyapun ovSolveGF$I$MapGF$I%Map2GF$I*MatrixAddGF$I2MatrixExponentialGF$I/Matri xFunctionGF$I.MatrixInverseGF$I5MatrixMatrixMultiplyGF$I+MatrixNormGF$ I,MatrixPowerGF$I5MatrixScalarMultiplyGF$I5MatrixVectorMultiplyGF$I2Mi nimalPolynomialGF$I&MinorGF$I(ModularGF$I)MultiplyGF$I,NoUserValueGF$I %NormGF$I*NormalizeGF$I*NullSpaceGF$I3OuterProductMatrixGF$I*Permanent GF$I&PivotGF$I*PopovFormGF$I1ProjectionMatrixGF$I0QRDecompositionGF$I- RandomMatrixGF$I-RandomVectorGF$I%RankGF$I6RationalCanonicalFormGF$I6R educedRowEchelonFormGF$I$RowGF$I-RowDimensionGF$I-RowOperationGF$I)Row SpaceGF$I-ScalarMatrixGF$I/ScalarMultiplyGF$I-ScalarVectorGF$I*SchurFo rmGF$I/SingularValuesGF$I*SmithFormGF$I*SplitFormGF$I8StronglyConnecte dBlocksGF$I*SubMatrixGF$I*SubVectorGF$I)SumBasisGF$I0SylvesterMatrixGF $I/SylvesterSolveGF$I/ToeplitzMatrixGF$I&TraceGF$I*TransposeGF$I0Tridi agonalFormGF$I+UnitVectorGF$I2VandermondeMatrixGF$I*VectorAddGF$I,Vect orAngleGF$I5VectorMatrixMultiplyGF$I+VectorNormGF$I5VectorScalarMultip lyGF$I+ZeroMatrixGF$I+ZeroVectorGF$I$ZipGF$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "#QRDecomposition" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "Qnorm, Rnorm := QRDecomposition(A); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6$-I'MatrixG6$%*protectedGI(_syslibG6\"6#/I$%idGF(\"5 %pI%=eH%)oW=-F$6#/F+\"5u0V=eH%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "Qnorm . Rnorm;" }}{PARA 11 "" 1 "" {XPPMATH 20 "-I'Ma trixG6$%*protectedGI(_syslibG6\"6#/I$%idGF'\"5a5V=eH%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "GS := GramSchmidt([Column(A,1), Col umn(A,2), Column(A,3)]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "7%-&I'Vector G6$%*protectedGI(_syslibG6\"6#I'columnGF)6#/I$%idGF)\"5u " 0 "" {MPLTEXT 1 0 15 "#printlevel:=3:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "Q := Matrix(GS);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 11 "" 1 "" {XPPMATH 20 "-I'MatrixG6$%*protectedGI(_syslibG6 \"6#/I$%idGF'\"55&\\^o&H%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "# Reversely engineer R from A and Q\n" }{MPLTEXT 1 0 58 "Diag := Transpose(Q) . Q; R:=Diag^(-1) . Transpose(Q) . A;" }}{PARA 11 "" 1 " " {XPPMATH 20 "-I'MatrixG6$%*protectedGI(_syslibG6\"6#/I$%idGF'\"5q) \\^o&H%)oW=" }}{PARA 11 "" 1 "" {XPPMATH 20 "-I'MatrixG6$%*protectedGI (_syslibG6\"6#/I$%idGF'\"5I3:&o&H%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 4 "Q.R;" }}{PARA 11 "" 1 "" {XPPMATH 20 "-I'MatrixG6$%*pr otectedGI(_syslibG6\"6#/I$%idGF'\"5]4:&o&H%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "b := <5,7,29,63>;" }}{PARA 11 "" 1 "" {XPPMATH 20 "-&I'VectorG6$%*protectedGI(_syslibG6\"6#I'columnGF(6#/I$% idGF(\"5q5:&o&H%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "unp rotect(D):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "D := Transpos e(Q) . Q;" }}{PARA 11 "" 1 "" {XPPMATH 20 "-I'MatrixG6$%*protectedGI(_ syslibG6\"6#/I$%idGF'\"558:&o&H%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "D^(-1);" }}{PARA 11 "" 1 "" {XPPMATH 20 "-I'MatrixG6$% *protectedGI(_syslibG6\"6#/I$%idGF'\"5I9:&o&H%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "btilde:=D^(-1) . Transpose(Q) . b;" }} {PARA 11 "" 1 "" {XPPMATH 20 "-&I'VectorG6$%*protectedGI(_syslibG6\"6# I'columnGF(6#/I$%idGF(\"5!R_^o&H%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 2 "R;" }}{PARA 11 "" 1 "" {XPPMATH 20 "-I'MatrixG6$%*prot ectedGI(_syslibG6\"6#/I$%idGF'\"5I3:&o&H%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "xhat:=R^(-1) . btilde;" }}{PARA 11 "" 1 "" {XPPMATH 20 "-&I'VectorG6$%*protectedGI(_syslibG6\"6#I'columnGF(6#/I$% idGF(\"59[9&o&H%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "Lin earSolve(Transpose(A) . A, Transpose(A) . b);" }}{PARA 11 "" 1 "" {XPPMATH 20 "-&I'VectorG6$%*protectedGI(_syslibG6\"6#I'columnGF(6#/I$% idGF(\"59g9&o&H%)oW=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }} }{EXCHG {PARA 0 "> " 0 "" {XPPEDIT 19 1 "" "%#%?G" }}}} {MARK "0 0 0" 0 }{VIEWOPTS 1 1 0 2 1 1805 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }