with(LinearAlgebra):A := Matrix([[3,-1,1],[3,6,2],[3,3,7]] );

NiM+SSJBRzYiLUknUlRBQkxFR0YlNiUiKDcsMSItSSdNQVRSSVhHRiU2IzclNyUiIiQhIiIiIiI3JUYvIiInIiIjNyVGL0YvIiIoSSdNYXRyaXhHNiRJKnByb3RlY3RlZEdGOUkoX3N5c2xpYkdGJQ==

b := Vector([1,0,4]);

NiM+SSJiRzYiLUknUlRBQkxFR0YlNiUiJzdjYS1JJ01BVFJJWEdGJTYjNyU3IyIiIjcjIiIhNyMiIiUmSSdWZWN0b3JHNiRJKnByb3RlY3RlZEdGN0koX3N5c2xpYkdGJTYjSSdjb2x1bW5HRiU=

z := evalf(LinearAlgebra[LinearSolve](A,b));

NiM+SSJ6RzYiLUknUlRBQkxFR0YlNiUiKTsjMyskLUknTUFUUklYR0YlNiM3JTcjJCIrST54M04hIzY3IyQhK2A1VW9CISM1NyMkIitvdCUqeWxGNSZJJ1ZlY3Rvckc2JEkqcHJvdGVjdGVkR0Y8SShfc3lzbGliR0YlNiNJJ2NvbHVtbkdGJQ==

E := LinearAlgebra[DiagonalMatrix]([3,6,7]);

NiM+SSJFRzYiLUknUlRBQkxFR0YlNiUiJ1tzJyotSSdNQVRSSVhHRiU2IzclNyUiIiQiIiFGMDclRjAiIidGMDclRjBGMCIiKEknTWF0cml4RzYkSSpwcm90ZWN0ZWRHRjdJKF9zeXNsaWJHRiU=

B := A-E;

NiM+SSJCRzYiLUknUlRBQkxFR0YlNiUiKT9wdUgtSSdNQVRSSVhHRiU2IzclNyUiIiEhIiIiIiI3JSIiJEYvIiIjNyVGM0YzRi9JJ01hdHJpeEc2JEkqcHJvdGVjdGVkR0Y4SShfc3lzbGliR0Yl

C := -1.0*MatrixInverse(E).B;

NiM+SSJDRzYiLUknUlRBQkxFR0YlNiUiKGc5LCItSSdNQVRSSVhHRiU2IzclNyUkIiIhRjAkIjM6TExMTExMTEwhIz0kITM6TExMTExMTExGMzclJCEzKysrKysrKytdRjNGL0Y0NyUkITNbJkc5ZEc5ZEclRjNGOkYvSSdNYXRyaXhHNiRJKnByb3RlY3RlZEdGPkkoX3N5c2xpYkdGJQ==

d := 1.0*MatrixInverse(E).b;

NiM+SSJkRzYiLUknUlRBQkxFR0YlNiUiKUsheTAkLUknTUFUUklYR0YlNiM3JTcjJCIzOkxMTExMTExMISM9NyMkIiIhRjQ3IyQiMyhSciZHOWRHOWRGMSZJJ1ZlY3Rvckc2JEkqcHJvdGVjdGVkR0Y7SShfc3lzbGliR0YlNiNJJ2NvbHVtbkdGJQ==

x[0] := Vector([0,0,0]);

NiM+JkkieEc2IjYjIiIhLUknUlRBQkxFR0YmNiUiKTcmb0skLUknTUFUUklYR0YmNiM3JTcjRihGMUYxJkknVmVjdG9yRzYkSSpwcm90ZWN0ZWRHRjVJKF9zeXNsaWJHRiY2I0knY29sdW1uR0Ym

n := 10;

NiM+SSJuRzYiIiM1

for k from 0 to n-1 do x[k+1] := C.x[k]+d: end do: x[n],z;

