close all;
plot(T,tot_S1_c);
title('STAT1 protein in cytoplasm (total)')
figure(2)
plot(T,tot_S1_pc);
title('Phosphorylated STAT1 in cytoplasm')

figure(3)
plot(T,tot_S1_n);
title('STAT1 protein in nucleus (total)')

figure(4)
plot(T,tot_S1_pn);
title('Phosphorylated STAT1 in nucleus')

figure(5)
plot(T,Y(:,9)*nuc_scale);
title('STAT1 homodimers in nucleus')

figure(6);
plot(T,Y(:,10)*nuc_scale);
title('STAT heterodimers in nucleus')

figure(7)
plot(T,(Y(:,2)+Y(:,4)+Y(:,8))*cyt_scale);
title('STAT2 protein in cytoplasm (total)')

figure(8)
plot(T, (Y(:,4)+Y(:,8))*cyt_scale);
title('Phospho STAT2 in cytoplasm')

figure(9)
plot(T, (Y(:,6)+ Y(:,10))*nuc_scale);
title('STAT2 protein in nucleus (total)');

figure(10)
plot(T,(Y(:,10))*nuc_scale);
title('Phospho STAT2 in nucleus');

figure(11)
plot(T,Y(:,15)*cyt_scale);
title('IRF1 mRNA')

figure(12)
plot(T,(Y(:,16)+Y(:,17))*cyt_scale)
title('IRF1 protein in cytoplasm');

figure(13)
plot(T,(Y(:,18)+Y(:,19))*nuc_scale)
title('IRF1 protein in nucleus')

figure(14)
plot(T,Y(:,11)*cyt_scale);
title('STAT1 mRNA');