I have the following functions:
P[t_] := P[t] = P[t-1] +a*ED[t-1]; ED[t_] := ED[t] = DF[t] + DC[t]; DF[t_] := DF[t] = b (F - P[t]); DC[t_] := DC[t] = c (P[t]