|
MATHEMATICA
|
N0=4 v=Table[1, {n, 1, N0}] a0=2*IdentityMatrix[N0] a0[[N0, N0]]=1 a1=Table[If[n-m+1==0, -1, 0], {n, 1, N0}, {m, 1, N0}] a2=Table[If[m-n+1==0, -1, 0], {n, 1, N0}, {m, 1, N0}] a=a0+a1+a2 (* allometric matrix N-K>exp{1]*) anm=3*MatrixPower[a, 3]-MatrixPower[a, 2] digits=25 aa=Table[MatrixPower[anm, n].v, {n, 1, digits}] aa1=Table[aa[[n, 1]], {n, 1, digits}] anm=3*MatrixPower[a, 3]-MatrixPower[a, 2]
|