4734.sagews

L=matrix([[0,1,1.5,0,0],[0.9,0,0,0,0],[0,0.8,0,0,0],[0,0,0.7,0,0],[0,0,0,0.5,0]]);L

[0.000 1.000 1.500 0.000 0.000] [0.900 0.000 0.000 0.000 0.000] [0.000 0.800 0.000 0.000 0.000] [0.000 0.000 0.700 0.000 0.000] [0.000 0.000 0.000 0.500 0.000]

L^30

[ 1411.0818246953 2057.977947814 1612.6777054060 0.00 0.00] [ 967.60662324363 1411.081824695 1105.7169122365 0.00 0.00] [ 589.71568652617 860.0947762165 673.93721653758 0.00 0.00] [ 314.50403438420 458.6677561870 359.45288620923 0.00 0.00] [ 119.81762873641 174.7244635467 136.89531169637 0.00 0.00]

X=matrix([[10000],[20000],[15000],[10000],[8000]]);X

[10000] [20000] [15000] [10000] [ 8000]

L^30*X

[7.94605427843305e7] [5.44834564098910e7] [3.32081106376568e7] [1.77101887607211e7] [6.74609523374534e6]