Advanced Search

Question: Converting the very old version codes to Modern version codes

Question:Converting the very old version codes to Modern version codes

q41b3 20 Maple
linearalgebra linalg
June 27 2019
0

Could you help me converting this old version code to modern version code(Maple 2017)?
 

> 

restart
> 

N := 2; A := -N; B := N
> 

q := .3; p := .5; sa := .9; sb := .1; r := 1-p-q
> 

dimP := 2*N+1
> 

P := Matrix(dimP, dimP)
> 

P[1, 1] := sa; P[1, 2] := 1-sa; P[dimP, dimP] := sb; P[dimP, dimP-1] := 1-sb
> 

for i from 2 to dimP-1 do P[i, i-1] := q; P[i, i] := r; P[i, i+1] := p end do
> 

P

Matrix(%id = 18446745940352174910)

 (1)
> 

# change this part code to the modernversion with(linalg)
> 

J := diag(`$`(1, dimP))
> 

d := matrix(dimP, 1, [`$`(1, dimP)])
> 

b := matrix(dimP+1, 1, [`$`(0, dimP), 1])
> 

A := transpose(augment(P-J, d))
> 

linsolve(A, b)

linsolve(transpose(augment(Matrix(%id = 18446745940387027294), d)), b)

 (2)
> 

``


 

Download mdernVersion.mw
Please Wait...