<v, Krein array>
<2600, [ 25, 3718/207, 56875/5819, 1, 7150/4761, 325/253 ]>

P =
[   1   69  759 1771]
[   1   43  187 -231]
[   1   19  -41   21]
[   1   -3    3   -1]

Q =
[       1       25      299     2275]
[       1  1075/69    247/3 -2275/23]
[       1   425/69  -533/33 2275/253]
[       1   -75/23    39/11 -325/253]

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

    [ 0 69  0  0]
    [ 1 24 44  0]
    [ 0  4 44 21]
    [ 0  0  9 60],

    [  0   0 759   0]
    [  0  44 484 231]
    [  1  44 294 420]
    [  0   9 180 570],

    [   0    0    0 1771]
    [   0    0  231 1540]
    [   0   21  420 1330]
    [   1   60  570 1140]
]

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

    [         0         25          0          0]
    [         1   1250/207   3718/207          0]
    [         0  7150/4761 31250/2277 56875/5819]
    [         0          0    325/253   6000/253],

    [          0           0         299           0]
    [          0    3718/207    16250/99   29575/253]
    [          1  31250/2277  87022/1089 568750/2783]
    [          0     325/253    3250/121 753792/2783],

    [              0               0               0            2275]
    [              0               0       29575/253      546000/253]
    [              0      56875/5819     568750/2783 131913600/64009]
    [              1        6000/253     753792/2783    5508750/2783]
]