dtgsyl example program

This example solves the generalized Sylvester equations

\begin{matrix} A R - L B & = α C \\ D R - L E & = α F, \end{matrix}

where

A = (\begin{array}{r} 4.0 & 1.0 & 1.0 & 2.0 \\ 0 & 3.0 & 4.0 & 1.0 \\ 0 & 1.0 & 3.0 & 1.0 \\ 0 & 0 & 0 & 6.0 \end{array}), B = (\begin{array}{r} 1.0 & 1.0 & 1.0 & 1.0 \\ 0 & 3.0 & 4.0 & 1.0 \\ 0 & 1.0 & 3.0 & 1.0 \\ 0 & 0 & 0 & 4.0 \end{array}),

D = (\begin{array}{r} 2.0 & 1.0 & 1.0 & 3.0 \\ 0 & 1.0 & 2.0 & 1.0 \\ 0 & 0 & 1.0 & 1.0 \\ 0 & 0 & 0 & 2.0 \end{array}), E = (\begin{array}{r} 1.0 & 1.0 & 1.0 & 2.0 \\ 0 & 1.0 & 4.0 & 1.0 \\ 0 & 0 & 1.0 & 1.0 \\ 0 & 0 & 0 & 1.0 \end{array}),

C = (\begin{array}{r} - 4.0 & 7.0 & 1.0 & 12.0 \\ - 9.0 & 2.0 & - 2.0 & - 2.0 \\ - 4.0 & 2.0 & - 2.0 & 8.0 \\ - 7.0 & 7.0 & - 6.0 & 19.0 \end{array}) and F = (\begin{array}{r} - 7.0 & 5.0 & 0.0 & 7.0 \\ - 5.0 & 1.0 & - 8.0 & 0.0 \\ - 1.0 & 2.0 & - 3.0 & 5.0 \\ - 3.0 & 2.0 & 0.0 & 5.0 \end{array}) .

Example program for dtgsyl