function [An, Bn, Cn, Dn] = changebase(A, B, C, D, N)
    An = inv(N)*A*N;
    Bn = inv(N)*B;
    Cn = C*N;
    Dn = D;
end