Este documento presenta un algoritmo computacional para derivar la teoría de perturbaciones gauge invariantes lineales en espacios-tiempo cosmológicos anisotrópicos del tipo Bianchi I. Nuestro código se basa en los paquetes de álgebra tensorial y , dentro de la infraestructura computacional de escrito en . El algoritmo se basa en una formulación Hamiltoniana o de espacio de fases, y proporciona una forma eficiente y transparente de aislar los grados de libertad gauge invariantes en los campos de perturbación y obtener el Hamiltoniano que genera su dinámica. La restricción a espacios-tiempo de Friedmann-Lemaître-Robertson-Walker es directa.