NiQtSSdSVEFCTEVHNiI2JSIpMzckKUgtSSdNQVRSSVhHRiU2IzclNyMkIjNKcWxsJioqUXldJCEjPjcjJCEzI0cxW2hwaCNwQiEjPTcjJCIzST9AdTFYLHlsRjMmSSdWZWN0b3JHNiRJKnByb3RlY3RlZEdGOkkoX3N5c2xpYkdGJTYjSSdjb2x1bW5HRiUtRiQ2JSIpOyMzKyQtRik2IzclNyMkIitJPngzTiEjNjcjJCErYDVVb0IhIzU3IyQiK290JSp5bEZLRjc=

Example 2, p 310A := Matrix([[10,-1,2,0],[-1,11,-1,3],[2,-1,10,-1],[0,3,-1,8]]);

NiM+SSJBRzYiLUknUlRBQkxFR0YlNiUiKUtnJDMkLUknTUFUUklYR0YlNiM3JjcmIiM1ISIiIiIjIiIhNyZGMCIjNkYwIiIkNyZGMUYwRi9GMDcmRjJGNUYwIiIpSSdNYXRyaXhHNiRJKnByb3RlY3RlZEdGO0koX3N5c2xpYkdGJQ==

b := Vector([6,25,-11,15]);

NiM+SSJiRzYiLUknUlRBQkxFR0YlNiUiKC93PyItSSdNQVRSSVhHRiU2IzcmNyMiIic3IyIjRDcjISM2NyMiIzomSSdWZWN0b3JHNiRJKnByb3RlY3RlZEdGOUkoX3N5c2xpYkdGJTYjSSdjb2x1bW5HRiU=

z := LinearSolve(A,b);

NiM+SSJ6RzYiLUknUlRBQkxFR0YlNiUiKWskcChILUknTUFUUklYR0YlNiM3JjcjIiIiNyMiIiM3IyEiIkYuJkknVmVjdG9yRzYkSSpwcm90ZWN0ZWRHRjdJKF9zeXNsaWJHRiU2I0knY29sdW1uR0Yl

E := DiagonalMatrix([10,11,10,8]);

NiM+SSJFRzYiLUknUlRBQkxFR0YlNiUiKTs6RUgtSSdNQVRSSVhHRiU2IzcmNyYiIzUiIiFGMEYwNyZGMCIjNkYwRjA3JkYwRjBGL0YwNyZGMEYwRjAiIilJJ01hdHJpeEc2JEkqcHJvdGVjdGVkR0Y4SShfc3lzbGliR0Yl

B := A-E;

NiM+SSJCRzYiLUknUlRBQkxFR0YlNiUiKWc+bkstSSdNQVRSSVhHRiU2IzcmNyYiIiEhIiIiIiNGLzcmRjBGL0YwIiIkNyZGMUYwRi9GMDcmRi9GM0YwRi9JJ01hdHJpeEc2JEkqcHJvdGVjdGVkR0Y4SShfc3lzbGliR0Yl

C := -1.0*MatrixInverse(E).B;

NiM+SSJDRzYiLUknUlRBQkxFR0YlNiUiKSVbOzMkLUknTUFUUklYR0YlNiM3JjcmJCIiIUYwJCIzLysrKysrKys1ISM9JCEzNSsrKysrKys/RjNGLzcmJCIzOyI0NDQ0NDQ0KiEjPkYvRjckITMyRkZGRkZGRkZGMzcmRjRGMUYvRjE3JkYvJCEzKysrKysrK11QRjMkIjMrKysrKysrXTdGM0YvSSdNYXRyaXhHNiRJKnByb3RlY3RlZEdGREkoX3N5c2xpYkdGJQ==

d := 1.0*MatrixInverse(E).b;

