sys = BioSystem();

sys.AddConstant('k_tc',0.01)
sys.AddConstant('k_mdeg',0.003);
sys.AddConstant('k_tl',0.1);
sys.AddConstant('k_pdeg',0.03);

dDNAdt = sys.AddCompositor('DNA',10);
dmRNAdt = sys.AddCompositor('mRNA',0);
dGFPdt = sys.AddCompositor('GFP',0);

sys.AddPart(Part('translation',[dDNAdt dmRNAdt dGFPdt],[...
    Rate('0') ...
    Rate('k_tc*DNA-k_mdeg*mRNA') ...
    Rate('k_tl*mRNA-k_pdeg*GFP')]));

[T,Y] = sys.run_pulses([...
    Pulse(0,'DNA',10), ...
    Pulse(12000,'DNA',0)]);

figure();
plot(T,Y(:,sys.CompositorIndex('DNA')),...
    T,Y(:,sys.CompositorIndex('mRNA')),...
    T,Y(:,sys.CompositorIndex('GFP')));
legend('DNA','mRNA','GFP');
xlabel('time');
ylabel('concentration');