follicle_model_BT.m

 load('follicle_model_system_matrices.mat');
 sys = ss(A,B,C,D);
 tic;
 [sysb,G,T,Tinv] = balreal(sys);
 time = toc;
 W = transpose(T(1:500,:));
 V = Tinv(:,1:500);
 H = G(1:500);
 save('~/follicle_model_BT_basis.mat','V','W','H','time');