octave> P=driftmatrix(3,0,0) P = 1.00000 0.33490 0.08779 0.01562 0.00137 0.00002 0.00000 0.00000 0.40188 0.26337 0.09375 0.01646 0.00064 0.00000 0.00000 0.20094 0.32922 0.23438 0.08230 0.00804 0.00000 0.00000 0.05358 0.21948 0.31250 0.21948 0.05358 0.00000 0.00000 0.00804 0.08230 0.23438 0.32922 0.20094 0.00000 0.00000 0.00064 0.01646 0.09375 0.26337 0.40188 0.00000 0.00000 0.00002 0.00137 0.01562 0.08779 0.33490 1.00000 octave> [X,L]=eig(P) X = 1.00000 0.00000 -0.59720 -0.03290 0.42135 -0.10911 -0.24305 0.00000 0.00000 0.21852 0.19739 -0.48154 0.43644 0.57290 0.00000 0.00000 0.25059 -0.49346 -0.30096 -0.54554 -0.10416 0.00000 0.00000 0.25616 0.65795 0.00000 -0.00000 -0.45137 0.00000 0.00000 0.25059 -0.49346 0.30096 0.54554 -0.10416 0.00000 0.00000 0.21852 0.19739 0.48154 -0.43644 0.57290 0.00000 1.00000 -0.59720 -0.03290 -0.42135 0.10911 -0.24305 L = 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.83333 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.01543 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.55556 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.09259 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.27778 octave> plot(1:7,X(:,1));
octave> plot(1:7,X(:,2));
octave> P=driftmatrix(3,0.1,0.1) P = Columns 1 through 6: 5.3144e-01 2.0307e-01 6.4535e-02 1.5625e-02 2.4301e-03 1.6138e-04 3.5429e-01 3.7082e-01 2.2417e-01 9.3750e-02 2.5185e-02 3.1816e-03 9.8415e-02 2.8214e-01 3.2446e-01 2.3438e-01 1.0875e-01 2.6134e-02 1.4580e-02 1.1449e-01 2.5046e-01 3.1250e-01 2.5046e-01 1.1449e-01 1.2150e-03 2.6134e-02 1.0875e-01 2.3438e-01 3.2446e-01 2.8214e-01 5.4000e-05 3.1816e-03 2.5185e-02 9.3750e-02 2.2417e-01 3.7082e-01 1.0000e-06 1.6138e-04 2.4301e-03 1.5625e-02 6.4535e-02 2.0307e-01 Column 7: 1.0000e-06 5.4000e-05 1.2150e-03 1.4580e-02 9.8415e-02 3.5429e-01 5.3144e-01 octave> [X,L]=eig(P) X = Columns 1 through 5: -2.4287e-01 4.5634e-01 4.3095e-01 3.6663e-01 -2.3072e-01 -3.8043e-01 1.5873e-01 4.8360e-01 -2.8270e-01 5.1486e-01 -4.3834e-01 -3.4088e-01 2.8358e-01 -5.3448e-01 1.7028e-02 -4.5639e-01 -5.4838e-01 5.3338e-16 1.7542e-16 -6.0233e-01 -4.3834e-01 -3.4088e-01 -2.8358e-01 5.3448e-01 1.7028e-02 -3.8043e-01 1.5873e-01 -4.8360e-01 2.8270e-01 5.1486e-01 -2.4287e-01 4.5634e-01 -4.3095e-01 -3.6663e-01 -2.3072e-01 Columns 6 and 7: -3.2898e-02 -1.0911e-01 1.9739e-01 4.3644e-01 -4.9346e-01 -5.4554e-01 6.5795e-01 -2.0070e-15 -4.9346e-01 5.4554e-01 1.9739e-01 -4.3644e-01 -3.2898e-02 1.0911e-01 L = 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.53333 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.80000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.28444 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.11378 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00405 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.03034 octave> x=X(:,1)/sum(X(:,1)); octave> plot(1:7,x);
octave> P=driftmatrix(3,0.1,0.0) P = 1.00000 0.37715 0.11765 0.02768 0.00410 0.00024 0.00000 0.00000 0.39933 0.30253 0.13589 0.03686 0.00439 0.00005 0.00000 0.17618 0.32413 0.27795 0.13824 0.03296 0.00121 0.00000 0.04145 0.18522 0.30322 0.27648 0.13184 0.01458 0.00000 0.00549 0.05953 0.18607 0.31104 0.29663 0.09841 0.00000 0.00039 0.01021 0.06089 0.18662 0.35596 0.35429 0.00000 0.00001 0.00073 0.00830 0.04666 0.17798 0.53144 octave> [X,L]=eig(P) X = 1.00000 -0.92303 0.72750 -0.03290 0.26499 0.11360 0.48372 0.00000 0.17240 -0.49690 0.19739 -0.71446 -0.47833 -0.73309 0.00000 0.18451 -0.37505 -0.49346 0.36353 0.68760 -0.16798 0.00000 0.17635 -0.17937 0.65795 0.39075 -0.23917 0.25897 0.00000 0.16065 0.00550 -0.49346 -0.18363 -0.32885 0.30604 0.00000 0.13786 0.14683 0.19739 -0.27673 0.33483 0.04604 0.00000 0.09125 0.17150 -0.03290 0.15555 -0.08969 -0.19371 L = 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.90000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.67500 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00820 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.18225 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.05468 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.40500 octave> plot(1:7,X(:,1));
octave> P=driftmatrix(3,0.0,0.1) P = 0.53144 0.17798 0.04666 0.00830 0.00073 0.00001 0.00000 0.35429 0.35596 0.18662 0.06089 0.01021 0.00039 0.00000 0.09842 0.29663 0.31104 0.18607 0.05954 0.00549 0.00000 0.01458 0.13184 0.27648 0.30322 0.18522 0.04145 0.00000 0.00122 0.03296 0.13824 0.27795 0.32414 0.17618 0.00000 0.00005 0.00439 0.03686 0.13589 0.30253 0.39933 0.00000 0.00000 0.00024 0.00410 0.02768 0.11765 0.37715 1.00000 octave> [X,L]=eig(P) X = 0.00000 -0.09125 -0.17150 -0.15555 -0.08969 0.03290 -0.19371 0.00000 -0.13786 -0.14683 0.27673 0.33483 -0.19739 0.04604 0.00000 -0.16065 -0.00550 0.18363 -0.32885 0.49346 0.30604 0.00000 -0.17635 0.17937 -0.39075 -0.23917 -0.65795 0.25897 0.00000 -0.18451 0.37505 -0.36353 0.68760 0.49346 -0.16798 0.00000 -0.17240 0.49690 0.71446 -0.47833 -0.19739 -0.73309 1.00000 0.92303 -0.72750 -0.26499 0.11360 0.03290 0.48372 L = 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.90000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.67500 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.18225 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.05468 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00820 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.40500 octave> plot(1:7,X(:,1));
octave>