<v, Krein array>
<2001, [ 160, 7935/58, 161/16, 1, 368/29, 2415/16 ]>

P =
[  1 400 960 640]
[  1  55  -6 -50]
[  1  -3  -6   8]
[  1 -35  90 -56]

Q =
[      1     160    1725     115]
[      1      22 -207/16 -161/16]
[      1      -1 -345/32  345/32]
[      1   -25/2  345/16 -161/16]

L = 
[
    [1 0 0 0]
    [0 1 0 0]
    [0 0 1 0]
    [0 0 0 1],

    [  0 400   0   0]
    [  1 107 204  88]
    [  0  85 175 140]
    [  0  55 210 135],

    [  0   0 960   0]
    [  0 204 420 336]
    [  1 175 504 280]
    [  0 210 420 330],

    [  0   0   0 640]
    [  0  88 336 216]
    [  0 140 280 220]
    [  1 135 330 174]
]

L* =
[
    [1 0 0 0]
    [0 1 0 0]
    [0 0 1 0]
    [0 0 0 1],

    [        0       160         0         0]
    [        1   1287/58   7935/58         0]
    [        0    368/29 63683/464    161/16]
    [        0         0   2415/16    145/16],

    [             0              0           1725              0]
    [             0        7935/58 21970635/14848      55545/512]
    [             1      63683/464 22098699/14848      50393/512]
    [             0        2415/16     755895/512      50025/512],

    [        0         0         0       115]
    [        0         0 55545/512  3335/512]
    [        0    161/16 50393/512  3335/512]
    [        1    145/16 50025/512  3703/512]
]