clear;clc;
G = rand(3,3)
cG = poly(G);
GJeff = companjeff(cG)
% G = [1 2; 3 4];
% GJeff = [5 6; 7 8];
T = transitionjeff(G,GJeff)
Tinv = inv(T)
Check = T*G - GJeff*T % Should be 0
G
inv(T)*GJeff*T
%T*GJeff*inv(T)