NiM+SSJkRzYiLUknUlRBQkxFR0YlNiUiKSM0aTkkLUknTUFUUklYR0YlNiM3JjcjJCIzKTMrKysrKysrJyEjPTcjJCIzJUhGRkZGRkZGIyEjPDcjJCEzMysrKysrKys2RjU3IyQiMysrKysrKyt2PUY1JkknVmVjdG9yRzYkSSpwcm90ZWN0ZWRHRj9JKF9zeXNsaWJHRiU2I0knY29sdW1uR0Yl

x[0] := ZeroVector(4);

NiM+JkkieEc2IjYjIiIhLUknUlRBQkxFR0YmNiUiKTt1XUwtSSdNQVRSSVhHRiY2IzcmNyNGKEYxRjFGMSZJJ1ZlY3Rvckc2JEkqcHJvdGVjdGVkR0Y1SShfc3lzbGliR0YmNiNJJ2NvbHVtbkdGJg==

n := 30: for k from 0 to n-1 do x[k+1] := C.x[k]+d: end do: x[n];

NiMtSSdSVEFCTEVHNiI2JSIpL1V0Si1JJ01BVFJJWEdGJTYjNyY3IyQiM1NXXSsrKysrNSEjPDcjJCIzKTRhIioqKioqKioqKio+Ri83IyQhMi1iTSoqKioqKioqKioqRi83IyQiMm9rMSoqKioqKioqKioqRi8mSSdWZWN0b3JHNiRJKnByb3RlY3RlZEdGPEkoX3N5c2xpYkdGJTYjSSdjb2x1bW5HRiU=

