En este artículo se presenta una extensión del algoritmo del método de Gauthier, que soluciona la búsqueda de realización mínima multivariable partiendo de matrices de transferencia cuadradas. El algoritmo incorpora previamente la teoría de fracciones coprimas, desarrolladas con matrices de Silvester y factorización qr. Debido a que las fracciones coprimas tienen una especial relación con las matrices en fracción polinomial, se muestran sus diferencias, analizándolas independientemente. Se plantean las características generales y se nombran las funciones desarrolladas para hacer hincapié en los caminos de búsqueda de la fracción coprima que no son únicos, así como tampoco su representación en espacio de estado. Para la demostración se utilizó un sistema dinámico multivariable, donde se comprueban la eficiencia y las limitaciones del algoritmo elaborado, con base en funciones realizadas con la Toolbox de control de Matlab®.