% .. code1equations.m
% .. rate equations for code1

function dPdt = code1equations(t,P,kProt,dProt)

dPdt = [kProt - dProt*P(1)];