TTdSMApJNFJUQUJMRV9TQVZFLzEwNjAxMTJYLCUpYW55dGhpbmdHNiI2IltnbCEiJSEhISMqIiQiJCIiJEYnRichIiIiIidGJyIiIiIiIyIiKEYmCg==TTdSMApJM1JUQUJMRV9TQVZFLzU0NTYxMlgqJSlhbnl0aGluZ0c2IjYiW2dsISMlISEhIiQiJCIiIiIiISIiJUYmCg==TTdSMApJNVJUQUJMRV9TQVZFLzMwMDA4MjE2WColKWFueXRoaW5nRzYiNiJbZ2whIyUhISEiJCIkJCIrST54M04hIzYkIStgNVVvQiEjNSQiK290CiUqeWxGLEYmCg==TTdSMApJM1JUQUJMRV9TQVZFLzk2NzI0OFgsJSlhbnl0aGluZ0c2IjYiW2dsISIlISEhIyoiJCIkIiIkIiIhRihGKCIiJ0YoRihGKCIiKEYmCg==TTdSMApJNVJUQUJMRV9TQVZFLzI5NzQ2OTIwWCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjKiIkIiQiIiEiIiRGKCEiIkYnRigiIiIiIiNGJ0YmCg==TTdSMApJNFJUQUJMRV9TQVZFLzEwMTE0NjBYLCUpYW55dGhpbmdHNiI2IltnbCciJSEhISMqIiQiJDgwMDAwMDAwMDAwMDAwMDBCRkUwMDAwMDAwCjAwMDAwMEJGREI2REI2REI2REI2REIzRkQ1NTU1NTU1NTU1NTU1ODAwMDAwMDAwMDAwMDAwMEJGREI2REI2REI2REI2REJCRkQ1NTU1NTUKNTU1NTU1NUJGRDU1NTU1NTU1NTU1NTU4MDAwMDAwMDAwMDAwMDAwRiYKTTdSMApJNVJUQUJMRV9TQVZFLzMwNTc4MDMyWColKWFueXRoaW5nRzYiNiJbZ2wnIyUhISEiJCIkM0ZENTU1NTU1NTU1NTU1NTAwMDAwMDAwMDAwCjAwMDAwM0ZFMjQ5MjQ5MjQ5MjQ5MkYmCg==TTdSMApJNVJUQUJMRV9TQVZFLzMzMjY4NTEyWColKWFueXRoaW5nRzYiNiJbZ2whIyUhISEiJCIkIiIhRidGJ0YmCg==TTdSMApJNVJUQUJMRV9TQVZFLzI5ODMxMjA4WColKWFueXRoaW5nRzYiNiJbZ2wnIyUhISEiJCIkM0ZBMUY1Q0I3MzUxQ0Y0OEJGQ0U1Mzk4QzMwCjdGNDQ0M0ZFNTBDQjVBMEJEMDZGMkYmCg==TTdSMApJNVJUQUJMRV9TQVZFLzMwMDA4MjE2WColKWFueXRoaW5nRzYiNiJbZ2whIyUhISEiJCIkJCIrST54M04hIzYkIStgNVVvQiEjNSQiK290CiUqeWxGLEYmCg==TTdSMApJNVJUQUJMRV9TQVZFLzMwODM2MDMyWCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjMSIlIiUiIzUhIiIiIiMiIiFGKCIjNkYoIiIkRilGCihGJ0YoRipGLEYoIiIpRiYKTTdSMApJNFJUQUJMRV9TQVZFLzEyMDc2MDRYKiUpYW55dGhpbmdHNiI2IltnbCEjJSEhISIlIiUiIiciI0QhIzYiIzpGJgo=TTdSMApJNVJUQUJMRV9TQVZFLzI5NzY5MzY0WColKWFueXRoaW5nRzYiNiJbZ2whIyUhISEiJSIlIiIiIiIjISIiRidGJgo=TTdSMApJNVJUQUJMRV9TQVZFLzI5MjYxNTE2WCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjMSIlIiUiIzUiIiFGKEYoRigiIzZGKEYoRihGKEYnCkYoRihGKEYoIiIpRiYKTTdSMApJNVJUQUJMRV9TQVZFLzMyNjcxOTYwWCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjMSIlIiUiIiEhIiIiIiNGJ0YoRidGKCIiJEYpRihGCidGKEYnRipGKEYnRiYKTTdSMApJNVJUQUJMRV9TQVZFLzMwODE2NDg0WCwlKWFueXRoaW5nRzYiNiJbZ2wnIiUhISEjMSIlIiU4MDAwMDAwMDAwMDAwMDAwM0ZCNzQ1RDE3CjQ1RDE3NDZCRkM5OTk5OTk5OTk5OTlBODAwMDAwMDAwMDAwMDAwMDNGQjk5OTk5OTk5OTk5OUE4MDAwMDAwMDAwMDAwMDAwM0ZCOTk5OTkKOTk5OTk5OUFCRkQ4MDAwMDAwMDAwMDAwQkZDOTk5OTk5OTk5OTk5QTNGQjc0NUQxNzQ1RDE3NDY4MDAwMDAwMDAwMDAwMDAwM0ZDMDAwMAowMDAwMDAwMDA4MDAwMDAwMDAwMDAwMDAwQkZEMTc0NUQxNzQ1RDE3NDNGQjk5OTk5OTk5OTk5OUE4MDAwMDAwMDAwMDAwMDAwRiYKTTdSMApJNVJUQUJMRV9TQVZFLzMxNDYyMDkyWColKWFueXRoaW5nRzYiNiJbZ2wnIyUhISEiJSIlM0ZFMzMzMzMzMzMzMzMzNDQwMDIyRThCQTJFCjhCQTJGQkZGMTk5OTk5OTk5OTk5QTNGRkUwMDAwMDAwMDAwMDBGJgo=TTdSMApJNVJUQUJMRV9TQVZFLzMzNTA3NDE2WColKWFueXRoaW5nRzYjJSV6ZXJvRzYiW2dsISMiISEhIiEiJUYnCg==TTdSMApJNVJUQUJMRV9TQVZFLzMxNzM0MjA0WColKWFueXRoaW5nRzYiNiJbZ2wnIyUhISEiJSIlM0ZGMDAwMDAwMDAwNThCRTNGRkZGRkZGRkZGCkY2QjMwQkZFRkZGRkZGRkZGMTlCODNGRUZGRkZGRkZGRUI3OEJGJgo=